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The effect of the GB in n-type Ge on hole diffusion, GB-phctoelectric phenomena 
and transients was investigated. Boundaries of two kinds were found to exist. GB of the 
first kind are characterized by enhanced hole recombination and the absence of photo- 
electric phenomena. GB of the second kind revealed no recombinative action and, moreo- 
ver, the diffusion length of the holes as measured in the plane of lineage were, usually, 
considerably in excess of the vulume diffusion values. Such GB were the site of intense 
photoelectric effects. The conclusion is reached that a n-pt-n structure (corresponding 
to high potential barriers) is formed on the GB of the second kind, and that their proper- 
ties are determined by a specific “feed-in, feed-out” effect. GB of the first kind may 
be considered to correspond to n-p-n structures, or to dislocation lineage. 


I. Introduction 


Recently, interest has centered on the investigation of the electric properties of grain 
boundaries (GB) in semiconductors. This is due both to practical points of view (the 
possibilities of adapting GB effects to the construction of electronic devices) and to 
the information to be expected therefrom on the part played by dislocations in semi- 
conductors. 

Least of all is known as to the effect of GB on injected minority current carriers. 
Low-angle (dislocation) grain boundaries in n-type germanium are known to act as 
effective recombination centres for the holes (Vogel et al. 1954, McKelvey 1957). 
Some unexpected anomalies were detected by Mataré (Mataré and Reed 1956, Mataré 
1959) in measurements of diffusion length near GB of large angles of misfit. Moreover, 
photovoltaic and transistor effects on GB were observed (Pearson 1948, Mataré 1954, 
1955, 1956). No satisfactory theory of these effects has been proposed hitherto. 

It was the aim of the present author to investigate basic electronic processes occurr- 
ing on GB in the presence of injected minority current carriers. Investigation covered 
a large number of natural GB present in n-Ge. Some results of the experiments dis- 
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cussed, together with a preliminary interpretation, have already been published in the 
form of notes (Figielski and Sosnowski 1958, Sosnowski 1959). The interpretation of 
the present paper differs to some extent from the one proposed at the time. 


2. Experimental Part 


2.1. Samples 


Investigation was carried out on the natural GB in n-type Ge. Samples in the 
form of parallelepipedes of 10 x 2 x 1 mm? were cut from polycrystalline Ge ingots 
carefully purified by repeated zone melting. Such samples presented single planes 
of lineage. The planes were perpendicular with respect to the longer axis of the sample. 
Some measurements were carried out on polycrystalline plates presenting many bound- 
aries. 

Specific resistivity of germanium was 5—40 Q cm. Lifetime of the minority 
carriers ranged from 2 to 100 u sec, according to the sample. The grain boundaries were 
detected by etching with a CP-4 solution. Tin electrodes were welded onto the oppo- 
site faces of the sample. Previous to each measurement, the surface was purified by 
etching in 30% H,Oz, rinsed in distilled water, and dried carefully. 

X-ray investigation revealed the GB of the samples to be generally complex bound- 
aries (5 degrees of freedom). 

In measuring the mean diffusion length, photoelectromotive force and photoconduc- 
tivity, the sample was illuminated with a point light spot of 604 minimum diameter. 
The sample was placed upon a micromanipulator mounted on a microscope support, 
so that it could be shifted with respect to the light spot. The shift was measured with 
a micrometer limit gauge. 


2.2. Effect of GB on Diffusion of Minor Carriers 


Bicrystals were used for investigating the shape of the diffusion curves by the 
Haynes-Morton method (Valdes 1952), for three distinct cases: 


a) with the collector situated at some distance from the GB, and the light spot 
moving in a direction perpendicular to the plane of lineage (Fig. 1a); 

b) with the collector situated near the GB, and the light spot moving along the 
plane of lineage (Fig. 1b); 

c) with both collector and light spot far from the GB (Fig. Ic). 


For various GB, two kinds of diffusion curves were found to occur. Typical 
curves of the first kind are shown in Fig. 2. Curve 2 a was obtained with the light spot 
moving perpendicularly to the GB. In the boundary region, a discontinuity of the col- 
lector signal accounting for the decrease in hole concentration on traversing the line- 
age plane is observed. The curve presents a different slope on either side of the boundary, 
according to the different lifetimes of the holes in either half of the bicrystal. Curve 
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2b is that of the same crystal for the light spot moving along the GB. Here, the slope 
considerably exceeds those obtained on either side of the boundary. The curve-forms 
discussed are not very sensitive to homogeneous background illumination of the sample 


with natural light. 


(o) 


light spot 


collector 


f 


trail of the light 


Fig. 1 


The foregoing results point to enhanced hole recombination on such boundaries. 
From the discontinuity of the diffusion curve on the GB, the rate of recombination 
S per unit area of the lineage plane can be evaluated. It is readily shown that S is 


given by 
2K—1 1 
son (4-7) 
with K = In ”1 
Pa 
wherein: = is the ratio of the hole concentrations on either side of the boundary, 
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L,, L, are the mean diffusion lengths on either side of the GB, respectively, and D, 
denotes the diffusion constant of the holes. 

The index 1 serves for denoting that part of the bicrystal upon which the collector 
is situated. 

For boundaries of this kind, the formula yields values of S of the order of 10?—10* 
cm/sec. l ; 

Boundaries exhibiting such properties will be termed boundaries of the first kind. 

Most of the samples investigated exhibited entirely different diffusion curves. 
When the light spot moved across the GB, an anomaly of the collector signal of the 
type shown in Fig. 3 usually appeared within the region of lineage. The anomaly was 
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Fig. 3 


most distinct at low collector bias. With intense background illumination of the sample, 
the diffusion curve regained its regular shape, and the diffusion length computed 
therefrom approached volume value. The width of the anomaly (defined as the distance 
between the local extrema on the curve) was found to be almost equal to the diameter 
of the light spot. The anomaly proved to be associated with the change in potential 
of the collector due to the photovoltaic effect on the GB (cf. 2.3). For comparison, Fig. 4 
brings the graphs of the photoelectromotive force and collector signal as measured 
for the same boundary. i 

The effect of superposition of the photo-EMF on the signal proper could be elim- 
inated by enhancing the sensitivity of the collector e.g. by forming the contact elec- 
trically or by raising the bias voltage. Characteristically, on elimination of the effect 
of photo-EMF, curves of this type no longer revealed a noticeable discontinuity in the 
peak boundary region indicative of a higher rate of hole recombination. 
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Fig. 5 shows the diffusion curves obtained in shifting the light spot along the GB. 
These reveal a considerably lesser slope as compared with similar curves from the 
bulk of the crystal, corresponding to greater mean diffusion length as measured along 
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the boundary. The region of greater diffusion lengths extended to either side of the 


boundary by a width of the order of L. 
Intense background illumination of the sample reduced the diffusion length 


observed to its normal bulk value. A similar decrease (leading to values that were even 
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lower than normal) was also found to occur if a voltage of about 1 V was applied to the 


sample for biasing the boundary. l l 
Boundaries such as those described above will be termed boundaries of the 


second kind. 


2.3. Photovoltaic Effect on GB 


The photovoltaic effect across the leads of the bicrystal on illumination with 
a light spot was measured by a compensation method. The effect was found to be 
present in the case of boundaries of the second kind only. A typical graph of the photo- 
EMF as versus the distance of the light spot from the GB is shown in Fig. 6 (in semi- 
logarithmic scale). As the light spot approaches the GB, the photo-EMF increases, 
the slope corresponding to the mean diffusion length for the holes. When the light 
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spot crosses to the other side of the boundary, the photo-EMF changes its sign. The 
illuminated part of the bicrystal always achieves a potential that is negative with res- 
pect to that of the part in the dark. The width of the region wherein inversion of the 
sign of photo-EMF occurs is limited by the diameter of the light spot. 

Average photo-EMF observed was of the order of several millivolts!. At very 
intense illumination, the figure for many samples rose to 20—30 mV (saturation). 


2.4. Light Modulation of Current across GB (Photoconductivity) 


Variations of the bicrystal resistivity arising from illumination of the lineage 
boundary were measured with a bridge arrangement. Considerable variations of the 
conductivity were observed to occur in the case of boundaries of the second kind only. 


For boundaries of the first kind, the effect did not noticeably exceed volume photocon- 
ductivity. i 


1 It may be worth noting that the photovoltaic effect generated by the same light spot in 
n-p junctions in germanium (junctions obtained by the alloying method, t = 15 sec) was 10 — 50 mV, 
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Fig. 7 shows typical variations of the sample’s conductivity versus the position 
of the light spot on the crystal (applied bias — about 1 V). The slopes of the photo- 
conductivity curves are found to correspond rather satisfactorily to the diffusion lengths 
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for the holes. Hence it can be inferred that the direct factor modulating the boundary 
conductivity resides in the migration of a hole into the region of the plane of lineage. 
Moreover, the photocurrent at illumination of the GB was measured versus the 
voltage applied to the bicrystal. The curves are of the type shown in Fig. 8. These 
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are characterized by an obrupt peak discontinuity of the photocurrent intesity occurr- 
ing at a given threshold voltage. At voltages exceeding the latter, the photocurrent 
exhibits saturation (in some samples, the current was found to decrease very slowly 
with the voltage, within this region). The value of the photocurrent throughout the 
saturation region was generally some tens of times greater than that of the short current 
in the photovoltaic effect. For comparison, values of the photocurrent and of the short 


current are given for three average samples. 


Designation of GB Short current, Saturation photocurrent, 
in uA in uA 
GB 3 0,1 25 
GB 4 0,6 18 
GB 5 0,2 14 


Fig. 9 shows a family of static current-voltage characteristics, the parameter being 
the boundary illumination. 
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2.5. Additional Experiments on Injected Carriers 


Bicrystals exhibiting strong photoconductivity were used for investigating the 
effect of modulating the GB resistivity with current carriers injected through a point 
contact. The measuring arrangement resembled the one used in photoconductivity 
measurements; however, in the present case the light spot was replaced by an emitter 
at about 0,1 mm from the boundary. Sine voltage of a frequency of 1 kHz was applied 
to the emitter. The variable component of the current through the sample was meas- 
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ured versus the voltage applied to the bicrystal, at steady emitter current. Characteris- 
tics resembling those of photoconductivity were obtained (Fig. 10). Hence it seems 
reasonable that the mechanism of the photoconductivity effect and that of the transis- 
tor effect on GB are identical. 
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A curious photoelectric effect arose when a metal point (collector) of bias nega- . 
tive with respect to the crystal was placed on the GB. The light spot was near the 
collector and the lineage plane. The photocurrent flowing through the bicrystal was 
measured versus the collector voltage, for a steady value of the voltage across the 
sample (Fig. 11). 

At first, as long as the values of the voltage V, are small, the photocurrent retains 
its normal, constant value. At a critical value of the voltage on the collector, a peak 
appears in the photocurrent, in the shape of a sharp maximum. As the voltage on 
the bicrystal is augmented, the maximum shifts towards higher voltages and its 


amplitude is seen to augment. 


2.6. Current-Voltage Characteristics, and. Transients 


Current-voltage characteristics of boudaries of the second kind were investigated 
by the pulse method. The samples were placed in forevacuum. Rectangular voltage 
pulses of a duration of 10 u sec were applied to the sample at a repetition rate of several 
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In Fig. 14, a typical GB current-voltage characteristic is plotted in logarithmic 
scale. Small circles denote the values J, of the steady state current, whereas little 
crossse indicate the amplitude V, of the pulse front. Three characteristic ranges 
can be distinguished on the cuvre: 


A) for voltages below 1 V, the current through the . 
sample grows almost linearly. At first, no peak is 
observed. 


B) from about 1 V up to 10—20 V, there is relative 
saturation of the current, /, increasing very 
insignificantly with the voltage. Here, the initial 
voltage peak appears. Its amplitude increases as 
the square root of the voltage applied to the 


bicrystal.” 


C) for higher voltages, the current I, gradually 
increases, until breakdown occurs. The break- 
down voltages usually exceeded 100 V. 


In samples that had been kept for a time in 
air or in an atmosphere of water vapour, slow tran- 
sient effects of a relaxation time of the order of 0,1 
sec were observed (with an arrangement resembling 
the one of Fig. 12). These were probably due to 
the presence of water vapour channels in the region of lineage. Subse- 
quent to pumping away in forevacuum, the samples no longer exhibi any channel 
effects. 

Let us now consider in some detail the short-time transient related to charging 
of the GB barrier. The system equivalent to the measuring circuit is shown in Fig. 15. 
R denotes the measuring load, C — the input capacitance of the synchroscope (about 
40 pF), whilst r represents the sum of the output resistance of the generator and that 
of the sample (without GB). Herein, the GB is replaced by the resistance R, and capa- 
citance C, in parallel. In the steady state, for sufficiently high voltages, the charge 
Q, bound on the GB (i.e. by the capacitance C,) is proportional to V. As the voltage 
pulse is applied to the sample, the process of investment of electron states on the GB 
sets in. In the present equivalent circuit, this corresponds to charging of the capacitance 
C,. Let us assume that in the process of charging R, C, < RC. This is justified by the 
fact that, as the amplitude of the voltage applied exceeds a minimum value of the order 
of 1 V (Taylor et al. 1952), the potential barrier on the GB at first is destroyed, so 
that R, = 0. The rate of charging of both C, and C is then given by one and the same 


OE aie ee See 8 ey ge, iL Se 
. 2 The dependence Vp ~ V3 is exact up to 10 V. Deviation therefrom at higher voltages was due 


to the use of a high-resistance voltage divider on the generator output, affecting the parameters of the 
measuring circuit. 4 
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It follows that if the constant t, of charging of the barrier is much smaller than 
t, the amplitude of the initial voltage peak yields a measure of the charge bound on 


PAS ng eee 


Electronic Processes at Intercrystalline Barriers in Germanium 619 


the GB. Since the relation V, ~ V* has been found to hold even at t = 10-7 sec, 
the conclusion is that it takes no more than 10-8 sec to charge the barrier; most 
probably, the process requires a much shorter time. 


2.7. Résumé of Experimental Results 


Experimental results point to a twofold effect of the GB on the non-equilibrium 
current carriers. This fact has already been accounted for by the distinction between 
boundaries of the first añd second kind. 


Recombination is the characteristic of the effect of GB of the first kind. This 
is proved both by the discontinuity in the decrease of hole concentration on crossing 
the plane of lineage and by the shorter diffusion length as measured along GB of this 
type. The behaviour of these boundaries can be said to be normal, as any lattice defects 
present within the region of lineage, such as e.g. dislocations, are known to constitute 
highly effective centres for the recombination of minority carriers. 


As opposed to these, GB of the second kind present no recombination properties. 
Indeed, such boundaries exhibit an increase in the diffusion length as measured along 
the plane of lineage. This increase is but seeming, since it vanishes as a result of back- 
ground illumination. 

A second special characteristic of GB of the second kind consists in the intense 
photovoltaic effect thereon. i 

Still another property of lineages of this kind consists in the modulation of the 
GB conductivity by light. Experimental results point to holes entering the barrier 
region as the cause of this effect; on the other hand, the possibility of a process consist- 
ing in the liberation by photons of electrons bound on the GB is negl'gible. Modula- 
tion of conductivity cannot be explained on the ground of minority carriers (holes) 
passing through the potential barrier, since, were this the case, the photocurrent 
at saturation would be comparable to the short current in the photovoltaic effect, 
whereas the former is found to be stronger by one order of magnitude at the least. 
This points to a process of current gain occurring on the GB. The characteristic 
shape of the dependence of the photocurrent on the voltage applied to the bicrystal 
suggests that gain on the GB occurs only if voltages exceeding a critical value are 
applied. 

An attempt will be made at explaining these properties of GB of the second 
kind in subsequent sections of the present paper on the basis of a simple model of 
the GB. 

The effect of a negatively polarized point on barrier photoconductivity modula- 
tion, as discussed in Section 2.5, is as yet obscure, and its explanation requires fur- 
ther investigation. 

The question arises of the difference in structure between GB of the first and 
second kind. X-ray investigation of the GB is at present being conducted by J. Auley- 
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tner and M. Lefeld, and the results are to be published separately. It can already be 
said, however, that GB of the first kind are related to simple low-angle lineages 
(0 < 1°), whereas GB of the second kind are representative of lineages throughout 
a wide range angles of misfit. In addition, this latter group of boundaries comprises 


non-coherent twins. 


Discussion 


3.1. General Remarks; Notation 


The GB plane will be assumed to contain a number of acceptor type quantum 
levels partially invested by electrons when at thermodynamic equilibrium. The local- 
ized negative charge resulting from the electron produces a potential barrier in the 
region. of lineage (Fig. 16 a). The discussion will be restricted to an n-type crystal, 
exhibiting identical physical properties on either side of the GB. In this case, the 
potential barrier is symmetric with respect to the plain of lineage. 

Let us denote the height of the barrier at equilibrium by @p. If the boundary 
is biased, the barrier becomes asymmetrical (Fig. 16 b). Let ®, and Ø, donote the 


Fig. 16 


potential discontinuities on the positively and negatively polarized sides of the GB 
respectively. The potential difference applied V = ®, — @,. | 


The following notations will be useful: 
e — elementary charge, 


k — Boltzmann constant, 
T — Kelvin temperature, 
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n, pı — electron and hole concentration on side 1 of barrier, 

Ng, Po — electron and hole concentration on side 2 of barrier, 

No» Po — electron and hole concentration at either side of barrier, at thermodynamic 
equilibrium, 

n,» P — electron and hole concentration in central part of barrier region, 

n> Po — electron and hole concentration in central part of barrier region, at thermo- 
dynamic equilibrium, 

n; — intrinsic concentration of current carriers, 

Jp1»Jpo — hole current densities on sides 1 and 2 of barrier, respectively, 

JnoJp — densities of electron and hole currents through barrier, 

D — mean thermal velocity of electrons, 

Ha» Hp — electron and hole mobility within bulk of crystal, 

t, — lifetime of holes 

D — diffusion constant of holes, 

L, — mean diffusion length of holes. 


3.2. Currents through Barrier 


For germanium of not too high resistivity, and if the voltage applied is not conside- 
rable, the electron current flowing through the barrier can be expressed similarly as 
in the diode theory of the barrier layer (see, e.g. Heinish 1949). The diode theory 
holds provided the mean free path J of the carriers is not smaller than the barrier 


width W, the exact condition being given by 


ae 

e2@D 
Generally, this inequality is satisfactorily fulfilled in the case of germanium. ` 
The flux of electrons per unit area of the lineage plane directed from region 1 


to region 2 is 


i 


no exp (— B®,) 


| cr 


the electron flux in the opposite direction is 


My exp (— B®.) 


D| Sr 


hence, the resulting electron flux through the barrier is given by 


ji = — emg exp (— BO) [exp (— B V) m 
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For low voltages, V < B7} 


in = e7 "o exp (— BDz) BV (la) 
Eq. (1) does not account for the tunnel effect, electron reflexion from the barrier, 
and recombination in the GB region. Whereas the first and second can be shown to 
play no essential part here, the third of these effects will be included in the considera- 
tions of the present investigation. 

If the predominant part of the fall in the voltage applied to the sample occurs 
on the barrier, the hole current through the GB is determined primarily by generation 
recombination processes outside the barrier and can be computed as in the case of 
an n-p junction (Shockley 1950). The hole diffusion currents on either side of the 
barrier are 


; D, | Pb | 
=e— 1 — — exp (— B® 
Jpı D Po Po P (— £ê 1) 
; D, | Pb | 
=e ——exp (— fð.) — 1 
Jp2 TA Po Po p ( B 2) (2) 


respectively, p, being a function of the voltage applied. 
The condition of continuity of the hole current in the plane of lineage can be writ- 


1 
ten in the form j,, = jpg + Jp» Wherein re j, = Ris the number of holes recombining 


(or generated) per unit time in the GB plane, referred to unit area. Recombination 
(generation) within the space charge region of the barrier will be neglected. 

The rate of recombination of pairs of current carriers in the steady state can be 
quite generally expressed by 


R= r(np — nf) (3) 


the parameter r is dependent on the mechanism of recombination, and is generally 
a function of many variables (such as the concentration of trapping centers, hole and 
electron concentration, hole and electron capture cross sections etc.). For recombina- 
tion on GB, the respective carrier concentrations at the centre of the barrier region 
should be substituted in the foregoing expression: 


R = r( py — 73) 
With n, = n, exp (— B®,), we have 
ee Dia 
jr=eR= ern; E exp (— B®,) — i| (4) 
o ~ ee 


With eqs. (2) and (4), and from the continuity condition we have the following con- 
dition for p, £ 


Po exp (BDz) <p, < 2po exp (BD,) (5) 


—— 
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Hence, in a number of cases, it may be assumed with a sufficient degree of accuracy 
that 


Ps = Po &xP (BP2) 
Moreover, it follows from eqs. (3) and (5) that generation cannot occur in the plane 
of the GB, since here always R > 0. 
The parameter r in eq. (3) can be expressed by a quantity whose physical meaning 
is directly evident, namely, the rate of recombination of the electrons, referred to 
unit area of the Jineage plane, S,. 


n 


To this aim, let us rewrite eq. (3) as follows: 
R = r( p, Any + Mp APs) 
with A n, = n, — np» AP, = Pe — Pro (6) 


Since An, and Ap, are of the same order of magnitude, whereas for sufficiently high 
barriers p, > np» eq. (6) approximates to R Srp, An,. The quantity rp, herein is 
identifiable as the electron recombination rate S, in the lineage plane. Thus, eq. (4) 
now yields l 


j= e Sn = 2 exp (— B®.) — 1| (7) 


Although eq. (7) is a rather gross evaluation, it is quite adequate for the present 
qualitative considerations. 


3,3. Photovoltaic Effect 


We now proceed to discuss the condition for phato-EMF generation in the case 
_of homogeneous illumination of one half of the bicrystal (Fig. 17). 


space charge | 
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Assume light to generate pairs of current carriers in region 2 at the rate g. The 
steady state continuity equations for the holes in the regions on either side of the 


barrier are 


P1— Po ay 


lain 
z i 


wis ape ales ae een 8 
g Tp 5 1V Jp2 (8) 


with j,, = —eD, grad py, Jp- — eD, grad pp 
and the boundary conditions 
P\(— 0°) = Po Pa) = Po + BT 


pı (- K) =p, 4 exp (— BV) 0) 


Neglecting recombination within the barrier region, we have, moreover, 


in (- 2) =i = =ip (10) 


Finally, by eqs. (8), (9) and (10), the hole current through the barrier on illumination 


1S 


D, (Po + gtp) [exp (— BV) — 1 
koa aed an 


If the illumination is not too intense, the variation in electron concentration due 


to light may be neglected, yielding n, = ng = nọ. The electron current j, through 
the barrier is then given by eq. (1). 

Clearly, if the current leads of the bicrystal are open, j, + Jp = 0. This condition, 
together with eqs. (1) and (11), yields the EMF for weak illumination in the form 


kT gt ( Bubana y 
Pa e 2 > 12 
€7 Pe 2 D, pi 0) 
Eq. (12) differs from the similar expression for the photo-EMF generated at an n-p 
-1 
junction by the factor K, = £ + = = z] , which may be rewritten as K,= [1 + 
p Pb i 


+ I,/1,|, with [,/I, denoting the ratio of the electron and hole currents flowing through 
the barrier, for low bias of the boundary. 

The factor K, has the meaning of rate of separation of current carrier pairs on 
the barrier. K, increases with the percent participation of holes in the barrier’s conduc- 
tivity. 

The dependence of K, on the height Ø of the barrier is of the type illustrated 
in Fig. 18. A perceptible photo-EMF will be generated if the rate of separation ap- 
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proaches unity, i.e. if the height of the barrier exceeds a critical minimum value. For 
average samples (L, = 0.1 = 1 mm) this is the case at p, > (10% — 104)ny i.e. when 
the hole concentration within the inversion region of the barrier considerably exceeds 
the electron concentration outside (n-p* -n structure). 

Recombination within the GB region does not affect the form of the K, curve, 
but only tends to diminish the effective generated carrier concentration near the barrier. 
If there is lattice compression on the GB, the condition for EMF generation will be 
fulfilled at lower values of p,. 


g 
In 
= =a) = = 7 q Fe 
4 5 =2 7 0 205 log ip 
Fig. 18 


3.4. Transport of Non-Equilibrium Current Carriers over the GB 


Let us consider pairs of additional current carriers, generated e.g. by light, in 
diffusion throughout the neighbourhood of the boundary. Each hole entering the 
barrier region brings with it a positive charge. In the steady state, this charge has to 
be compensated, either by the capture of an electron in the boundary plane (recom- 
bination), or by the expulsion of another hole from this region (re-emission). 

Let Ø denote the lowering of the barrier’s height in the steady state, due ‘to 
the excess carriers. This lowering results in an increase in the current of holes flowing 
out of the barrier, 6j,., and in that of recombination of the electrons captured on 
the GB, 6), 

The probability for the re-emission of a hole from the barrier, K,, may be defined 


as 


K. = = - 
jpn + Or 


with 


In the case of a small perturbation of the barrier’s potential, 6@ < £1, we have, 


by eqs. (2) and (7), 


= L, TM iG 
E E (13) 
EW p Pb 


If K, = l, a hole entering the region of the barrier interacts electrostatically with 
the Ee holes in the region, causing one of them to be expelled immediately. The 
process resembles the “feed-in, feed-out” effect observed on a pt-n junction (Moore, 
Webster 1955). ; 

The shape of the dependence of the factor K, on the barrier height resembles 
that of the separation rate; albeit, the region wherein K, approaches unity can begin 
at lower values of the barriers. E.g., for S, = 105 cm/sec? L, = 0,l cm, we have 
K, 21, whilst p, > 100 no. : 

Thus, the “feed-in, feed-out” effect must occur on barriers presenting a conside- 
rable photovoltaic effect (GB of the second kind). 

As the probability K, for the recombination of a hole entering the barrier region is 

Ojr 
K, = gh =1-K, a4 
recombination will not occur in the case of such GB wherein the “feed-in, feed-out” 
effect predominates, and vice versa. 

The “feed-in, feed-out” effect, moreover, can yield an explanation of the increase 

in diffusion length along GB of the second kind. 


[px tnwersion layer 
Te othe OB 


Fig. 19 


Let us consider an inversion layer of the boundary, of effective width Wç and 
hole concentration p,. It will be assumed that a hole enters the layer at x = 0 (for 
simplicity, discussion is for the two-dimensional case, Fig. 19). Its charge lowers 
the barrier’s potential at some distance from the point x = 0. Such a perturbation 
of the potential produces, an ohmic current Jpx of the holes along the 0X axis, a diffu- 
sion current Jpz in the direction of OZ, related to the outward flow of holes (re-emis- 
sion); moreover and, a recombination current j,- The steady state continuity equation 


3 The value of S,, was assessed from the effect of photoconductivity, see Section 3.5. 
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for the holes in the inversion region of the barrier is 


Joz sanS ih ei 
Wos E div Jpx A We P 0 (15) 


wherein 
. do D 
Jpx = — &fp Pb Ta Jpz e T PoP? 


Ls denotes the mobility of the holes in the boundary plane, and p = ® — ®, is the 
deviation from equilibrium of the barrier potential (it is assumed that y < f~}). 
The solution of eq. (15) is of the form 


x 
P = Po EXP E z) 


r= | B war] a6 
Hp Po 


Since the concentration of the re-emitted holes in the region adjacent on the barrier, 
Ap = p — P,» is proportional to œ, we have 


Ap ~ exp (= a 
Thus, the parameter / has the meaning of the average distance at which re-emission of 
holes from the barrier occurs. The quantity 4 may be identified with the apparent 
diffusion length of the holes as measured along the GB. 
Assuming n, = 10%cm-8, Lp =10-tem, W= (100 — 1000) A, £s =100 


ó 


with 


and Mp = up, we have 
es 10 + 100 
L, 

Thus, for barriers of the “feed-in, feed-out” type, the apparent diffusion lengths 
measured along the GB can be many times greater than the bulk L, values. 

Experiments by the present author pointed to a decrease in 4 resulting from 
strong background illumination or from a bias voltage applied on the GB. This is — 
reasonable, since both background illumination and polarization reduce the ratio 


Pe , rendering “feed-in, feed-out” less effective. 
No 


3.5. Photoconductivity 


Holes generated by light in the neighbourhood of the GB can pass freely through 
the biased barrier, yielding the primary photocurrent 6j,,. On the other hand, 
they are able to modulate the height of the barrier and, consequently, the intensity 
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of the electron current flowing through the GB. Let us denote the perturbation in the 
barrier potential, resulting from the incoming carriers, by 6@. This perturbation 
causes the electron current to vary by the amount Ôj,- 
The quantity 
Ko le 
Ope ai Oj, 

characterises the current gain on the barrier. By eqs. (1) and (14), rate of gain the 
gain coefficient can be written as follows: 


K = xK, 
wherein x = ——, for applied voltages V > p=. 
AS, 


Hence, the total photocurrent flowing through the barrier, dj, is given by 
Oj = ja (l + xK,) 


The gain in current on the barrier depends on K,. If K, œ 0 (i.e. if the “feed-in, 
feed-out” process is effective), we have K,, = 0, and the total photocurrent through 
the barrier is equal to the primary current of light generated holes flowing into the GB 
region. As K, increases (e. g. due to voltage applied on the GB), the secondary elec- 
tronic current plays an increasingly important part in the photoconductivity of the 
barrier. At sufficiently high voltage, when each hole entering the barrier region 
undergoes recombination (K, œ 1), the rate of gain reaches its maximum value x. 
The photocurrent then attains saturation: 


Sjoat = Sion ( i a) (17) 

The mechanism proposed here explains the experimental voltage dependence 
of the photocurrent for GB of the second kind. The characteristic discontinuity in 
photocurrent intensity occurring at a given value of the voltage applied would thus 
be associated with a transition from a barrier of the “feed-in, feed-out” type to one 
of the recombination type. 

If the total photocurrent at saturation 6j,,, and primary photocurrent Oj, are 
known, eq. (17) yields the surface recombination rate for electrons in the plane of 
lineage, S,. The primary photocurrent can be reasonably assumed to be of the order 
of the short current in the photovoltaic effect. The values computed for S, on this 
assumption are of the order of 105 — 108 cm/sec. 


4. Conclusions 


GB can affect non-equilibrium current carriers in a twofold manner. The height 
of the barrier or, more exactly, the K, value, is the critical factor. 


In the case of a low barrier, i.e. for K, 0, recombination on the GB predomina- 
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tes. Such GB exhibit no appreciable photo-EMF, but should present photoconductivity 
(current gain on the barrier). 

In the case of sufficiently a high barrier, for K, = 1, the properties of the GB 
are determined by a “feed-in, feed-out” effect, consisting in electrostatic interaction 
between holes in the inversion layer of the barrier. The effect provides for immediate 
expulsion of a hole from the barrier region for every incoming surplus hole. In the 
case under consideration, transfer of non-equilibrium current carriers through the 
GB is no diffusion process. Such a boundary will fail to exhibit recombination 
properties. 


For average samples, the condition K, = 1 is fulfilled if Po ~ 100 i. e. if the 
No 
boundarie» have the n-pt-n structure. On such boundaries, considerable photo-EMF 
can be generated. 


Voltage applied on GB reduces the ratio = , and hence leads to a decrease in 
0 
K,. Thus, a threshold voltage exists at which transition from a “feed-in, feed-out” 
type barrier to a recombination type barrier occurs. 

The GB of the second kind observed in the course of the present investigation 
can be said to be identical with the high barrier GB of the “feed-in, feed-out” type. 
From a structural point of view, the present grain lineages of a large angle of misfit. 

According to the model proposed, boundaries of the first kind may be said to 
be related to low barriers of the recombination type. However, since no photoconductiv- 
ity was observed on such boundaries, neither did they present a distinct fall in voltage 
as current flowed through the bicrystal, it is thought that they should rather be related 
to a dislocation structure. In the case of dislocation low-angle boundaries, the various 
electron levels may be so far remote from.one another that no continuous potential 
barrier will arise within the plane of lineage. 

The fact that the photo-EMF generated on the GB is strongly dependent on the 
position of the light spot, throughout the region of sign inversion suggests, the con- 
struction of a photoelectric element presenting a high degree of sensibility to small 
displacements of the light spot. | 

Similarly, the effect of photoconductivity on GB can be employed for constructing 
a photoelement of high quantum yield. 

Photoelectric effects on GB probably play an essential part in the well-known 
photosensitive layers (infrared detectors). In particular, the present author (1959) 
attempted to explain the high, isotropic photo-EMF in directionally evaporated PbS 
layers in terms of photovoltaic potential differences generated on the intercrystalline 
barriers. The present analysis of intercrystalline barriers will, in general, also hold for 
barriers occurring on the surface of a semiconductor. E.g., for sufficiently high bar- 
riers, a surface “feed-in, feed-out” might be predicted. Certain experiments carried 
out by M. Jastrzębska inthis laboratory point to such a process as being the reason 
for the anomalously great diffusion lengths found in silicon. 
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QUENCHING OF PHOTOLUMINESCENCE OF SOLUTIONS 


By C. BOJARSKI 


I Department of Physics, Gdańsk Technical University 


(Received February 12, 1960; revised manuscript received March 18, 1960) 


Accounting for the effect of quencher molecules beyond the active sphere on the 
excited luminescent molecule, and for fluctuations in their concentration, the accuracy 
of the model implying the active sphere of Jabłoński’s theory of the quenching of photo- 
luminescence of solutions was improved. The expression obtained for the relative yield 
is compared with experimental results of other authors. 


At present, several theories on the quenching of photoluminescence of solutions 
exist (Vavilov 1928, 1936, Vavilov and Frank 1931, Sveshnikoff 1935, Förster 1949). 
The theory of the quenching of photoluminescence recently proposed by Jablonski 
(1954), in distinction of earlier theories by other authors, accounts i.a. for the fluctua- 
tions of the concentration of the quenchers (quenching molecules). Therein is intro- 
duced the so-called “many-shell model” of the luminescent centre, consisting of an 
excited luminescent molecule enveloped by monomolecular layers of solvent molecules, 
among which some quenchers may be present. In a given solution, the various lumi- 
nescent centers differ from one another by the quenchers being differently situated 
within the shells. The probability of finding a luminescent centre with a definite 
distribution of the quenchers in the separate shells is given by the generalized Smolu- 
chowski distribution (1904, 1907): ~ 


P (ky ko.. ka) = [J om a 3 (1) 

l=1 Is t 
wherein z denotes the number of shells assumed, k, the quencher population of the 
shell 1, and v, the-average number of quenchers in the same shell, with the relation 
vı = v. n, Wnerein n is the number of quenching molecules per unit volume and v, 
the volume of the shell /. Basing i.a. on the foregoing assumptions, Jabłoński gave 
the following expression for the mean relative yield in the case of solid solution: 


3 Si 
i = Ni P (kis kake ka) aE SER (2) 
‘eee 1 AF y + W, 24 kl 


(631) 
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where 7) denotes the yield of a solution containing no quenchers, y — the proba- 
bility per unit time of an electronic transition connected with the emission of photolu- 
minescence, W, — the probability of inner quenching, and W, is the probability 
of quenching by the k, quenchers situated in shell ?. 

The theory, originally restricted to a single shell only, resulted in the following 
very simple expression for the relative yield: 


Toa l— e” (3) 

No ig 
where v = vn denotes the mean number of quenchers in an active sphere of volume v 
at the concentration of n cm~? of quenchers. Subsequently the theory was improved 
by accounting for the dependence of the quenching probability on the distance R 
between the excited and quenching molecules!, and for quenching by quenchers 
situated outside the so-called active sphere. The relative yield of a solution as a function 
of the quencher concentration is now computed from the following expressions: 


i= >) Pe th (4) 
k=v 
wherein 
yk 
P, = Cimi (5) 


is the Smoluchowski distribution, and n, the yield of the group of centres with k 
quenchers within a sphere of radius R’; moreover, 


Nk Nk Nk—i 
i (6) 


wherein 


ay 

= arct aS 

Nk Euj : Z) 
Ok—1 


1 According to Förster (1949), the probability for quenching is inversely proportional to the sixth 
power of the distance R separating luminescent molecule and quencher in cases in which quenching is 
due to excitation energy transfer from luminescent molecule to quencher, thus 


Leia. vo 
r-~(¥}, ©) 


where t is the mean lifetime of the excited state of the luminescent molecule with quencher molecule 
concentration tending to zero. 


(7) 
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together with 


no _ 1l 
ne ° 


The last expression denotes the relative yield of a centre with no quenchers 
(k = 0); R, is the distance between excited molecule and quencher at which the probabi- 
lity for quenching equals that for emission, and R’ = 1,327 Ry, 04-1 = Ry (Np-1) ° 
This new version of the theory proved to be in good agreement with the experimental 
results, except for the experimental points corresponding to the highest concentrations 
of the quenchers (Jabłoński 1957, 1958; Głowacki and Pohoski 1959). According 
to Jabłoński (1958), a better approximation can be attained by applying the “many-shell 
model” of luminescent centre. In this paper, however, the problem under considera- 
tion will not be solved in the manner mentioned above; the method will consist in 
an improvement of the accuracy of the model with an active sphere. It seems that 
this aim can be attained by accounting not only for the effect on a luminescent mole- 
cule of quenchers situated outside the active sphere, but also for the fluctuations of 
their concentration by application of the generalized Smoluchowski distribution 
given by eq. (1). For this purpose, certain assumptions admitting of the computation 
of the relative yield as a function of the quenching molecules concentration by means 
of eq. (2) will be introduced. First of all, the dependence of the quenching probability 
on the distance of eq. (9) is assumed. The space surrounding an excited luminescent 
molecule is divided in concentric spherical layers of equal volume, which may contain 
quenching molecules. The volume of the first layer is chosen as that of the so-called 
active sphere assumed in the simplified model of a luminescent centre. With the 
assumption of such layers the “thickness” of the layer with the number z tends towards 
zero for z -> oo. The division of the neighbourhood of an excited luminescent molecule 
in layers as carried out here provides a way of considering fluctuations in quencher con- 
centration also beyond the active sphere, and the assumption of equal layers simplifies 
the calculations essentially. Obviously, such handling of the problem does not coincide 
with the many-shell model introduced by Jabłoński (1954). In.order to determine 
the probabilities w,,, it will be assumed that for a given layer the quenching probability 
is the same for all quenchers present therein, but that it is different for different layers. 


Thus, 
l= 


De ez, 
2 kd 10 
k= 0, 1 eae ( ) 


12, +4. 00 


LB desea LG 


wherein W, denotes the quenching probability by one molecule situated in the layer 
l, and k; is the number of quenchers in the same layer. By the assumption of v; = 
= v (l = 1, 2, ... z), where v, is the volume of the layer l, 


j Ry= Ry yi (11) 
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where R, denotes the radius of the outer surface limiting the layer l. By eqs. (9), (10) 
and (11), 


W. 
Wa = ki se l= ERA A (12) 


where W} is the quenching probability by one quenching molecule present within the 


first “layer”. Assuming VW, = y + Wy = = „with Tọ denoting the duration of photo- 


luminescence for >) k= 0, and by eqs. (12) and (2), the following expression for the 
i=1 


relative yield in the case of a luminescent centre with z layers is obtained: 
7 oe Se hs 1 
N ay aes ate 
gas 2 2 II kı! > ky (13) 


where y = v n is the average number of quenchers in volume v, and n their average 
number per cm. In principle, expression (13) can serve for computing the yield for 
a luminescent centre with an arbitrary number of layers z. Clearly, the greater is 
z the more tedious are the calculations. In practice, however, neglecting layers with 


numbers z > 3 does not influence the dependence of A as a function of v. Considering 
0 
the dependence of the quenching probability on the distance as determined by formula 


(9), it would seem obvious that only the quenchers in the nearest neighbourhood of 
the active sphere exert any considerable influence on the excited luminescent molecule. 
Fig. 1 represents a graph of the relative yield as a function of the logarithm v, as 
determined by eq. (13) for z = 3 (graph III). For comparison, Fig. 1 also contains 
graph I as given by eq. (3), being the case of z = 1, and graph II as given by eqs. (4), 
(5), (6), (7) and (8). Graph III obtained in the present case is shifted with respect to 
graphs I and II and presents a somewhat different slope, especially at high concentra- 
tions (high v) viz., it intersects graph II, decreasing more slowly than the latter. It 
should be stressed that eq. (13) contains only one constant, which may be determined 
by comparing the theory with experimental results. This constant is the volume v of 
the first “layer” of the luminescent centre. In Fig. 2, the expression of eq. (13) is 
compared with the experimental results obtained by Férster (1949), relating to the 
quenching of the fluorescence of trypaflavine in methanol by rhodamin B. The agree- 
ment is very good throughout the entire range of concentrations. The value of the 
radius of the active sphere is 68 A. Fig. 3 brings the dependence of the relative yield 
given by eq. (13) as a function of the logarithm v, and the experimental results obtained 
by Głowacki and Pohoski (1959) relating to the quenching of the fluorescence of 
solutions of trypaflavine and fluorescein in ethyl alcohol by rhodamin B. In this 
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Fig. 1 WE versus logarithm of v. 
No 


—----- graph I, from eq. (3). 
siie graph II, from egs. (4), (5), (6), (7) and (8). 
graph Ml from eq. (13). ; 
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HE =74 0 1 logv 


Fig. 3. A versus logarithm of v. 
n 


theoretical curve computed from (13). 
O, A — experimental points according to Głowacki and Pohoski, relating 
to the quenching of trypaflavine and fluorescein in ethyl alcohol by rhodamin B. 


case too, the agreement as between the experimental results and the theory is very 
good with the exception of the experimental points corresponding to the highest con- 
centrations of the quenchers?. The mean values of the radii of the active spheres are 
54 Å and 76 Å for trypaflavine and fluorescein, respectively. 
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(Received February 15, 1960) 


Measurements of the exponent of the density spectrum for the electron-photon 
component of the extensive air showers of cosmic radiation were made by means of a 48 
channel hodoscope connected with 6 counter sets. The results obtained show that within 
the limits of experimental error the exponent does not vary with the thickness of a lead 
absorber (in the range 0—27 mm). 


1. Introduction 


In the case of pure electron-photon cascades and a development parameter of 
s > 0.6, the theory predicts an excess of the number of photons over that of the electr- 
ons. From numerous experimental data on extensive air showers, the respec- 
tive electron-photon component for altitudes from mountain to sea level is known 
to behave as in a pure electron-photon cascade of development parameter s = 1.2—1.3. 
Hence, the prediction od cascade theory relating to the excess of photons over elec- 
trons should be expected to hold in the case of extensive showers also. 

The presence of photons in extensive air showers is detected from the transi- 
tion effect of the electron-photon component in a Pb absorber..The theoretically 
predicted excess of photons over electrons, and the theoretical transition curves for 
electrons P, (t) and photons P;(t) computed by Arley (1948) had lead to the conclusion 


that the transition curve for the electron-photon component R (t) = P, + 2 P: should 


exhibit a large maximum. P, and P, denote the probabilities for at least one electron 
to be present under a Pb absorber of thickness ¢ if a photon or electron, respectively, 
was incident thereon, and fje is the photon/electron ratio. Alas, all experimental 
investigations on the electron-photon transition curve have yielded an exceedingly 
small maximum. For this reason, early investigators concluded that a photon/electron 
ratio in extensive showers is much smaller than unity (Bassi and co-workers 1952, Millar 
1951), or one smaller than unity (Massalski 1954, Milone 1952a, 1952b, 1953a, 
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1953b, 1954, Babecki and co-workers 1957). In the course of investigations 
carried out at this laboratory (Dubinsky and co-workers 1957, Massalski and Oleś 
1958), attention was drawn to the steepness of the P, curve resulting from the soft 
part of the electron spectrum as being the reason for the low transition effect, and to 
the fact that photons contribute to the recording of the electron-photon component 
even at zero thickness of the Pb absorber above the counters. With these two factors, 
agreement was obtained between the theory and the experimental data. 

The foregoing authors, in measuring the photon/electron ratio in extensive 
air showers, had assumed a constant value of the exponent of the density spectrum 
under the absorber. Janossy and co-workers (1957) and Somogyi (1957) pointed 
to the possibility of this assumption being incorrect, i.e. to that of an exponent 
varying with the thickness of the absorber. This, in turn, would lead to a revision 
of the counter method of determining the photon/electron ratio. It was the aim of the 
present paper to investigate this problem. 


2. Apparatus 
In considering measurements of the transition curve, it is convenient to introduce 
the reduced density. Let denote the density of the particles incident on the absorbent; 
then, beneath a layer of the latter of thickness t, the reduced density of the particles 
amounts to 


e= Rt) =o |e. (t) + fro] : 


e 


Thus, the number of coincidences of n trays, each presenting the area Q, situated 
under a Pb absorber of thickness ¢, can be written as follows: 


T fK (1 ae e~ QeR)n do, N (1) 
0 


wherein Kg~‘’* do is the differential density spectrum of the showers recorded. Herein, 
the reduced density @ has been assumed to possess a similar spectrum 


0 7 OLD dp is Bu 


Eq. (1) or similar equations have been utilized in the papers cited above for determin- 
ing R and, subsequently, f/e. The number of coincidences of n uncovered trays can 
be written as follows: 


= f Ko™®+9 (1 — e2)" do, (3) 
. 0 
whilst, by eqs. (1) and (2), that of covered trays is: 


4 


RE f Ko-?t as e~ 20)n do. (4) 
0 
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Eqs. (1) and (4) hold provided the lead absorber does not affect the space angle of 
the measuring trays. This latter condition can be fulfilled if the trays are placed within 
a lead “well” with walls exactly bounding the active surface of the counters, and if the 
lead absorber is situated immediately above the measuring trays. In order to compare 
the coincidences at different area of the trays, as required for determining the exponent 
of the density spectrum under the lead absorber, the smallest surfaces to be accounted 
for, as e.g. single counters in the present case, must be shielded similarly. Trays of 
greater area are obtained by connecting in parallel several identical smaller ones. Thus, 
identical measuring conditions and space angle are provided for all the areas 
compared, from the smallest to the largest. The various counters are reached by local 
cascades from an area of the absorber approximately equal to the active area of the 
counter. These cascades hit only one counter simultaneously. The foregoing conditions 
provide for the device to record the reduced density. 

In the present investigation, 48 hodoscope G.-M. counters of an area g = 450 cm? 
each were used, these being situated in three uncovered trays A, B, C, and in three trays 
D, E, F shielded from below and laterally by a5 cm Pb layer and covered from above 
with a Pb layer varying in thickness from 0 to 27 mm. The counters in the trays A, B, C 


tray D tray A 
7 LLLLLLLLLL LLL), 
VA ZA — 
aaa A ———— 
va CAS 
4 [ae | 
| es | 
SS 
tray 8 tray F tray C 
YLLLLLIILLELD LY, 
—_—— ee ee | a ee ene | d 
e] was RAEE] Cee 
E GLL, < 
(aee | 
(a | 
———= 
_————— a | 
——li 
_— ee | 
——I 
ee 
(n | 


Fig. 1. Arrangement of counter trays. 


were situated 5 cm apart from one another, whereas the various counters in trays D, E, F 
were separated from one another by a Pb layer 5 cm thick. Fig. 1 illustrates the arrange- 
ment of the trays, and Fig. 2 — the vertical section through one of the A, B, C and 
through one of the D, E, F trays. Simultaneous measurements recorded ABC and DEF 
coincidences, which provided the master pulses for the hodoscope. Measurements 
with the A, B, C trays yielded the exponent of the density spectrum of the electron 
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component. Data from the D, E, F trays yielded the dependence of the exponent of 
the reduced density spectrum (i.e. of the density spectrum of the electron-photon 
component) on the thickness of the absorbent. Thanks to good statistics, the results 
obtained from the A, B, C trays served for comparison with those from the D, F, F tray 


measurements. 


Or O7G"'@ "OF OF CFO 


0O 70 20 3040 50cm 
t< 


Fig. 2. Vertical section through one of the D, E, F, and one of the A, B, C trays. 


3. Experimental results 


wy Dependence of the exponent of density spectrum on absorber thickness 


The exponent of the density spectrum can be computed from the data of hodoscope 
measurements, utilizing the formula 


Q (5) 
‘ ($) 


and combining halves of the A, B, C and D, E, F trays in various ways. In eq. (5), 
I(Q) and I 9 denote the number of coincidences of arbitrary multiplicity (in the 


present case, triple coincidences), for a tray of area Q or Q/2. Eq. (5) results from (3) 
or (4) and holds both for uncovered and absorber-covered trays if the conditions 
of Section 2 are fulfilled. Fig. 3 presents the values of the exponent y as computed for 
the trays A, B, C, and also for D, E, F at various values of the Pb-absorber layer. 
From the results obtained, the exponent of the spectrum, throughout the range of 
0 to 27 mm Pb, is seen to be practically independent of the absorber thickness. 


3.2. Dependence of the exponent of density spectrum on the number of counters hit, at 
various values of the absorber thickness 


Hodoscope measurements yield the number of triple coincidences P(m) for m 
counters hit in one tray and at least one counter hit in each of the remaining trays. 
This number can be expressed by the formula 


I(m, y) =K k fo = e-ae)m e—(M—m)ae (l =: em Mao)2 o7Y+D do, (6) 


0 
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Fig. 3. Computed values of density spectrum exponent y for trays A, B, C and trays D, E, F, versus Pb 
absorber thickness. 


wherein q denotes the active area of the G.-M. counter, M = 8 is the number of 


counters per tray, and ey — the number of combinations of hittings of m of the M 


counters. For the D, E, F trays, the quantities K and g should be replaced by K and 0 


respectively. 
Integration yields 


I(m, y) = Kq’I(— y)P(my), (7) 
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wherein T'(— y) is the Gamma-function, and 


m 


P(m, ») = (—1)" 9° (1) o [(M—1)’—2(2M—D’ + (3 M—)). 


1=0 


The values of the function T'(— y) P(m, y) were computed for the device applied i.e. 
for M = 8, m = 1, 2, 3, ... 8, assuming the values of y = 1.3 and 1.5, as the quantity 
y required ranged within these limits. The functions I(—y) P(m, y) intersect at one 
point my = m(M) independently of the value of y (Zawadzki 1957). Hence, the exper- 


imental curve p(m) can be normalized so that, for m = mọ, 


C p(m) = p,(™m) = T (—y) P (mo, y), 
with C denoting the normalisation constant. 

Fig. 4 presents the theoretical lg P(m . y) (—y) curves computed for the two 
values y = 1.3 and 1.5, and, moreover, the experimental lg p,(m) curves determined 
for DEF coincidences of the trays under a 27 mm Pb, absorber layer and for ABC 
coincidence. 

The curves obtained with the D, E, F Pb-covered trays are seen to present a form 
similar to that obtained with the uncovered A, B, C trays. Hence it may be concluded 
that shielding the counters with lead and placing a lead absorber 27 mm thick above 
them fails to affect perceptibly the dependence of the density oe exponent 
on the density ọ of the particles incident. 


3.3. Dependence of the exponent of the electron-photon component aney spectrum 
on the energy of the particles recorded 


The exponent of the density spectrum as obtained by the half-tray method at 

zero absorber thickness is somewhat higher (1.41) for the trays 4, B, C than for 
` D, E, F (1.38), the counters in the latter being separated from one another by lead 
partitions. The reason for this can reside in a photon transition effect in the lead parti- 
tions. At zero absorber thickness, the trays A, B, C record primarily the electron 
component, whereas the recordings of D, E, F contain a considerable photon contribu- 
tion. For the above reasons, the intensity of the coincidences recorded at zero absorber 
thickness by the D, E, F trays exceeds that of the 4, B, C ones. The per cent contribu- 
tion of the photons in the recordings of extensive showers by the D, E, F trays increases 
with the thickness of the absorber. On passing the maximum of the transition curve, 


this contribution is comparable to that of the electrons, since L = l in the expression 


for the reduced density 0 F a5 2 P.) ‘recorded under the absorber, whereas, on 


having passed the maximum, P, ~ ae The forms of P, and P, are determined by the 
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Fig. 4. Theoretical lg P/m, y/I"(— y) curves for y = 1.3 and 1.5, and experimental Ig p,(m) curves 
versus the number of counters hit. 


energy spectrum of the photons and electron recorded. An increase in absorber thick- 
nessis equivalent to passing to higher energies in the energy spectrum of the electron- 
photon component. In recording the reduced density under an absorber of varying thick- 
ness, we in fact record the electron-photon componentat various threshold energies of 
the device. The experimental result of a constant value of the exponent of the density 
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spectrum of the electron-photon component points to an electron-photon component 
density spectrum having, within experimental error, the same exponent as the electron 
component spectrum, and to an exponent independent of variations of the threshold 
energy of the particles recorded. 


3.4. Conclusions 


The results of Sections 3.1 and 3.2 show that, in recording the reduced density, 
the exponent of the density spectrum is found not to depend on whether the measure- 
ment has been carried out with or without absorber, and to be independent of the 
thickness of the latter. 

In order to compare the results obtained by the Budapest group with those of 
the present authors, the arrangement of the trays D, E, F was changed to agree with 
that of Janossy and co-workesr. Fig. 5 illustrates the new arrangement of the counters 
in the trays D, E, F. The density spectrum exponent computed for the device thus 
modified and utilizing the half-tray method at 25 mm Pb absorber thickness amounted 
to y = 1.12, a value resembling that of 1.18 obtained by Somogyi. Moreover, the 
lg p, (m) curve for the geometry of Fig. 5 differs considerably from the curves obtained 


0 10 20 30cm 
ee ee eee 


Fig. 5. Vertical section through one of the D, E, F trays for modified arrangement of counters 


formerly and points to a lesser value of the exponent y (Fig. 4). The divergence in the 
results obtained arises from the geometry of counters and absorber. When the coun- 
ters are separated from one another by a Pb partition and the absorber is situated 
immediately above the counters, a local cascade initiated by a single electron or photon 
hits but one counter. In this case the device records the reduced density. Experiment 
proves that, in this case, the exponent of the density spectrum does not vary with 
the thickness of the absorber i.e.. with the threshold energy of the particles recorded. 
On the other hand, if the device is sensitive to the number of particles in the local 
cascades, measurements reveal a density spectrum exponent that varies with the thick- 
ness of the absorber, since the local cascades contribute to an increase in the density 
of the showers recorded leading to a smaller value of the spectrum exponent. If the 
device is modified by removing the absorber to a greater distance from the counters, 
and‘if the counters are placed one beside the other, a local cascade hits more than 
one counter and the device yields a density in excess of the reduced density, as was 
the case of the Hungarian authors. The role played by the absorber may be said 
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to be that of a factor modifying the effective area of the counters or of the entire tray. 
The variation of the effective area depends on the thickness of the absorber, and will 
increase with the ratio of the areas of absorber and counter (or absorber and tray) 
and with the distance separating absorber and counters. 

The authors wish to thank Professor Dr M. Miesowicz and Docent Dr A. Zawadzki 
for their valuable discussions throughout the present investigation. The authors 
are sincerely indebted to their colleagues of the VI. Laboratory of the Nuclear Research 
Institute, Cracow Department, for their kind help. 
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The resistivity of a disordered binary alloy is composed of three terms: a vibrational 
term Q,, a disorder term gg and a mixed term @,,. We have devoted a special attention 
to the mixed term, bringing also certain improvements to the term Q, and ọ4. Finally, 
Om has been compared to g,, and gq, and deviation from the Mathiessen rule has been 
discussed. 


The electric resistivity of an alloy is due to the interaction of the conduction elec- 
trons with the thermal vibrations of the lattice, on one hand, and with the ions which 
disturb the periodicity of the ordered lattice, on the other hand. Thus the resistivity 
o is split into two terms: a purely vibrational term ọ, and a residual term @,, the index 
d of the latter-being meant to indicate that this term is due to the disordering of the 
lattice. The Mathiessen rule can therefore be written in the form @ = @, + 04- 

In fact, the resistivity of an alloy is greater than @, + 04 (Kohler, 1959). Consequ- 
ently, we may write 


0 = 0, + Oa + Om (1) 


where @,, is a mixed term which is due both to the thermal vibrations and to the disor- 
dered state of the lattice (Corciovei, 1959b). The various terms appearing in formula 
(1) have been studied by different authors. We do not find it necessary to discuss the 
purely vibrational term @, (Peierls 1955). The purely disorder term was the object of 
theoretical research, especially in recent times (Muto 1936, Murakami 1953, Krivoglaz, 
Matysina 1955, Gibson 1956, Corciovei 1959a). As for the mixed term, it has been 
less thoroughly studied than the former two. In the present paper, in which the nota- 
tions used ‘are taken from a previous paper of one of the authors (Corciovei 1959b), 
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We eee 


we have devoted special attention to the mixed term g,,, bringing certain improvements 
to the terms g, and @, too. Finally, ọ„ has been compared to ọ, and gg, and it has also 
been shown in how far deviations from the Mathiessen rule occur. 


1. The perturbation Hamiltonian 


We have chesen as the object of our study a binary alloy with an equal or almost 
equal number of type A and type B atoms, free of any impurities. We assume that 
the crystal lattice is such that it can be split into two identical sub-lattices a and B», 
which are shifted with respect to each other. When the order is perfect and if the num- 
ber of type A atoms is strictly equal to that of type B atoms, all the points of the 
a-sub-lattice are occupied by A-atoms, while all the points of the f-sub-lattice are 
occupied by B-atoms. 

We shall denote with j an arbitrary lattice point of the a-sub-lattice and with k an 
arbitrary lattice point of the B-sub-lattice. Finally we shall assume that each atom yields 
a single electron in the conduction band, the potential of this electron being U, (r) 
in the field of the A-atom, and U,(r) in the field of the B-atom. If the parameters 
Lj and u, which are defined as follows 


1 a 
Lj = {* 1 if j contains an atom of type f ua E > if k contains an atom of type B 


(2) 


are introduced and if u, and u, denote the deviations, which are time dependent 
and correspond to the vibrations of the ions located in points j or, respectively, in 
points k, from their equilibrium positions r, and r,, then the potential in which a con- 
duction electron. is moving is given by? 


U(r)=} o> [0 + 4) U,(r —r;— u) E) Ug (r — r; —u;)] =f 
ae pa [1 — uy) Ua (r — r, — u) + (1 + u) Up (r — r, — 4,)] (3) 


We shall consider as an unperturbed state the perfect order state, in which the 
atoms are in their rest positions. The potential corresponding to this state is obtained 


by substituting in (3), uj = 1, p, = l and u, =u, = 0 for all j and all k. It then 
becomes 


Uslr) = 31 Vale — 4) + 3) User — 1) (4) 


1) Such as for instance, the simple cubic lattice, the cubic centered lattice, and the NaCl type lattice. 
In the following we have usually assumed that the lattice is cubic centered, the a-sub-lattice being the 
one consisting of the corners of the cubes, while the B-sub-lattice contains the centers of the cubes. 
*) Vectors are denoted in boldface type. 
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This potential is periodical with respect to any displacement corresponding to the 
distance between two points of the first sub-lattice or of the second sub-lattice, respec- 
tively, thus: 


Uy [r — (r; —7;)] = U lr — (r, — ry)] = Uo (r) (5) 


For the monoelectronic wave function corresponding to the potential (4), the 


Bloch theorem may be applied. If we consider a single band, we may write the wave 
function y, corresponding to a particular propagation vector h® as follows 


Vn = V5 eihr u (r) (6) 


where u(r) is subject to the same periodicity rules as Ug (r) in (5). It is normed in an 
elementary cell® and its dependence on h has been neglected. 

The resistivity of a perfect periodical lattice is known to be null. The resistivity 
of the alloy considered is different from zero because of the deviation from periodicity 
due to the thermal vibrations and to the disorder, that is owing to the existence of 
a perturbation hamiltonian H?* = U — Up. Calculations carried out elsewhere 
(Corciovei 1959 a) show that HP® is the sum of three terms: a purely vibrational 
term H”, a purely disorder term H’, and a mixed term H”. We may, therefore, write 


He WY PH + Be, (7) 
where 
H = — Siu, W^ (r — r) — Diu, W? (r —7,) 
j k 
Ht = Yu) Uer) I Uer (8) 


G 2i (1 — pu, W (r —r;) + pa (1 — u) U, W (r —7;) 


the notations adopted being 
W4 = grad U, WP? = grad Up U= 3 (Ug — A U E (Wr — w4) (8) 


In the following we.shall calculate the matrix elements of the perturbation hamil- 
tonian HP® between the states (6), and thereafter, to find the resistivity, we shall 
apply the formula (Gibson 1956) 


| Sa Sha ie 
e = Keav (n)? Í [ARR |2 (L — cos 8) dQ, dQy 0) 


ee aie ae a 
3) h is a vector in the reciprocal space which takes N/2 values in a band, N being the total number 
of the atoms of the alloy, whose volume is assumed to be unit. 
4) There are N/2 elementary cells in all. 
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where v is the velocity of the electrons corresponding to the Fermi surface, 0 is the 
angle between h and h’, and the integration is carried out on the Fermi surface with 
respect to the solid angles coresponding to h and h®. 

For thermal vibrations we shall use the expressions 


l ; 1l ifr 
w=) vf) etigh) ==) milfs) eats) (10) 
yN fs yN fs 
where v, and V; are the vectors corresponding to the a- and ĝ -sub-lattices which give the 
direction of vibration for the propagation vector wave f and the vibration mode s®, 
while q(f, s) are the vibration amplitudes corresponding to f and s, assumed to be 


harmonically dependent on the time in the form exp (+ iwt), where œw = o(f, s). Thus, 
the three terms of the perturbation hamiltonian become 


oe 7 YY oaWar — 1) tri q (f, 9) — 
LS, cd 
a za Y Y vW — ry) cit g (Ff, 5) 
fs k 


moy ATena e (11) 
al 
m = — oor? i W — J eur; 9 S 


ait z 2 z (1 — ur) vaW (r — r) eitrk q (f, s). 


2. The calculation of pertubation matrix elements 


Let us denote with D the vector connecting two neighbours of the first order, 
and with D the magnitude of D, which is evidently constant irrespective of the direc- 
tion of D. We shall introduce the function 


YR (r) = Vi eikr yD m=/F ehr u(r + D) (12) 


which, owing to certain conditions of the type (5), depends only on the magnitude 
D of the vector D. With the aid of this function it will be possible to express certain 
integrals which appear in the matrix elements of the perturbation hamiltonian. We 


5) In deriving the formula (9) it was assumed that the energy bands are spherical, a hypothesis 
which was adopted also by us. The original formula of Gibson (1956) contains a misprint, which was 
corected in (9). 

6) There are 6 modes of vibration, of which 3 are acoustic and 3 are optical modes. It should be 
noted that the relations (10) are valid for perfectly ordered lattices. We shall assume that the relations 
(10) are also valid for deviations from order but in this case the mass of an atom in a or B, respectively, 
will depend on the order parameter. See formula (23’). 


im {re m 
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shall be employing the following notations for the integrals appearing in the calcula- 
tions: 


oe PR) Ua) cea Woo (As (rn) WA ir) Py (r) dr 
V Vv 

US [PRO UC) Vea de WA PEWA wR lr) dr (13) 
V yV 

E O O E WE = ROWE Y, (vr) de 
V V 


UBD — [ WP (r) Us(r) PR) de Wie = J ER (e) WP (r) YR (W) dr 


V V 
and 
Ds th Ee AO a e ie! AD a S BO AO 
U aw U Uhn) Unw = 2 (Unw — Uhn) Wiw Az Wir Sa Wir 
Dee BD ‘AD 
Wir a? [Win a Wir (14) 


‘the integrals being taken over the entire volume V of the alloy, which is equal to 


1 NQ, Q being the volume of an elementary cell. These integrals can generally be 
reduced to the integration over an elementary cell, whose center is at the origin of 
coordinates, if the potential decreases rapidly enough at the exterior of the cell. In 
this case, we can write, e.g.: 


f Piir) Ur) Yy (r) dr = f P(r) U(r) Y, (r) dr 
Vv Q 


and so on, for all the integrals. 
We can now calculate the matrix elements Hf; given by 


H3S = f P(r) (H° + H° + H”) V(r) de = Why + Hi + Hh (5) 
J 


The calculation is rather simple and follows the same lines as the one given in a previous 
paper (Corciovei 1959 b). It uses the periodicity properties of U(r). and the formulas 
(13). By introducing the vector K = h’ — hand denoting r; simply by j, and r, by k, 
we obtain i 


eee = N- (Z > v, W4 eK S AE 2 pa DAW ete) q (fs) 
Sar, 9 
Hi, = 3) 0 — h) Uw e — D3 — 1) Ui u 
j k 


fs J 
Formulas (16) contain scalar products of the vW type, where v stands for V, or Vp, 
while W is one of the vectors in (13). We shall use for an approximate relation currently 


He, = N+ [- EEL — 1) 0 Why E + DDL — ty) OW EN] g (f, s) 
, fi 
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used in conductivity problems. This relation can easily be justified in the present 


case, too. Indeed, we can write W,,, in the form”: 
2 p 
Way = F jl eikr u* (r) u(r) grad U(r) dr (16') 
ý 


Carrying out the integration by parts, we may write 


War = =a grad (e'Er u* (r) u(r)) U(r) dt 


the other term vanishing owing to the conditions imposed at the boundary of the 
integration domain. If we limit ourselves to the integration over a single elementary 
cell centered at the origin of coordinates, owing to the rapid decrease of U(r) outside 
the cell, we may write 


War = sa eKr u* (r) u(r) U(r) dt — — al eit grad [u* (r) u(r)] U(r) dr 


Q 


Assuming that the quantities [u* (r). u (r)] and U(r) are even with respect to 
each of the three coordinates and noting that exp (iKr) is practically constant within 
the elementary cell, the second integral vanishes. Thus we obtain the simple relation 


Wy = —iK | Y U(r) Py de 
2 


We can, therefore, use the following approximative relations 
W49 =—iKU?2 W22=-—iKU32. W%,—— iKU}, 
We = —iKUR? «Wie = — iKURD OW? = — iKU?, (17) 
It should also be noted that b exp [i (K + f)j] differs from zero only if K + f = 0, 


j 
when it is equal to NV/2, and the same is true for the sum taken over k. If we denote 


ie = —Kandp,=v,K p= vK and using (17) we get, instead of (16) 
= + [i N/2 N*] 2 ips Dak + Pg 4 wl qa (fo; s) 
Hir = 2 (= H) Uire = (1 — wy) Uks ann (18) 


Hm, = iN [SB OH) Pa Uiw EPI — D) D) (1 — pd py UR PMY (fs) 
so j fs k 


To calculate the resistivity, we must figure out the square of the matrix elements 
of the perturbation. At the same time we must take the average in time. If we denote 


7) When functions yd are present, function u(r) is replaced by uP (r) and the arguments remain 
unchanged, 
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this operation with ~ and if due account is taken of the dependence of the amplitudes 
q(f, s) on the time, we shall have evidently 


q(f, s) =0 (19) 
and 
a(f.s) (P's s') =A(H. 5) 5,76," ACh 5) =|a(f 5)? (20) 


For the sake of brevity, we shall not discuss in detail this averaging and in the 
following we shall use the relations 


rrd 
Taw a [Ain = Thaw’ aR Taw a6 Taw (21) 
where 
rrd d 
Tir = ees Taw = JAENA 


Tre = (Hgy Hay + Her Hiw) + Hgy Hie + Hir Hav) + 
E E TE pamm a 
+ (HA, He, + HiH) zl er) 


Owing to the relation (19) in Třw >» after the average is taken, only the second and 
the fourth terms are left, since the others vanish. Moreover, since Hyw » contains only 
the term corresponding to f, in the sum over f, it will be seen that in the second term 


of T%, also, only the term corresponding to fy = — K will be left after taking the 
average in time. Therefore, we may write 
ee 


hh’ — 4 X [Pe ea AF Peg URDI [Pa Us, ai Pp UEP i A (f; s) 
Ta e U N D A Lp) [Oh Uae BP + 
PREN J 
SE Uiw UD eEG) ae D D (1 E Hp) (ie Uy) k beac (22) 
k R 
N 
Tg = (L—R) D (Pa Uiw — Pe Um) (Pa Une + Pg Uae) + 
| AO BD\* Tm i2 
zt: (Pa Uik — Pp UP ) (Pa Uhr + Pp Urn) SAT. s) H |H | 
where |H”,,|? has the form 
r 3 {i 
a"? = : DY I5 $ a-u- a) et SF eC ae) 
s hE 
—>) > —#,) (l— u) INC) (og Un) eae ts US \epsl) en dal 
; k 


J 
+ >) 0 =e) (L — uw) (Ppa)? eH Nek) 4} (ff, 5) (22’) 
k F 
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and the factor (1 — R) appears as a consequence of the adoption of the long range 


order parameter, which is defined as 
2 2 
Rp ATNA M (23) 
j 


3. Calculation of the resistivity 


Using formulas (9) and (21) we may write 


f (Thy BY pee + Thr) (1 — cos 8) dOn dQ): = 0o + Od + Om 
(24) 


e= het : aye 


where the integration is performed for the solid angles corresponding to h and h, 
with the assumption that h? = h’2, their common value being given by the value 
of Fermi energy. In the following we shall assume that Up» is not dependent on the 
angle between h and h’. Since the vectors h and h’ in (24) have a fixed magnitude, 
the quantities Up,,and UP, may be regarded as constant. Consequently, the integrals 


originating form Tj, and Tm, — |H?%,.|? are of the form 
> ff PaPe A (for s) (1 — cos8) dQ, dQ, (25) 


where z relates to the contribution of the 6 vibration modes. In the case of the inte 


grals originating from T4,, and |H Hej certain exponentials appear under the integral 
sign; such cases, however, will ba treated in Appendices 1 and 2. 

The quantities p,= V, (h’ —h) and pz=v,(h' —h)contain the scalar products of the 
vector h’—h and the vibration vector V, or Vg. The vibration vector for a determined f 
has a determined direction for each of the three optical modes as well as for each of the 
three acoustic modes, namely: a direction parallel to f (longitudinal vibration) and two 
directions perpendicular to f (transversal vibrations) for both the optical and the 
acoustic cases®. However, in the present case, 


| f=h=-K=—W—h 


which means that — fọ has the same direction as h’ — h. Thus, the sum in s will 
contain only the contribution of the longitudinal vibrations, and we shall be able to 
write p, = VK and pẹ = v,K the contribution of the transversal vibrations vanishing 
because the scalar product of V, or Vg and h’ — h is zero. 

As for the term A (fy, s), if the significance of this term resulting from (20) 
is taken into account, we may write, aine to the quantum theory of thermal 


8) It should be noted that, in fact, parallelism or perpendicularity are HEE obtained only 
in certain particular cases. 
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vibrations: 


1 w 
5 Mao?(f, s) vz A(f, s) = -5 Mpo? (f, s) 3 A(f, s) = : oe 


where M, and M, represent the average mass of a point in the a- or f-lattice, respec- 
tively. It is evident 


that M,=}4(1+R8)M,+40—R8)™, (23’) 


where M, and Mp are the masses of the A- and B-atoms. In the case of complete 
order (R = 1), we have M, = M4, Mẹ = Mpg while in the case of complete disorder 
(R= 0), we have M, = M, = 4 (Ma + Mp). 

By the foregoing considerations, the integrals (25) have the form 


h oe dQ, dQy 
JMM; > ffe — cos v) Dida (eho — = i (25') 


where X relates to w (fo, s), which depends on s; indeed, the sum in s has two 
s 


longitudinal terms corresponding to the optical and acoustic vibrations. In a first 
approximation, the frequency of the optical vibration may be regarded as a constant 
and equal to a numerical value wọ, which does not depend on f and s. For the acoustic 
vibrations, we may assume a linear dependence of w on f; therefore w = cf, where 


~ cis the wave propagation velocity. Consequently, if the terms of the sum are written 


out, integral (25°) becomes: 


ET fj i K2(1 — cos ð) E FT) +o (eo 7 | dQ, dw 


(25°) 


In calculating an integral of this type (25”) we shall apply current methods (Gibson 
1956), as the scalar quantity K defined by the relation 


K? = (k — h)? = k? (1 — cosd) 


If the z-axis is considered as directed along h, we have dQ, = 27 sind dò = 
=2nK dK/h? and (1 — cos) = K?/2h?. Substituting these in (25’’), it is seen that 
no term dependent of h is left and we can, therefore, integrate by der which yields 
the factor 4r. INE the integral (25”) becomes: 


Anèh ia 1 i 1 S 
VMM "| ar ik aero | y° dy (25°) 


with the substitution Ae mys 
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We can now write the expressions for the resistivity” 
pee ree 32/3 
’ E Ən i Fe i eA | —gecle —ae + 
sv A \ he? v? (400)?] LYM, «=~ Me Wy (efolkT — 1) 


2 
1 y* dy 
ch (cfichy/kT — 1) 
0 


04 = [3x/he® v? (dox)?] 8x2 N[(1 — R?) (UP, + UR) — 27 (0 — R3) U3, UR 
AE ħe? v? (4n)? 2 VM, Lee VM, 


Uia 32/3 a se 
a VM; mia w PEO at (eħchvikT — 1) + en (20) 


where @,, is the contribution of (Hm 2, which, as shown in the appendix I, is negligible, 
while for calculating @, we have made use of an approximate expression of |H Aa, 3 
that is of 7¢,,, taken from another paper (Corciovei 1959) (for certain useful details 
see Appendix II). 

It is important to see what is the sign of ọ„. To this purpose we must consider 
that Uá, and UB, are approximately proportional to the masses M ‘4 and Mpg, respec- 
tively’). Here the sign present no interest, since we have to do only with such products 
as Ufa. Ue or ya, or UB We may, therefore, write 


Une  Unw _ 1 Une Ume al Une Uwe | Ms — My _ 


VM, VM, 2 YM, ~2 YM, YM, 


Mz — Ma 1 
-E ~ VM- VMs) VM+ VM. a A 
VM; ~ VM. VM; 
In the case of a more advanced order we have M,~ Ma and M, ~ ~ Mgand it 


is seen that (VM, — VM M4)? intervenes, which shows that ọ„ is positive. The positive 
sign persists for RA 1. 


®) With the hypotesis that the quantities U are real (spherical symmetry of the potential U ‘A> Up). 


In the complex case the alterations required are evident. The expressions (26) seem to be proportional ~_ 


to N. In reality, if the normalization of the Bloch functions (formula (6), see also (16’)) is considered, 
also N-? from the matrix elements intervenes and the final formulas are proportional to N-1. 

10) This is so because the potentials are approximately proportional to the nuclear charges, which 
in turn are proportional to the atomic masses. A more complete discussion is possible (Corciovei 1959b). 
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Using the expressions for ọ,, 03 and @,,, we may calculate the relative quantities 
0,,/0, and 0, /o,!Y. We obtain: 


S M ES UA 
Caoa e L ieee r ER) (27) 
Ov MaU aM bee 


l 


= -z [Ma *U nw — Mp Uw] [Ma Uin + Mp Un] x 
d 
2 
32/3 il yt dy 3 
k Cae) T il ehichy|kT — + h h? (1 — R) (28) 


0 


xX 2 2 ~ 
2 [(1 — R?) (Uiw + Uiw) — 2y (0 — R?) Unw Unn] 


In the case of high temperatures, Ach < kT, we may expand in series the exponen- 
tials exp (Achy/kT) and the ratio o„/ọq takes the form 


Om 1p A 4h (1 — R) kT 
Bag Eh (1 — R%) (US + Udy) — 2y (a — R®) Uw Un] 
(28) 


wu a 


where we have denoted with F the quantity 


T p i ra Le ns 
y Ma VM; VM. VM; 


In the case of low temperatures, we have 


2ħch|kT 


2 , 

1 yeah TIS : z+dz 

ch | eħchyikT—] ch \ħch e— l 
0 


0 


with the substitution hchy/kT =< The maximum of therfunction zi (exp ae 1) 


is near z) = 4, after which value the function decreases rapidly because of the exponen- 
tial. However, the upper limit of the integral 2 ħch|kT > 4, since hch > kT and we 
may take oo for the upper limit. 

The integral now equals I’ (4) ¢ (4), and we may write 


2a eels Jawa-R 
=e se ea ats) SOLON 8") 
ea 4 [A — R (Ufy + Udy) — 2y (0 — R?) Unw Un] 


where I is the Euler function and ¢ the Riemann function. These functions have been 


taken for the particular value 4. 


11) The ratio ọ„/0q results simply by dividing Q,,/0g bY Qn/@y- 
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Me EE ee 


The formulas (27) and (28) enable us to say something about the importance 
of the three resistivity terms, 0,, Og» @,, in the various ranges of temperature. Thus, 
for very low temperatures, the only important term is ọ4, the vibration term and the 
mixed term being very small as compared with it. This is due to the integrals in (28), 
which. decrease at least as the fifth power of the temperature at low temperatures. 
Thus, in this range 0, > 0, > @,,- At somewhat higher temperatures, o, and @,, incre- 
ase more rapidly than gj. 

Above the Curie temperature, @, remains practically constant”, while ọ, becomes 
predominant. The ratio @,,/o, remains very small, especially owing to the difference 
in the numerator of (27), while the ratio @,,/0, continues to increase with temperature 
also above the Curie temperature. 


APPENDIX I. Calculation of o, 


We may write 


: 3n Print 
Om = he®v? (An)?, |Hrr | (1 == GOS 8) dQ, dQy (I.1) 


where | HT, Hm i is given by (22’). 
With the same assumption as those made for formula (25), the integrals involved 
will be of the form 


I=) f f: Pq Bp ORES A(f, J dQ, dQ, (1.2) 
fs : 


We will also have 


2 h/Ma : 
vz A(f, s) = Sogeer oy (1.3) 


At lower temperatures, v? A (f, s) will have smaller values, since the exponential 
will be considerable. At higher temperatures, œ < kT and we may use the expressions 


v2 A (f, s) = kT/M, w? (1.4) 


If this expression is maintained over the entire temperature range, it will lead to an 
overestimation of the integral, should the temperature be low. If for optical frequencies 
we assume that wœ = wọ (constant), and for acoustic frequencies!®) w? = c? f2, then 
the term (I.4) is taken outside the integral, since here f is independent of h or W. 


12) The temperature dependence of Qq is investigated in a previous pepe by one of the authors 
(Corciovei 1959b). 


18) We neglect the dependence of c on s. - 


ER a e 


The Mathiessen Rule in Disordered Binary Alloys 659 


Therefore, we may write: 


kT T ; 
= TM, a? om f| Teper vp] EG-K) (1 — cos 8) dQ, dQy = 
@+'4B fs 
kT Ip = ia 
SENS — ORD J,, (¢) = EN Juno I.5 
VM. M; 0 ix (9) /M,™M; ” a(o) ea 


For calculating the integral J;,, whose significance results from formula (I.5), 
we cannot operate as before with the vector K, since now there occur scalar values 
of h and h’ with v, or j—Kk. In order to calculate the integrals, we shall adopt spherical 
coordinates, chosen as follows with respect to the Cartesian axis: the z-axis is directed 
along the direction of j—k, and the x- and y-axes located in a plane perpendicular to z. 
In this case, we may introduce the angles with respect to the z-axis and the angle of 
the projection of the vector considered on the xy-plane with respect to the x-axis, 
which we shall denote with g. If v, is the direction defined by the angles y and Qy, 
Vp is that defined by y’ and gp, h by y and @,, h' by y’ and Py» then, by making use 
of the spherical relations, we may write the relevant scalar products in the form 


V, h/v, = h [cos y cosy + sin ysin yp cos (p, — Po)l 
v,h'/v, = h [cos y cos yp’ + sin y sin yp’ cos (Py — Po)] 
Vghħ/vg = h [cos y" cos y + sin y’ sin y cos (p, — pol 
Vs hiv, = h [cos y’ cos y’ + sin y’ sin y cos (Py — Po) 
hh'/hh' = cos 0 = cos y cos y’ + sin y sin y’ cos (py — Ph) 
dQ, = sinypdyd p, 
dQ, = siny dy d oy 
K (j—k) = h [r; — 7] (cos y — cos y’) = AR;, (cos y — cos y^) (L6) 
By integrating over y, and ¢,, from 0 to 27, we have f 
Jy = 20h [A m (Lh +12 + Lh + +2 (312 —4hh—2hh +1 + 2) 
(1.7) 
where m = cos y. cosy’, n = sin y sin y’ cos (Po — Po) 
+1 


and l, = f x" exp (thR;,x) dx 
=i 
Assuming that V, has the same direction as Cy the average of mis 1/3 and that 
of n is 2/3; these averages have been taken for all possible directions of V, which was 
assumed to be isotropically distributed (for various f). . 


14) s takes three values (for optical and acoustic modes respectively), 

15) We shall assume that the longitudinal vibrations with respect to f play an important part, while 
the transversal vibrations havé no significant contribution to make. In this case, we may assume alles If 
These approximations are permissible since we are interested only in the order of magnitude of Q,- 
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We get!) 
Jy =i m h A + 41? + 31? — 24, hl (1.8) 


and, hence, we may write 


en = ara N a pz Da (1 — m) l — m) Six Ji rA = < 
a2) on Bhs eee Tce (L9) 


ep È (L — wx) (L — ur) See Jie ae 


where, as per formula (1.5), Spp is i 


rea ‘ 
Se= 2, & +5 zp) eif J-W) sig €10) 4 


Since we have obtained S;, from a product of the form |H7, Hi" | it follows that 
Sik will always appear too, and we may write, instead of formula (1.10) 


Se » (a $ ap) cos [f(j — k)] | (ul) 


We shall substitute to X a volume integral in the space f carried out over a spherical _ 


domain of radius x/a. Therefore 


cos[f (j — k)] dV = 


EA at eae 
s-h | (areal |; 


< 
g 


Ar OR 
E ER 
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spi 
On the other hand, if we calculate the integrals of the type J, = (ie exp (tAR,, x)dx, 
1— 
we have for large R the expression 


2x? h? 32 cos(2h Rjx) 11 
3 (Ry? a 


L 
Jir ~ 


It is obvious that the series X X S; J} is covergent and, therefore, it can by no 

j k 
means be proportional to N?, and is only proportional to N. This is due to the fact 
that for a given atom, we make firstly the summation over all the shell of neighbours, 
and the contribution of very distant neighbours is negligible. With more reason our 


series XZ (1 — u;) (L — Mg) Sip Iie which is the same series as above with the difference 
jk 


that a number of terms are missing, is at most proportional to I. It follows that e Om 
is of the order of 1/N, and therefore is negligible with respect to o,- 


APPENDIX II. Some integrals appearing in 03. 
Making use of a previous paper by one of the authors (Corciovei 1959), we find 


Ai 


for |H¢%,,|2 the following approximate expression: 


ryd 2 N 2 0° D? 2) 770 D 
where we have put 


i 
ee y Cam RS) 


where z represents The number of first order neighbours of an atom, and the sum 

applies to all its first order neighbours, D being the vector connecting the atom with 

its neighbours. š 
Corresponding to the relation (24), we have 


et e 9) dQ dQ, (11.3) 
Ce he® v2 (4n)? | hh’| ( — cos h h . 


a 
The only part depending on h, h’ in the expression of |Hiy/|? is yg- It follows that we 
have to do with integrals of the form 


f i) (1 — cosb) dQ, dQ, = 16x? . (11.4) 
and an integral of the form 


l6n*y = z ff Yg {1 — cos 8) dQ, dQy = > fel eKP (1 — cos 9) dQ, dw 
Š (11.5) 
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a 


where the integration was carried out as shown in Appendix I, with the difference 
that the z-axis is here directed along D. In the present case, 


+1 
= f x" exp (ih Da) dx 
— x 


As a result of the calculation, we obtain 


cos(2hD) sin(2hD) sin? (hD) 
r= |— age t T e aad 


where D = (V3/2)a may be substituted in the equation. For example with h = (V3/2) x 
X (xja) the approximate relation y = 0.5 is obtained. When h decreases, y increases. 
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DYNAMICS OF THE UNIVERSE 


By JAROSLAV PACHNER 
Physical Institute, Technical University, Praha, CSSR 


(Received March 23, 1960) 


The dynamical law of the universe is deduced from three simple and plausible 
principles without using the energy momentum tensor of continuous matter. — It is 
showed that the square of the Hubble’s factor may be interpreted as proportional to 
a certain density of matter which is transformed during the expansion of the universe 
into rest mass (in form of the evolution of new nebulae and stars) and into different forms 
of energy. The maximal radius of the universe is reached when all Hubble s mass is con- 
sumed. During the subsequent contraction the stars and nebulae are destroyed by chang- 
ing their rest mass into Hubble’s mass. In the moment of maximal contraction all the 
mass exists only in the form of the Hubble mass. This process repeats itself from infinity 
to infinity. The computed numerical values of the radii of the universe, of its total mass, 
and of its age are of the plausible orders of magnitude. — It is also proved that the 
proper meaning of the well-known cosmological constant is: density of the total matter 
kept on a constant value because of the continuous creation of matter from nothing. 
An example of such a universe with continuous creation of matter is the de Sitter universe. 


Introduction 


Some years ago the present author made his first attempt to clarify the 
well-known cosmological problem, but his investigation remained then unfinished. 
Inspired by the interesting paper of Lindquist and Wheeler [1] on the lattice universe, 
he resumed recently his work and obtained simple and plausible results on the dynam- 
ics of the universe, which proved, however, to be different from those of many 


previous papers. 


l. Principles of dynamics of the universe 


The dynamics of the universe may be founded upon three plausible principles: 

(1) The geometry of the universe is Riemannian. On the basis of purely geometri- 
cal considerations the spatial curvature of the universe K; has to be determined. 
The universe is assumed to be spatially isotropic; this assumption, based on astronomi- 
cal observations, is the necessary starting point of our considerations. 


(663). 
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(2) The universe is divided into N equal cells (with N finite or tending to infinity). 
The total mass within the cell causes a certain curvature of the space, the magnitude 
of which on its boundary has the value K. The curvature of the space has to be deter- 
mined by the laws of general relativity theory. 

(3) The curvature on the boundary of the cell equals the mean curvature of the 
universe: 


e (1.1) 


This is the dynamical law of the universe expressing the causal connection between 
the mass as the physical cause, and the geometrical curvature of the Riemannian space 
as its consequence. 

It seems perhaps odd to determine the curvature of the universe from the curva- 
ture K caused by the mass of only one cell, but in fact the integral effect of the mass 
of the remaining (N — 1) cells on the curvature is implicitly contained in the geometri- 
cally determined mean spatial curvature K, of the universe. 

We quote here the formulae for the curvature in the direction of the coordinate 
surface (x, ,) written under the assumption of a diagonal fundamental tensor [2]*) 


Rw, yv 
? 
Suu Evy 


lL [ote oe a\(a a)fa 
Ri T ce =e Tweed NaS a. 
ti 2 ( 9x, IXu ) Š (he tot a w) E 


By K; we denote the spatial curvature of the spatially isotropic universe, i.e., the 
curvature computed in the direction of the coordinate surfaces (r, 8), (r, y), and (8, p), 
and by Kj, the curvature calculated in the direction of the coordinate surfaces (r, t), 


(ð, t), and (p, t). 


K (dg X) = — (1.2) 


where 


2. Geometry of the universe 


We consider first the spherical universe with the metric 


R?(t) f dr? ; 
ds? = — T ee T r? (d8? + sin? ð dø?) | + dt?. (2.1) 
Ro 
By the transformation 
r= R Sin x, ER (2.2) 
where Rọ is a constant and using the substitution 
R OSR AG); : (2.3) 


1) The minus sign in formula (1.2) is added, because we use the metric (—,—,—,+), whereas 
Rashevskir [2] assumes (+,+,+,—). 
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we have 
ds? = — R? A? (dy? + sin?y (d 8? + sin?) d y)] + R? d t°, (2.4) 
where 


Ons ye T Osen U S2z. (2.5) 


The curvatures K,, K,,; are now given by the relations 


1 WANA 1 deed 2A 
Se l(a a a Ku= REA a?’ i) 


The spatial volume of one cell (strictly speaking: the hypersurface of one cell) is deter- 
mined by the formula 


X n 2n 
a fel [V — det g,,¢ 440d 
00 0 
= n R} 4? (x) (2 % — sin 2 7) Gian 253) (2.7) 
The volume of the whole universe (y = 2) is then 
V = 27? R? A? (T). (2.8) 


In the pseudospherical space the metric is as follows 


R? (t) dr? : ; 
alee 2 (dg2 2 2 2 
st E Fee ga +'sin? } dy?) | + dt?. (2.9) 
Ro 
With the help of the relations 
r = R, sinh y, ERAL (2.10) 
R ( = R, 4 (1), (2.11) 
we get 
ds? = — R? A? [d y? + sinh? y (d 9? + sin? Od p) +. Rp dz, (2.12) 
where | 
0Syso, OS0Sa, O59 82H. (2.13) 
The curvatures are now ; Dr 
LAPA ioe 
So Ria E ie Baten = —— —— 2.14 
Ky = R? pe a A f Kyi RRA ace ( ) 
and the spatial volume of one cell 
v = n R? 4? (1) (sinh 2 y — 2 y) 3 (2.15) 


The volume of the whole universe is, of course, infinite. 
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Each of both spaces contains two special cases. In the spherical space we have 
the Einstein universe of constant radius (i.e., with A = 1), and the space of constant 
curvature K; = K,, = 1/R? (with A = cosh t). The pseudospherical space contains 
as a limiting case the pseudoeuclidean space (in the very late stages of expansion 
when A->oo and d A/dt > 1), and the de Sitter universe of constant curvature Ky = 
= Kı = 1/R? (with A = sinht), as may be easily proved by the transformation 


p= — ae, ee rs (2.16) 


with the help of which we obtain from Eq. (2.12) the metric of de Sitter’s universe 
in the known form [3] 


ds? = — R? {dw? + sin? w [d¢2 + sin?¢ (dd? + sin? bd ¢?)}}. (2.17) 


3. Physics within the cell 


Lindquist and Wheeler [1] showed that also with a relatively low number 
of cells, N = 600, of their lattice universe the radius of the universe calculated by 
their method agrees to 1.2% with the well-known result of Friedman [4] deduced 
under the assumption of a perfect spatial isotropy of the universe. We may be thus 
sure that it is possible to choose sufficiently large regions of the universe on the 
boundary of which the gravitational field is with a sufficient accuracy spherically 
symmetrical. This field is then the same as if the total mass within the cell would 
be concentrated in one point in the centre. On the boundary of the cell we may 
use Schwarzschild’s metric. 

Each of our cells is not isotropic, but merely spherically symmetric. If we wish 
to state a causal connection between the curvature caused by the mass within the cell 
and the mean curvature of the universe, the choice of the curvature K in the direction 
of the coordinate surface (®, p) as determining the spatial curvature K, of the uni- 
verse seems to be the most (and the only) natural one. 

The curvature on the boundary of each of the equal cells computed by Eqs. (1.2) 
and (1.3) is then given by the formula 


K = 2mjr’ (3.1) 


where m denotes the total mass within the cell which causes the curvature of the space. 
By r we understand the radius of the cell equaling R,A sin x. 
In the spherical space, we put Eq. (2.7) into (3.]) and find 


K = (8x P/3) o, (3.2) 
o being the density of the total mass and P a numerical factor given by the formula — 


3 27 — sin 2% 
RIN) 7.4 a (3.3) 
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The angle y corresponding to the radius of the cell is determined by the condi- 
tion that the volume of the universe is divided into N cells: 


N= Vp. (3.4) 
The values of P(N) for N = 5, 8, 16, 24, 120, and 600 can be easily computed with 
the help of Table IV of paper [1]. If N tends to infinity, y approaches to zero, and 
from Eq. (3.3) we have 

P(oo) = 1. (3.5) 


We must now make an assumption concerning the density of the total mass 
in the universe. We write 


Q(t) = cı 45/4? (x) + IT (2). (3.6) 
Inserting this expression into Eq. (2.8), we determine the total mass of the universe 
M = 27? R? [A o, + A? (x) I (v)]. (3.7) 

Since 
A, = Afr), (3.8) 


t, being an arbitrary moment, ọ} denotes the mean density corresponding to all 
the forms of matter which cause the curvature of the space in the same moment 1. 

The total mass in the universe remains unchanged as long as JJ (r) = 0. This 
function describes the creation of matter from nothing if Z(t) > 0, and its annihila- 
tion into nothing (or the creation of a negative mass) if I(t) < 0. We introduce 
this function in order to recognize later the proper meaning of the cosmological 
constant. Considering for a moment ọ, = 0, and letting 


I(t) = A/8xP = const., (3.9) 


we see from Eq. (3.6) that the density of the total mass remains now constant because 
of the continuous creation of positive (A > 0) or negative (A < 0) mass. 


4, Dynamics of the universe 


The dynamical law of the spherical’ universe will be obtained by inserting Eqs. 
(2.6) and (3.2) with (3.6) into (1.1): 


RERAN ad 8a lea? | re, 
(44) +a] learo ore 
We restrict our investigation to the function J (t) given by Eqs. (3.9) and find 


A 
C-A+ ZARA 


(4.2) 
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with the dimensionless constant 
C = (82P/3) ale Ro 1: (4.3) 


Eq. (4.2) in a slightly different notation is the same as Eq. (7) deduced by Friedman [4]. 
The constant 2 from Eq. (3.9) is thus the well-known cosmological constant. 
Its real physical meaning is: density of the total mass in the universe kept on a con- 
stant value because of the continuous creation of matter from nothing. Since the 
dimension of a density in geometrical units is (length)~*, the proper meaning of the 
cosmological constant remained hidden for so long a time. 
In the universes of constant curvature with 


A =cosh t,, (= sinh: 7, (4.4) 


respectively, matter should be created in such a degree that the total density remains 
exactly constant, as follows from Eqs. (1.1) and (3.2). In the beginning of the expan- 
sion, the amount of the mass in the former is finite and different from zero, and in 
the latter (i.e., in de Sitter’s universe) it is zero. The radii of both universes grow 
to infinity. The stationarity of de Sitter’s universe is but a fiction caused by a suitable 
transformation of coordinates (see Eqs. (2.16) and (2.17)). 

More than 24 centuries ago Parmenides expressed the view that all that exists 
has never been created out of nothing and can never be ulterly annihilated. Up to 
the present the hypothesis of creation or annihilation of matter is strange to any 
scientific theory, and we exclude therefore the function JI (t) and the cosmological 
constant A from our considerations. 

The integral in Eq. (4.2) with A = 0 may be evaluated by elementary methods: 


t —42=VA(C— A) + CarccosV A/C (4.5a) 


The integration constant has been so chosen that the time is measured from the begin- 
ning of the expansion. The parametrical representation of Eq. (4.5a) is 


A = 4 C (1+ cos £) = C cos? 4 &, (4.5b) 
r — 4n = 4 C (E+ sin 4). 


In analytical geometry the function (4.5a, b) describes a cycloid. 
Let us now revert to Eq. (4.1) without the physically meaningless function JT (r). 
At first, we make the important assumption that the mean total density 9, in the 
moment T, which causes the curvature of the space is composed of two parts: 


0: = On + On (4.6) 


The first of them, og, contains the density of rest-mass and of other forms of matter 
(e.g., the density corresponding to the thermal energy, radiation ect.), while the other, — 
Cp corresponds to the kinetic energy of the cosmical expansion and may be therefore - 


called the Hubble’s density. 
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The first term on the left-had side of Eq. (4.1) is the square of Hubble’s factor 


of cosmical expansion h, for 


A (2 aA) (188) 9 
agile Rh dies (4.7) 


This term cancels in Eq. (4.1) the above mentioned Hubble’s density given by the 
relation 


07, = (3/87 P) h?. (4.8) 
The radius of the universe at the present time is then 
R, = V3/8x Pop. (4.9) 
From Eq. (4.7) we determine also 
dA/dt = R,h. 
With the help of Eqs. (4.5b), we find 
tet = + Venler (4.10) 
and further the maximal radius of the universe 
3/87 P 
er e ea EE Oa (4.11) 
eR Vor . R 


From this equation the density in the moment of maximal expansion can be easily 
computed to 


2 
o= oR ( ET 5 ` ; (4.12) 


The total mass in the universe is constant 


Vor OR 4P 

During the expansion there exists a steady transformation of Hubble’s mass into 
rest mass (in form of the evolution of new nebulae and stars) and into other forms 
of energy. In the moment of maximal expansion Einstein’s universe with zero Hubble’s 
mass is reached. After that a subsequent contraction begins, during which the stars 
and nebulae are destructed by changing their rest mass into Hubble’s mass. Also 
other forms of energy are transformed into Hubble’s mass, till in the moment of 
full contraction all mass exists only in the form of Hubble’s mass. Immediately Hubble’s 
mass begins to transform again into rest mass and other forms of matter. This process 
repeats itself from infinity to infinity. The red shift of spectral lines emitted by the nebu- 
lae testifies that we are, living in a period of evolution of the rest mass. This agrees 
with the astronomical observations of new nebulae and stars by Ambartsumian. 


- i 
M = 2n? R? o, = 27? R3 0) = 217/128 P* on + OH _ 3% Re (4.13) 
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We now determine the time of the half-period of pulsations of the universe, 
i.e., from its maximal contraction to maximal expansion. Since the constant C in 


Eqs. (4.5a, b) equals unity, this half-period Tọ is 


Te ie E oe (4,14) 
Ver. QR 


In accordance with Eq. (4.10), the time that will elapse from the present days to the 
moment of maximal expansion is given by the formula 


V or OR OR OR 


In our calculations geometrical units have been used throughout. We now trans- 
cribe the resultant relations into the usual CGS-system (G being the Newtonian 
gravitational constant, and c the velocity of light): 


Op = (3/8 PG): k = LB Oh [ojem?] (4.8a) 
cV3/8PG  4.0x 1018 4.2 x 10-5 
= = n CMIS S 
Vor V Por V Por 
cV/3/8xPG or + oH _ 4.0x 101 or + on 
Vor eno Poe Por 
42x10” or + Qr t 0H 
V Por OR 
_ ÈV2ZIa128P* G? or + en _ 11x10% a or + on 
Vor QR VP? or 
_ 0.55 x 108 or + on 
s VBR on — ~8R 
r, — V3a]2PC or + on _ 21x 10° or + on 
Vor QR V Por OR 
6.7% 10 or T OH 
qy V Pre. oe 


__ V3/8x% PG y= 4 ORT OH OR + OH Ae V2) ai 
Eora OR OR OR 
2 1.34 x 103 OH R+ On 
as Ys +O STi S fel = 
4.25 x105 E z H a | 


RIE [light-years] (4.9a) 


Ro = [cm] = 
[light-years] (4.11a) 
[g] = 

Mo (4.13a) 


[sec] = 


[years] (4. 14a) 
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Assuming @p = 5 X 107% g/cm3, k= 06 x 1078 see=1, P=1 we obtain 
peo 10 lems) 0) = 9:4 x 10-2 g/cm, 
eon e ie TO sc OES 
M =3.7 x 10g = 19 x 10°Mo, 
T, = 22 x 10° years, T #3 19 x 108 years, 


The computed values are of a plausible order of magnitude. Only the time (T) — T*) 
that elapsed from the beginning of the expansion to our days is rather short, probably 
because of the too high assumed Hubble’s density oy. 

Let us refer once more to Eq. (4.1). If instead of Eq. (4.6), we assumed here 
that Hubble’s mass did not produce the curvature of space and we had 


01 = Or, (4.16) 


the amount of rest mass and of other forms of matter would remain constant during 
expansion and contraction of the universe and Hubble’s mass would be stored during 
the expansion as potential energy of gravitating masses receding from each other 
in order to be transformed again into kinetic energy during the subsequent contrac- 
tion. This is, however, an idea of Newtonian mechanics which contradicts the prin- 
ciples of relativity theory. 

A first inspection of Eq. (2.14) shows now that the pseudospherical space without 
steady creation of matter could be formed either by a negative mass causing negative 
curvature of space, or under the assumption that in the beginning of the expansion 
the density gy is greater than gp and do not influence the curvature of space. Because 
there are no empirical proofs for a real existence of such a negative mass and relati- 
vity theory rejects the idea of a potential energy, the pseudospherical space as well 
as the pseudoeuclidean one as its limiting case are excluded from our investigation. 


5. Dynamics of the universe by the method of Lindquist and Wheeler 


We give here a short report on this interesting method [1]. According to it the dyna- 
mics of an infinitesimal test body placed exactly at the boundary of a Schwarzschild cell 
describes the dynamics of the whole universe. The physical considerations and mathe- 
matical treatment in the paper in question seem, however, rather complicated and 
will be simplified here without changing their proper meaning. 

From the differential equations of a geodesic, governing the motion of the test 
body, we need merely the fourth, and the expression for the interval, as the first inte- 
gral of the equations of motion 


dT dv dr dT 
ee ee er) S wl 
ds? dr ds ds ° ee) 


ele os 
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T being the time coordinate, s the proper time of the particle, and 
e = ] — 2mjr — ihre (5.3) 


In contradistinction to Lindquist and Wheeler, we take into account once again 
the term with the cosmological constant. 
Integrating Eq. (5.1), we have 


dljds = Ce”; (5.4) 


C, being an integration constant. Putting dT/ds from Eq. (5.4) into (5.2), we obtain 
with the help of Eq. (5.3) 


(dr[ds)2 = C2 —1 + 2 mr + 2472. (5.5) 


Hence 


(5.6) 


At first sight the similarity of Eqs. (5.6) and (4.2) is obvious if we assume 


EN (5.7) 
Putting C? from Eq. (5.5) into condition (5.7), we find 
E ae ney E 


After multiplying this relation by the mass of the test body, it acquires a simple 
physical meaning. On the left-hand side we have now the kinetic energy.of the test 
body, and on the right-hand side two terms which remind us of the Newtonian poten- 
tial energy and of a quasi;elastic potential energy determined by A and caused, as 
showed above, by continuous creation of matter. From the relativistic standpoint 
the mass m in (5.8) consists of two parts. The first, mp, corresponds to the rest mass 
and to other forms of matter of the central body of Schwarzschild’s cell, the other, 
My, to the Hubble mass-contained within the cell. The latter is of such a magnitude 


that the equation 
1 (ai\* _ ma 
De \ ds} 8 er 


is fulfilled (we neglect here the cosmological constant as physically meaningless). 


The condition (5.7) is thus satisfied as long as the mass mp in the universe is positive. 
If we had 


1g C <0, (5.9) 


we should have to do with a pseudospherical space which has been excluded from 
our considerations because of the above mentioned reasons. 
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We have rather departed from the paper of Lindquist and Wheeler. If we deter- 
mine the integration constant C4? from the condition that at the maximal expansion 
of the cell (r = ry) the test body comes for a moment to rest, we have (putting A = 0) 


GE = 2 mir, (5.10) 


and Eq. (5.6) becomes identical with Eq. (30) of paper [1]. It is also not difficult 
to prove now not only the similarity, but the identity of Eqs. (5.6) and (4.2). 

The differential equation describing the dependence of the time coordinate T 
on the radius of the cell is obtained from Eqs. (5.2) and (5.4): 


(dr/dT)? = e”” — e"/C,?. (5.11) 


Inserting into (5.11) Eq. (5.3) with A = 0 and (5.10) and integrating, we get just 
equation (29) of paper [1], namely 


Go — 2m) o -H7 E. 7 
T= Sa eae ed) arc cos Ro (5.12) 


fra Vr (o — 2m) + V2m (ro — 7) l 
Vr (o — 2m) =E 2m (o — 7) 
We have thus proved that there exists a full agreement of the results obtained by 


Lindquist and Wheeler and those deduced here. The present paper gives, however, 
a deeper insight into the dynamics of the universe. 


Conclusions 


Since the relations obtained in the present paper differ essentially from those 
of many papers on the cosmological problem where the field equations with the 
energy momentum tensor of continuous matter have been applied, we may conclude 
that the introduction of this tensor into the non-linear field equations of general 
relativity theory demands a careful examination in each particular case. 

From the finite radius of the universe, we may not conclude that it is the only 
existing one. We should suppose the existence of many closed universes so embedded 
into the” cosmical space of higher number of dimensions than four ’’ that their hyper- 
surfaces do not intersect each other. Since there exists no physical interaction under 
them, they are incapable of being observed, but this does not signify that they do 
not exist. 
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In this work a scintillation layer-detector with high detection efficiency for fast 
neutrons and low efficiency for y-rays is described. The detector is designed for the 
neutron measurements in the 5—20 MeV energy range. Its operation is based on the 
difference in range between protons and electrons of the same energy. 


Introduction 


In order to solve different problems in nuclear physics and elementary particle 
physics, high efficiency neutron detectors for measurements against a y-ray back- 
ground are desirable. At the end of 1958 some papers [1—3] appeared concerning 
the construction of neutron detectors in which discrimination of the Compton elec- 
trons from the recoil protons was based on the difference in the decay time of light 
pulses produced by these particles. 

In this paper another detector of the same type, a layer-detector, for effective 
counting of fast neutrons against a y-ray background is presented. After this detector 
has been constructed and examined, a brief communiqué on the use of similar de- 
tectors appeared [4], its data however, were not described in detail. 


Principles of work and construction 


The layer detector is based on the use of the difference between the proton and 
electron ranges of the same energy). The detector for counting neutrons of energy 
E, (E, < E, < E, max) consists of several layers of plastic scintillators A mm thick. 
Let us take E, as electron energy which produces scintillation pulses of the same 
amplitude as that from the recoil protons of energy E,. Thus a thickness A should 


1) The possibility of Go layer detectors for discriminating between short-range and long-range 
protons was previously examined in our laboratory. 


(675) 
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be chosen which is less than the range of the electrons of energy E,, but greater than 
that of the recoil protons of maximum energy Ey max 

Under these circumstances the electron of energy greater than Æ, may pass into 
an adjacent layer or traverse several layers. In the detector the adjacent layers are 
separated by intransparent foils, so that the light pulses from all the even layers are 
collected on the one photomultiplier and from all the odd ones on the other. If the elec- 
tron energy is sufficient to enable it to pass into the adjacent layer then in both photo- 
multipliers the pulses will be simultaneously produced. The electronic circuit discards 


Fig. 1. The geometry of the detector layers, light-guides, and photomultipliers. 


the coincidence pulses and allows the single pulses produced in any of the two photo- 
multipliers to reach the pulse-height analyser. This enables the separation of the 
recoil protons from the electrons produced by gamma rays to be made. The layer 
detector consists of 28 discs made from plastic scintillator (on polystyrene base with 
2% peterphenyl and 0,2% NPO added) 80 mm in diameter and 4 mm thick. The k 
value was determined from the above mentioned relationship for E, = 5MeV and 
E,, max = 20 MeV. The detector consists of two identical parts set one behind the 
other. In each part the light from the even discs is directed by means of light-pipes 
made from plexiglass to two photomultipliers FEU-29 and from the odd discs to 
two other FEU-29 photomultipliers (Fig. 1). In order not to let the light from the 
even discs fall on the photomultiplier which “sees” the odd discs (and vice versa) 
the side surfaces of the discs were divided into 4 equal parts and the two opposite 
surfaces were covered with aluminium foils. 

The detector has 8 light-tubes and 8 photomultipliers, each being put into an_ 
iron shield. There was no oil between the scintillators and the light-guides. This 
caused some slight decrease in amplitude at the output of the amplifier but 
simplified the work with the detector. 
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Detector electronics 


The H.V. supply applied to the photomultipliers was controlled separately for 
every of them and amounted approximately 1000 V. Allthe four photomultipliers 
which viewed the even or odd detectors had a common 2 kQ loading resistance. 
The schematic diagram is shown in Fig. 2. The pulses coming from each group of 


Fig. 2. Block diagram of the detector circuit: PU], PU2, — preamplifiers, U, U, U U, — amplifiers, 
S. S. — coincidence circuit, D. L. — delay line. 


photomultipliers, through the pre-amplifiers PU, and PU, and along the 50 meter- 
-long cables RK-50, were send into the amplifiers U, and U, (with a controlled gain 
of 1—1.5), which were designed to equalize the amplitudes among the groups of the 


Output 


6I71P 


6I71P 


———— 


= 0.97, amplification coefficient R = 0.4 for an amplitude 


Fig. 3. General diagram of {he summator: S 
A<0.6 V. 
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photomultipliers. The maximum amplitude of the pulses from Co-60 at the input 
of U, and U, was 0.1 V. From the amplifier outputs the pulses were led to the summa- 
tion device and into the amplifiers U}, U, (at the gain of 400). The coincidence 
circuit with the resolving time 0.4 usec gave the veto pulse which was formed by the 
amplifier with saturation and was put into the closing device. The veto pulse length 
was 0.6 usec. The summator was characterized by a summing coefficient S = 4142/41,2 
where A, z was the amplitude at the summator output, when the pulse A appeared 
at one of the inputs. 4,,, was the output pulse amplitude when the pulses 4/2 
appeared at the two inputs. 

In this work a summator was used (Fig. 3) whose coefficient S was 0.97 in the 
working range of the amplitudes. 

After the summator the pulses were additionally amplified and delayed, so that 
they came to the closing device at the center of the veto pulse. The single pulse pro- 
duced in any of the two photomultipliers passed through the closing device and 
came to the amplitude analyser. 


The basisc characteristics of the detectors 


When the coincidence circuit is turned off and the summing coefficient S = 1, 
the layer detection does not differ from that of a single crystal and has a high counting 
efficiency for neutrons as well as for y-rays. When the coincidence circuit is turned 
on the counting efficiency for neutrons is slightly lower, but the counting efficiency 
for y-photons is strongly decreased. In order to characterize the decrease in the count- 
ing efficiency of neutrons and y-rays, the K coefficient 


K = N-/Nt 


is introduced, where N- is the number of pulses at a certain level of the integral 
discriminator (coincidence circuit turned off), and N* is the number of pulses at the 
same discriminator level (coincidence circuit turned on). 

It should be noted that the discrimination coefficient K for y-photons when 
they fall along the detector axis is higher than that for the isotropic distribution. All 
the detector parameters were measured when the source was placed on the detector 
axis at a distance of 30 cm from the centre of the detector. This was due to the fact 
that the detector had to count neutrons and discriminate y-photons under oo 
geometrical conditions. 

The order of magnitude of K may be evaluated by simple reasoning. The estima- 
tion of the discriminating efficiency at E, < 10 MeV can be made by taking into 
consideration the Compton scattering process in the scintillation layer. (In a layer 
4 mm thick, the average electron-energy loss is 1.5 MeV). 

If the level of the closing device corresponds to 2 MeV, asis the case with the above 
detector, then all the Compton electrons of energy > 2MeV will be discarded. If the 
counting level corresponding to 2 MeV is fixed in the analyser, the higher energy 
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y-rays can be counted only by the double Compion scattering process in the discs 
of identical parity. This happens only if in both scattering processes the electrons 
do not cause the operation of the closing device and the sum of the energies of both 
Compton electrons is higher than 2 MeV. This results in a decrease of the discrimina- 
tion coefficient to some tenths. 

It should be noted that the exact determination of K is difficult because then 
the multiple scattering process must be taken into account, as well as the fluctuations 
in the ionization losses and the electron bremsstrahlung. It is clear, however, that 
the thinner is the layer, the greater is K for a given energy of y-photons and the thicker 
is the detector the greater is the probability of multiple scattering and correspondingly 
the smaller is K. 

The measurements of K for different counting levels were made for electrons 
and y-photons of energies 1.33 MeV (Co®) and 2.75 MeV (Na) respectively and 
4.5 MeV (p + B1!— C12)2). The results of the measurements are presented in Fig. 4, 


10? 


70 


15 Ep, MeV 
AE 5 10 Eo, MeV 


Fig. 4. Discrimination coefficients K, and K, plotted against threshold values of registration: py 
threshold of recoil proton registration, E, — threshold of Compton-electron registration. 
So se Sy a a a a 

2) This reaction as well as that of t+ d —> He + n were made by means ofa Cockroft-Walton-type 
accelerator of energy 300 keV. 
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where the K values versus counting levels are also shown. It is interesting to deter- 
mine the ceofficient K for electrons whose energies are sufficient to enable them to 
pass through all the discs in the detector, i.e. for energies > 20 MeV. 

From the measurements it resulted that the discrimination was higher than 1000. 
The discrimination coefficient is also given in the case when the veto pulse is caused 
by coincidence pulses from two photomultipliers looking at the same disc. This 
coefficient at the level of 1 MeV (Co®) was shown to be 100, but at the level of 2 MeV 
(Po-Be) it amounted to 1000. 

It follows from the preceding discussion that when the coincidence circuit was 
turned ona slight rejection of neutrons was observed in the layer-detector. This process 
was caused by a passage of a certain number of the recoil protons to the adjacent 
discs, as well as by double neutron scattering in discs of different parity. 

For the measurements of the discrimination coefficient K, for neutrons and for 
efficiency determination of the layer-detector, the integral spectrum of the recoil 
protons from 14.8 MeV neutrons produced in the reaction t + d—> He + n was 
measured, The recoil proton spectrum, when the coincidence circuit is turned on, 
is shown in Fig. 5, and the values of K, for different counting levels of the recoil 
protons is given in Fig. 4, The value of K, is seen not to be higher than 1.5. 
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Fig. 5. Integral spectrum of recoil protons from 14.8 MeV neutrons. Energies of recoil protons are 
indicated by arrows. 


The layer-detector efficiency when the coincidence circuit is turned on at the 
Jevel E, = 5 MeV was shown to be 15%. When the coincidence circuit was turned 
off, the efficiency at the same level was 20%, the difference being consistent with 
the estimated value of the probability of scattering in the detector of the same size. 
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The energy resolving power was determined on the basis of the recoil proton 
spectrum from the neutrons of energy 14.8 MeV (Fig. 5) as well as from the amplitude 
distribution of pulses of the relativistic particles which passed through the detector 
along its axis. 

The line half-width was 20—25%. The resolving power depended chiefly on 
the poor conditions of light collection on the thin discs, and on the difference in gain 
of the photomultipliers caused by not quite precise setting of their voltages. 

The inhomogeneity of light collection from each disc is characterized by 


A ae 2 (ines = Imin) 
7 Tmax ote Tenia 


where Tnax and Iain are the maximum and minimum current values of the photo- 
multiplier which were obtained with the irradiation of different parts of the detector 
surface. This value, when two photomultipliers looking the disc were employed with 
no reflector on the disc, was as high as 9%. 


e100 ee 


The measurements in which only one photomultiplier looked at the disc gave 
25%. When the disc was covered by a thin layer of Al by evaporation in vacuo this 
value fell to 19% and 85%, respectively. The above discrimination-coefficient de- 
scribed well the properties of the detector only when the summation coefficient was 
equal to 1. In fact, a decrease in S caused a decrease in the amplitudes of the summed 
pulses in the spectrum obtained with the coincidence circuit turned off, so that the 
shape of the spectrum close to its upper limit was caused only by single pulses. There- 
fore the discrimination coefficient in this case in the neighbourhood of the maximum 
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‘Fig. 6. K, dependence on the summation coefficient S. 


amplitude decreases to l and, vice versa, when S > 1, the discrimination coefficient K 
close to the maximum amplitude runs to infinity. The dependence of K on S is given 
in Fig 6, where the results of the measurements of K are shown (with Co®) near 
to the maximum amplitude for different values of S. In the electronic circuit of the 
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detector the summation device with S = 0,97 has been employed, so that the given 
values of K, and K, are somewhat lowered. The relationship between the discrimination 
coefficient and the gain U, and U, was also investigated; it was represented by a satur- 
ated monotonic curve. The working point was chosen at the beginning of the saturation. 
By comparing the layer detector of this work with that given in [4], it may be shown 
that the value K of the latter at the counting level of 2.2 MeV electrons (E, = 2.75 
MeV) is approximately greater by a factor 2. This is connected with the fact that the 
thickness of each layer in the detector in [4] is less by a factor 1.2 and the over-all 
thickness by a factor 2.7 than in the detector used in the present work. All this indi- 
cates that the detector [4] has a neutron efficiency twice as low as ours. 

What concerns efficiency of neutron detection against the y-ray background given 
in [2], there are only few data available for a detector of great thickness. K is probably 
greater than for the layer detector. The weak point of such a detector lies in a long 
dead-time (5 x 10-6 sec). This leads to a decrease in the neutron counting efficiency 
for high intensity y-photons and electrons. 

The layer detector can be employed for much higher intensities. The decrease 
in the efficiency of the neutron counting is caused by the length of the veto pulses 
and repetition frequency of pulses. In the detector presented above the veto pulse 
length was 0.6 usec, but it can be made as low as 10 8 sec, according to the length 
of the light pulse in the scintillator. 

The layer detector is most efficient in discarding the y-rays along its axis, however, 
it may succesfully be employed for decreasing the y-ray background distributed 
isotropically. For example by turning on the coincidence circuit, which was located 
in the building for measurments of the synchrocyclotron of the Joint Institute of 
Nuclear Research, the external background was lowered by a factor of 20 at the level 
corresponding to the 3.0 MeV electrons. 

A neutron detector of this type can be constructed as a ine proportional counter. 
In this case the amplitudes of protons and electrons of the same energy differ slightly. 
When the neutrons are counted at a certain level, K, in this case would be much higher ` 
than that for the same level when the plastic scintillator is employed. The proportional 
layer-counter could be particularly useful for low energy neutron counting. 

Thanks are due to N. W. Sizov for his help during the work with the Cockroft- 
Walton-type accelerator, as well as to D. K. Akimov.and V. A. Zapevailo for their 
advice on the construction of the electronic part. 
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ON THE STATE OF A FERMI-SYSTEM WITH CORRELATION 
OF PAIRS OF PARTICLES WITH PARALLEL SPINS 
Il. THERMODYNAMICS 


By ZYGMUNT GALASIEWICz* 
Joint Institute for Nuclear Research at Dubna (U.S.S.R.), Laboratory of Theoretical Physics 
(Received April 6, 1960) 


The thermodynamics of the “anomalous” state of a Fermi-system (in which the cor- 
relation of pairs of particles with parallel spins is taken into account) has been investigated. 
The transition temperature T, to the “anomalous” state has been found, as well as the 
temperature dependence of the specific heat for T ~ 0 and T ~ T,. For T= T, the 
specific heat has a jump. In addition formulae for the paramagnetic susceptibilities for 
T ~ 0 and T = T, have been obtained. 


Introduction 


In paper [1] (being part I of the present paper) the possibility of an “anomalous” 
(non-superconducting) state of a Fermi-system is studied. This state is connected 
with the creation of pairs with parallel spins. Creation occurs if the electron-electron 
interaction is attractive; we find that in the expansion of the interaction term in spheri- 
cal harmonics only the coefficients with odd indices give a contribution. 

If in the expansion of the interaction in spherical harmonics only terms with 
even indices J,,, are retained (for instance Jp in [2]), then the compensation equation 
for the pairs with antiparallel spins has two solutions: a trivial one and a nontrivial 
one (giving the superconducting state). However, the compensation equation for pairs 
with parallel spins has only the trivial solution. 

If in the expansion of the interaction in spherical harmonics only terms with 
odd indices J,,,,, are retained (paper [1] — I), then the compensation equation 
for the pairs with antiparallel spins has only the trivial solution. On the other hand, 
the compensation equation for pairs with parallel spins has two solutions: a trivial 
one and a nontrivial one (giving the “anomalous” state). 

The solutions of the compensation equation and of the secular equation for the 
collective oscillations are in the case J,,, even functions of p[3,4] (the superconducting 
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state, longitudinal oscillations) in the case J,,,, odd functions of p[1,5] (the “anoma- 
lous” state the oscillations of the type of the transversal oscillations). 

In papers [4], [5] it is proved that for the general form of the interaction term 
the collective oscillations split into two branches: for pairs of particles with antiparallel 
spins and for pairs with parallel spins. The terms of type J, give a contribution to 
oscillations of the pairs with antiparallel spins only (the solution of the secular equa- 
tion the even function of p therefore the longitudinal oscillations), the terms of type 
Ton4.1 give a contribution to the oscillations of the pairs with parallel spins (the solu- 
tion of the secular equations the odd function of p therefore the oscillations of the 
type of the transversal oscillations). In paper [1] it is shown that if 3 7, > Jy the energy 
of the “anomalous” state is lower than the energy of the superconducting state and 
certain characteristics of this state are studied (for T = 0). 

Here we want to examine some characteristics of this state for temperatures 


TH 0. 
1. The compensation equation 


The compensation equation has the form (see [1]) 


2 (p) (5, o) — EO cG, o) =0, () 


where 
< 1 p pad ra Il >, 
C(b, o) = — y $ M(B, — p: — Pp) DG, o) (2) 
p 


and 
® (p, o) = u (p) v (p, 0) (1 — 2 n3,), 
F(p) = v? (p, o) (1 — ng.) + u? (p) nze» (3) 


l AiR I) iA a Med > M PERSI sy 
O(p) = e(p) + 7 ‘`F [2I (p, p's B's P) — 1 (b: P's P B')) F(p’) (3a) 
E 


. . Es . . 
since the expectation value of aż, az, is not zero but nz, [4]. The functions {u, v} 


are the coefficients of the Bogolubov-transformation for transition from the Fermi- 
op + . are 
operators azs, a3, to a3,, a,; these fu, v} satisfy the condition 


u? (p) +v? (p, 0) = 1. (4) 
Similarly as in paper [1], we put 
I(p,—p3— psp’) = pa L, (Ps pP’) P, (cos y) ~ L (p, p’) cos y 
= Í (pp) (cos O cos O’ + cos g sin O sin O’), i 
C (p, o) =0 c08. 8% (p) -v | (5) 
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where O, O' are the angles between the vectors p, p’ and the polar axis, and ọ the 
angle about this axis. 


Eq. (2) after taking into account (1), (3), (4) can be written in integral form, 
as follows 


PptaA x 
deg LPR cos? O’ sin @’ d@’ dp’ i } 
Bi rei) J an STO I E ENEA O 
where 


l S 
ne = n = ————, Q(p) = Vy cos®O F t): (2) 


eMPT 4. 1 


Hence we get finally for y as function of temperature the equation 


f aes p Sere 
Jae x? dx 
a) te) ae Ë 


(8) 
af St th = == 
x? dx oO 
V2? —y rat 
where 
Peon dn pr 
Oy, = TE > O dE n2 E' (9) 


this means that we assume the interaction to be attractive. 
3. The dependence of y on the temperature and the transition. temperature to the 


“anomalous” state 


After a change of variables Q = y xch p, we get similarly as in [6] 


agge ain 
y (T) a ae 1 
1 00 
is sfe dx 3 ("+K A j (10) 
as 5 mil 
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Hence for T ~ 0, making use of the asymptotic formulae for Bessel-functions 


with great arguments, we obtain 


v (0) 18% T p e-WOslT dy —9 / Z y (0) ( ) 
yD) =V50 J” ora =oy/ 30 (VF) (5m 


T (o) Ga 
p (T) = p (0) e XTIvO” (lla) 


where ® (z) is the error function. 


Now we want to find the critical temperature T, (for T > T, y(T) must vanish): 
Moreover, we want to get y (T) for T ~ T,. For this purpose we make use of a suitable 
representation of the Bessel function Kg (z) [7], [6]. Thus we obtain 


1 co 
n 2O als nT we dadon 1 ae 
vl) =) me ak | >, s 


ps 


Poe |p 02 


where In y = 0.577, .y = 1.8. 
Putting in (12) y = 0, we find 


y (0) = 2o e'h(e-Hlei)3 = eth == (13) 
Hence 
wd = 2.45 7 (13a) 
Komal ues F- 
T _12@)(v()\’, 526) (v(n\* 
naz Ol O) 4 £8 Te, he aa 


where ¢ (n) is the zeta function. 
Therefore for T ~ Th 


y (T) Tc — T 
LAGU 
if i Te 


From (13a) and (15) it results that the function y (T)/T, decreases from 2.45 
to 0 in the temperature interval (0, T). 


(15) 
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4. Entropy and specific heat 


Now we want to obtain the formulae for the temperature dependence of the 
specific heat, especially for T ~ T,. Therefore we must first get the formula for the 
entropy. 

We start from the expression 


S=—2 » [mg In ng + (1 — ng) In (1 — ng)). (16) 
P 


Inserting (7) into (16) and passing to integral form, we obtain 


1 (eo 
Si dn 1 ðn 
0 yx 


After changing the variables similarly as for (10), we have 


1 oo 
dn y 
s= 4-7F T x? dx D (—1)”+1 K; bes } (18) 


mji 


Hence for T ~ 0, making use of the asymptotic formulae for Bessel functions 
with great arguments, we get 


A 2 
mmg Poast a] o 


Now we have to find s (T) for T ~ T,. Since for T> T, p(T) > 0, we write (17) 
in the form 


sl(0 s1 (0 sll (0 sv) , 
s(y) = O, i sy a yet O, 
= 2/3 n? < r|1— ey + z zo (2); (20) 


Having found s (T), we can obtain the temperature dependence of the specific 
heat C(T) for T ~0 and T ~ T,. We denote the specific heat in the normal state 
fore! = 1", by | i 

o an 
Cn ()=3 = na TE Te 


For temperatures near T = 0, we have 


Idon ol 5 (5) a 
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however, for temperatures T ~T, (T < T) 


CAB) _ T — T, 
CAT) ae Te 
From (22) we see that for T= T, the specific heat has a jump, since 
C, (T)/C, (T,) = 1.83. 


(22) 


5. The total momentum of elementary excitations 


Let us consider the total momentum P of elementary excitations in the case 
of macroscopic motion with constant velocity ù 


P= Dk (niy + mag) (23) 
k 


Now n = n (Q — ki). For small velocities 


PREG, 2 Bie) ae (24) 


In integral form, we have finally 


1 ue 

~~ 2 dn a oe xa > 2 dn 
— pF Z fa x ) + uy 2x7] dx 5Q(x, 6) dl, (25) 
ò —0o0 


where ù pis the component of vector u perpendicular to the quantization axis of elec- 
tron spin and ù — the parallel component ¢ of u. We see from (25) that the vectors 
p and u are not parallel and we cannot write p = N, m ù and define the number N, of 
normal electrons in the “anomalous” state. 

Consider p for T ~ 0 and for T ~ T,. After changing the variables 


¢ = yx sh p 
we get 


dn 2y ch gy dp 
= PF OR T a (a0 a a ad a (ore e eiry CO 
Hence for T ~ 0 


P=Py +P Py (27) 
where 


Ad 3 gean ne. iges 
i= gob Seve (| a= NS ye (| mü, (28) 


Pery S 
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m is the mass of the electron, N the number of electrons in unit volume, 


1 2 2 dn oe 87 pr 
m 3°" GE 30r) (29) 


since in formula (9) E’ ~ pz/m. 


From (28) we see that the perpendicular component of p dominates for T ~ 0; 
it is three orders of magnitude greater than the parallel component. 
In order to obtain (26) for T ~ T,, we may use the following indentities 


l 
3 eich dopisane, sie T T 3T 
Ak e| ri (evxchgiT +. 1) (e—vacholT +. 1) 6 T v\ (29) 
OT \T 


and 
Of ( 
Se oT 
+ fu mii (evechelT + 1) (edt + ]) 52 [¥)\ , (30) 
vaca 
where f(T) 
F(T) = r faf TI TO (31) 
For T ~T, 
ex T 7 4 9258.95 e 
f(T) = eae eo (32) 
Hence A for T ~ T, 
Y 3 
a +5752) (z-z s= aul) m 63) 


Finally p for T ~ A (T< T) 


» 1 = > 1- ree) T k 
P = Nm E + T (llu — 8u) T, | 5 (34) 


We see that even for T ~ T, the parallel component of p is smaller than the per- 
pendicular one. 


6. Paramagnetic susceptibility 


From paper [l], we have for the paramagnetic susceptibility 


ae y (T) 
D= E E ma = i (35) 
where we get y (T) from (10). 
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For temperatures T ~ 0, we obtain from (lla) 


T \3 
x ( gyre =) ae [2 — 1.96 e- Se: e7? (5m) lk (36) 
For temperatures T ~ T,; we get from (15) 
ws LON a —6/es aS) 
XL) ae TE £ 3.25e T, : (37) 


With the elementary excitations considered here (type J,,,4) we could explain 


the Reif experiment [8] which gives the T-dependence of the paramagnetic suscep- 
tibility, if we put @,; ~ 3.5. 


It is a pleasant duty to thank Prof. N. N. Bogolubov for proposing this problem 


and helpful advice and D. V. Shirkoy for valuable discussions. 
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SCATTERING OF SLOW NEUTRONS IN FERRIMAGNETIC 
CRYSTALS 


By Jerzy KociNsk1 
Ferromagnetic Laboratory, Institute of Physics of the Polish Academy of Science, Poznan 
(Received April 25, 1960) 


Cross sections for scattering of slow neutrons by spin waves in ferrimagnetic crys- 
tals have been calculated. The type of ferrimagnetics with spins pointing in the same 
direction while at absolute zero temperature has been treated. The formula for the scat- 
tering cross section given by Halpern and Johnson for the ferromagnetic case has been 
generalized to lattices with two kinds of spins. By means of it the cross sections for various 
scattering processes have been calculated in the manner used by Maleev in the ferro- 
magnetic case. 


1. Introduction 


The scattering of slow neutrons by spin waves in ferrimagnetics has as yet not 
been treated theoretically, though the subject has been investigated in the ferromagne- 
tic case by many authors, at last by Maleev (1957). 

We shall investigate this scattering for ferrimagnetics with two kinds of spins, 
pointing in the same direction while at absolute zero temperature. We begin with 
defining the states of the ferrimagnetic domain in the manner introduced by Dyson 
(1956 a) for ferromagnetics. We express then the spin operators of the scatterer in 
terms of harmonic oscillator operators. This permits to extend the general cross 
section formula, given by Halpern and Johnson (1939) for ferromagnetics, to the 
ferrimagnetic case. We shall derive explicit cross section formulae for various scatter- 
ing processes and discuss the cases of elastic scattering and single magnon scatter- 
ing. 


2. Spin waves in ferrimagnetics 


We will investigate neutron scattering in ferrimagnetic crystals of simple cubic 
or body centered cubic lattice, with two kinds of spins s, S. These spins are attached 
to the lattice points in such a way that the nearest neighbours of epingstdte S spins 
and vice versa. In our previous paper (1960), physical spin-wave states in ferrimagnet- 
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ics of the form 


|A> = M (4,9) (a 3f + C, SHOS (1) 
A 


have been introduced. To investigate neutron scattering, we have to express spin-wave 
states and spin operators in harmonic oscillator terms. To this end, in analogy to 
Dyson’s treatment of the ferromagnetic case (1956 a), we at first define the set of 
states 


ju = LET (2 8) (uy G9 TT 2S) HH (ug!) (S$ 10> 2) 


J 
for a ferrimagnetic consisting of two magnetic sublattices. Every sublattice consists 


of N atoms, the atoms j of one sublattice have spins s, the atoms k of the other sub- 


lattice have spins S. The integers u; attached one to each sublattice point j can take 


the values 0, 1,... 2s, while u, the values 0, 1,... 2S. These states are arthogonal and 
the formulae 


<v u> =F, 6,, Fu =F (u,) I F (u), (3) 
j k 
where 
F(uj)=1- e re 2 = F(u) =1- (1-3) r fr ae 
(4) 


determine the normalizing factors. The states |A> and |u> are orthogonal unless 
2y + Dee = DAL (5) 
j k A 
This means that one unit of reversed spin (of one or other sublattice) corresponds 
to one spin wave, as in the ferromagnetic case. 
Now we pass over to the ideal spin-wave model. In this model to the sites j and 
4 of the sublattices are attached harmonic oscillators with states labelled by integers 


he taking values from 0 to oo. The oscillators have creation operators qt, ĝt and 


aA operators Nj» ô, satisfying the relations 
l,» 7,1 = ni > h] = 0, lj,» na] = 6; ja? (6) 
ln 3.) = [ÔE OET = 0, > lro DL] = 4,5 
ng ty = U; » oF oo (7) 


The integers u; , u, are connected with the u fai of the physical model (2) by the rela- 
tions 


4 = 2+ 2M | (8) 
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The meaning of (8) is that the sum of the quantum numbers of the oscillators of either 
of the two ideal sublattices is equal to the sum of units of reversed spins in both sublat- 
tices of the model (2). We shall need this type of correspondence between the ideal 
and physical models in defining the ideal spin-wave states. A complete set of states 
for the whole ideal system is given by 


ju) = [E (DE GN T D BEN] | 0). (9) 


These states are orthogonal and normalized. 

For both sublattices there is one reciprocal lattice consisting of N vectors 1 which 
satisfy (i) the periodicity condition that 7 is identified with 2 + n u for nè = N, and 
any two rh ù of the lattice, and (ii) the condition that (n/2z) i .jand (n/22) À - kare 
integers for all the sublattice vectors j, k. 


In the ideal model we define in the reciprocal lattice the complete, orthogonal 
and normalised set of states 


|4) = IT [D Gf) D5? (ÊF) | 0) (10) 
= H ((a,)~* BE) | 0) TT 10,4 (7) | 0) = | a) | b), (10a) 

where &* and B} are defined by 
ât =N Dep G4.) Hy it = N-e (iiA pr (1) 


The a}, Bt are creation operators for harmonic oscillators, whose states are labelled 
by the integers a} b}; they satisfy the relations 


[a,, at] Sa yne [By Br] aa Days (12) 

Ot a = Ay, By Piz b,. ~ (13) 

We establish a one-to-one corespondence among the physical states |4 > and the 
ideal states |A) in such a way, that to the physical state (c, S+ CSF) |0> there 
corresponds the ideal state at ĝt |0). That means that a pair of harmonic oscillators 


attached to a reciprocal lattice point represents the ideal spin wave. The integers 
ap b, are equal and satisfy the relation 


a, = b1 = Ar (14) 


where A, is the number of spin waves with wave vector À. 
In the linear approximation the Hamiltonian of the system in terms of the ideal 


spin-wave model has the form 


=F +bE Lt ep te + Bt BD (15) 
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with spin-wave energy 
&=3Js {(s +S) [n (Yo = Ya) yol E (G= s)? [7 (Yo Eg Ya) =~ li? +455 yi3 (16) 


which has been calculated in the previous paper (1960). For small values of A we have 
the expression 


a, — Lalas? £5% Fob? + 2557084) po 


6645) (16a) 


In the following we shall use the spin-wave energy expression with only nearest neigh- 
bour interaction and drop the term with 7. Then we have 


era Ja yoô? A?, (17) 


where J, is the exchange integral between nearest neighbours, yọ — the number of 
nearest neighbours of an atom, ô — the distance to the nearest neighbour. 

We shall need the volumes of the unit cells of the sublattices. For the simple 
cubic lattice the basic vectors of sublattice are the face diagonals of the cubes and 
the unit cell volume is equal to 2 a?; for the body centered cubic lattice the sublattices 
are simple cubic and their unit cell volume is aè. 


3. Spin operators expressed in terms af harmonic oscillator operators 


Dyson (1956 b) has given a general method of transferring any operator O of the 
physical model to the ideal model. Following Maleev (1957) we write the transferred 
spin operators in the form- 


f= —s Yes (if 7) [3y9 — (SN) Dts Sl (18a) 
$+ = (2s N-D} z exp (~i. j&i, — (18b) 
p eeN Dorp (— ii J) [@_, — (2s N)? D hrotu abe], (180) 
$=- Sgap (—i 7-8 [89 — (SN): Dhabi, 09) 
S$ SINN Deap if. i) Br, (190) 
S= 2SN} z exp (— iÀ. K (B_,—(2SN)7 Z Phroru Bx Bal (190) 


for both sublattices. The z — axis is determined by the direction of magnetization 
_ of the ferrimagnetic domain. The above expressions are used in calculating the cross 
sections. 
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4. Scattering cross section of ferrimagnetic spin lattice 


Halpern and Johnson (1939) have found the cross section formula for magnetic 
scattering of neutrons in ferromagnetics. We shall generalize this formula for the 
ferrimagnetic case. Thus the Halpern and Johnson formula (8.22) for the wave func- 
tion y, takes for a ferrimagnetic the form 


M\t/ M GER | ee og 
ee (i) a 2, p MARA, (20) 
n' A’ 


where M is the mass, k and k’ represent the values of the wave vector before and after 
the scattering, R the position of the neutron; n’ is neutron’s spin function and (A’) 
the domain’s state after the scattering. The magnetic interaction potential has the form 


A 


n Arh? 1 fei: Be s 
prea M | Fi (q) (n'a'| Den 8; (Sn — € . Sne) |na) + 


+F (Q t'l Le ari O md) | (21) 


where rọ = e?/mc? is the classical electron radius, $, — the neuti on’s spin operator and 
y — the value of its magnetic moment in nuclear magnetons; q = = k — F is the scat- 


tering vector and E gla; the states |a), |b) are given by (10a), andF? (q), F&(q) 
are the formfactors for the two-sublattices, given by 


FE (a) -f dry} 3 tas ai aal w F$ (Q= -f ac yi È A i z Yrs 
(22) 
the sums z being over all the electrons of an atom j or k with total spin S or 
S; ey. and p, are correspondingly the spin operator and _Position vector of the 
electron v, and p,,, are the wave functions of atoms i, k. For the operators $, 
5. the expressions (18) and (19) are used. 
By means of (20), when a neutron is scattered from a state with wave vector E to 


one with k’ the differential cross section per atom has the form 


do Lik: M \2 Kn A’ ie 
dO ~ 2N E ea) Dara cat ey 
nA’ 


See exp (—2W/) is s the Debye-Waller factor. The flux of incoming neutrons has been 
assumed normulized to unity. In the following we assume unpolarized neutron beams 
so that calculating the scattering cross section (23) we have to average over the direc- 
- tions of neutron spin. 
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5. Cross section formulae for various scattering processes 


By virtue of (21) the cross section formulae for ferrimagnetic crystal differ only 
by numerical factors from those for the ferromagnetic case. Thus in the subsequent 
part of the investigation we may closely follow the treatment of Maleev (1957). 

The form of the operators (18) and (19) implies the following scattering processes: 
elastic scattering and four types of inelastic scattering; with creation or annihilation 
of a single magnon, with creation of one magnon and annihilation of another magnon, 
and with creation of one magnon and annihilation of two other magnons. By means 
of (23) and (21) we obtain the corresponding cross-section expressions for monocrys- 
tals: 

Elastic scattering: 
doy l 


So = ay Be E) 


pa eins 
x & (9) (2 a wh <ar] + SF$ @(2 ay <ar) | (24) 


Here and in the subsequent expressions we have replaced 6, by a, by virtue of (14). 
Inelastic scatterings: 

do 1 gels 

SEIE y e-i- 

7 


dO QNP k 
x (ĊF? (q) + SPF§ (9)? (Cau>r +1) (25) 


dod á Les eS 
J 


x <|> F* (q) až Koai 5 + St Ft (q) af lez Les a) D 


k ary 2e—2W (1 + e$) x 


x 


A+ pb AF u 
29) 
do41,— leeks „> > >|? 
a l = NS T roy?e-2W (1 == e?) > eS ee x 
J 
x (FE (a) + F$ (9))? <an (a, + 1)>r (27) 
dorr aral 


k 
dO = CN) a roy2e—2 (1 + ez) x 


2 
Y eilg—at oto) +j 
7 


x (SEFË (g) + SF} (q)? <la, + 1) aye) 7. (28) 
The symbol <...> 7 means an average taken at temperature T so that <a u>r is 
equal to Planck’s function n(u) = [exp (e/kT — 1)]-2. 


For the present we shall discuss the elastic scattering (24) and the single magnon 
scattering processes (25) and (26). 
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6. Elastic scattering 


In the formula (24) for the elastic scattering cross section, the expression 
z <a;>r/N (s + S) represents the relative deviation of spontaneous magnetization 


at temperature T from the magnetization value at absolute zero temperature. For low 
temperatures it is given by 


ae a vee ee! a H SAT L Cd. (29) 
N(s+S) s+S\2xT, 2]) -s+S ) 
where C(%) is the Riemann zeta function and 
_ TssSyov 
oe 3(s + S)k re 


is defined in analogy to the approximate Curie temperature for the ferromagnetic 
case (Dyson 1956a); v = ô? vp ê and v is the volume of the unit cell of a sublattice. 
In the following we shall neglect in (24) the square of (29). Since for large N 


(ef 4]? = (2m)? Nog? 3 ô G — 2), (31) 
j T 
where 7 means 2% times the reciprocal lattice vector (Weinstock 1944), expression 
(24) takes the form 


doa _ L gyre Q — e F(A [L — 2 CCD) 
x (27) v9! K 6(q — 7) (32) 
with 
F; (4) = (s Ft (q) + S F$ (Q)? 
Fy (9) = Is F, (q) + SFs (9) + (s + 5) F? (0) F3 @IA (0) - (33) 
To obtain the cross section formula for polycrystal, we integrate (31) over the 
directions of 7 and get 


do, 1 2 = Er 
70 mao Olin. F, (0) D — 2A) G(T) x 
x (1 — e?) e724 (2 sn = — z) : ie) 


F . : . 2 
. where g, is the number of various reciprocal lattice vectors with length t. For ez 
the geometrical conditions of scattering yield the expression 


0 2 
s= [s gme TEE (35) 
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where @ is the scattering angle satisfying the condition 2k sn O/2 = q, and ¢ is the 
angle between k and the z — axis determined by the direction of domain magneti- 
zation; ® is the polar angle of i , measured in the plane perpendicular to i, from the 


direction parallel to the vector É x 2. 
Integrating (33) over ©, ®, we obtain the total polycrystal cross section 


1 = T E 
o= apa 0r Ca) a Y E FO) [1 2A (2) G(T)] e2 x 


TS 2k 


xı + cos? ¢ + — i ay — 3 cos? J ; (36) 


7. Single magnon inelastic scattering 


Using (30) and neglecting 2’ <a,>,/Nin comparison with 1 in (26), the expres- 
sions (25) and (26) may be written in the form 


fe 


Beds IN JA roy? Cm aw atl + ez) F(q) [ro az D ar a Y 
x (27)? vo 6G F ji +7) (37) 
with 
FQ) = (GH @ +S FS @)?. (38) 


Since the spin wave theory is valid for u ô <1, we may replace g by 7 in all the terms 
of (36) except in the 6-function. The 6-function expresses the condition of momentum 
conservation 


kak rs iO), (39) 
The vectors $, k’ and jt satisfy also the condition of energy conservation 


; 2M 2M __sS 
k'a — BB = So ey FF 3G 3 FS) 12708 mai. (40) 


where the constant a> 1. 

There is no essential difference between the expression (37) for single magnon 
scattering, and Maleev’s formula (12) for the same scattering in ferromagnetics. 
Instead of the numerical factor S F? (q) in formula (12) of Maleev, we have in (37) 
the numerical factor (1/2) F (q), with F (q) defined by (38). To make Maleev’s formu- 
lae suitable for our case we have to replace the factor SF? (q) by (1/2) F (q), and the 
ferromagnetic expressions for Curie temperature T, and spin wave energy £, by cor- 
responding ferr‘magnetic expressions (30) and (17). “By adapting the Guttering formu- 
lae to polycrystals, we must yet remember that the function G (T) is defined in the 


ati 
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ferrimagnetic case by (29). Since Maleev’s discussion of scattering holds exactly for 
the ferrimagnetic case, we shall not repeat it here. 

I should like to thank Professor Szczepan Szczeniowski for critical discussions 
on the subject of the paper and for reading the manuscript. 
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SPIN-SPIN RELAXATION TIME IN FERROMAGNETICS 
By Janusz MorkowskI 
Department of Ferromagnetics, Institute of Physics of the Polish Academy of Sciences, Poznan 
(Received April 25, 1960) 


Spin-spin relaxation time caused by pseudo-dipolar coupling between spins is 
calculated. It appears that for spin waves excited in ferromagnetic resonance in metals, 
having wave vector magnitudes of the order of inverse skin-depth, the relaxation time 
is determined practically by processes of splitting one spin wave into two (and confluence 
of two into one) only, leading to relaxation time of the order of magnitude 10-8 sec for 
nickel (at room temperature and magnetic field strength of about 5000 Oe). 


§ 1. Introduction 


In the first part of this paper (Morkowski 1960, hereafter quoted as I) a formula 
for the spin-spin relaxation time has been derived and some of its consequences qis- 
cussed. In the present paper we calculate the relaxation time for spin waves with 
non-zero but small wave vectors. 

As a starting point for our calculations we take the formula (5.20) of I: 


AN H (1.1) 
where | 
7 = ioe {2A}, Am (Gp — Ge) Ô (En + E» — €o) òl +d—-B + 
4 Ad Az (ae + a2 +1) 8(, — & — &) 8(4 =¥ — @)} (1.2) 
and 


I 24r p—y x—y iy, ae Vat x 
Oh Hee + Yy ) ( ux Er ah) 
lid HUT 
x [a2(@o + a° +1).— aay] 6 (en +x — & = &) x 
; x ôl +% —¥—-7) (1.3) 
(701) 
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with the same meaning of the symbols as in I (compare equations (3.23), (3.24), 
(5.6), (5.11), (5.12), (3.29), (3.28) of I). 

Our intention is to calculate the values of the spin-spin relaxation time which 
are characteristic for the conditions existing at microwave resonance in ferromagnetics. 
In part I we have calculated the relaxation time for spin waves with zero wave vector, 
which are the only ones excited by the microwave magnetic field in perfectly homo- 
geneous and non-conducting ferromagnetic materials. It is well known, however, 
that in usual conditions spin waves with non-zero wave-vectors may be excited. Kittel 
and Herring (1950) have shown that in conducting ferromagnetics non-homogeneneous 
modes are excited by microwave fields. In spin-wave theory language it means that 
are excited spin waves with non-zero wave vectors u, whose magnitude equals approxi- 
mately u œ 1/6, where ô is the skin-depth of the microwave field (Kittel and Herring 
1950). In presence of other inhomogeneities (numerous crystal imperfections, influ- 
ence of boundaries) nonhomogeneous modes corresponding to spin waves with non- 
-zero wave vectors may be excited also in non-conducting ferromagnetic media (com- 
pare for instance Clogston et al. 1956). Here we are interested mainly in ferromagnetic 
metals; thus we shall calculate the spin-spin relaxation time for spin waves with u œ 1/6 
(usually ô is great compared to the lattice constant a thus ua œ a/6é is small compared 
to unity). 

In evaluating the relaxation time we shall adopt the usual (compare for instance 
Van Kranendonk and Van Vleck 1958) approximations, mentioned also in I. We outline 
them shortly here. Both exchange and pseudo-dipolar interactions will be confined 
to nearest neighbours only (with coupling constants J and D respectively). Only 
long spin waves — compared to the lattice constant a — will be considered, 


ie. |A| a < 1. Sums over reciprocal lattice vectors A will be evaluated by transform- 


ing them into integrals. Approximations allowed by the condition ua œ ajô <1 
will be used. 


§ 2. Relaxation time as determined by 3-wave processes 


The contribution (1.2) to the inverse relaxation time, determined by 3-wave 
processes of splitting one spin wave into two others (and reverse processes of conflu- 
ence of two spin waves into one) will be dicsussed here in detail. 

The coefficients A+ in the nearest neighbour approximation and for long spin 
waves are (c.f. I): 


K 
A = q Ja? E [*(%x + iy) + Az (îs + iA,)], (2.1) 


_where q equals 0, D/J,% D/J for simple cubic, body centered cubic, and face centered 
cubic lattices, respectively. Summations over one index, say 6, may be performed 
in (1.2) explicitly on account of the existence of the ô (ù + p — o)-factors. Sum- 
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mation over v is replaced by integration. In order to perform this integration we intro- 
duce polar coordinates in v space with polar axis parallel to the vector 22). Introducing 
the quantity 


1 l il 
Sous aed oD dQ, , (2.2) 
B 


we get after somewhat laborious calculations 


ed 
1. 1 l 
PR ea reo * 


J n/2u ~ 
E dya 
3 
x (3 n? + ut — pv? + v’) ody + Ka vir, (2.3) 
E Lyp 
where 
eee D 
= To nh a ie SE (2.4) 


n = L/w and s equals 0, 2/9, 1/16 respectively for simple, body centered and face 
centered cubic lattices (L and w have the same meaning as in I, i.e. L =.2 upH, œ = 
= J S a? with ug — Bohr magneton, H — magnetic field strength in oersteds, J/2 — 
“exchange integral” and p = (kT)-1, e = BL). 

The second integral in (2.3) is not given explicitly as it is different from zero only 
if condition y? > 2 7 is fulfilled. In common experimental situation and for u œ 1/6 
this inequality cannot be satisfied, thus the integral in question vanishes”). In order 
to simplify the calculations leading to (2.3) we have replaced for integration the first 
Brillouin zone by a sphere with radius R, having the same volume as the real one. 
The result is thus not accurate for values of u close to n/2R. This is, however, of no 
importance as in practice values of u are greater than 7/2R and 1/1? changes rapidly 


with u. 


1) This transformation to polar coordinates may be realized for instance by the equations: 


Se ey ae e oy + 2 vys 


Hilo Ho 


ù 
PSL E o Vy = 


tga vpt ngs m= Ut u); ry = v sin B cos p, 
lai H 
vy = v sin ® sin p, vy = v cos Ò 


2) Condition u? > 2n, would be realized only for magnetic fields smaller than about 10 Oe with 
pio. 
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For u < 7/2R Eq. (2.3) evidently predicts VT = 0, thus, in particular, spin 
waves with zero wave vector cannot be relaxed by 3-wave processes, as it was mentio- 
ned in I. Merely in order to obtain this result we have paid some attention to the inte- 
gration limits. For significant values of u (say n aju smaller than about 2) and for 
not very high temiperatures we can safely pass with R to infinity and it will be done 
below. 

Confining ourselves to the case u? < 2 , we introduce a new integration variable 
E = Bw v? in the first integral of (2.3). Integration will be performed by well known 
method of expansion of Bose factors into geometrical series. 


The result is 


: =Ap* {Gey | czw" sT = liek n| x 
x [Z (u) — Z,(uy)] + (Bo)? [2 (1/24)? — u°] x (2.5) 


x [Z3 (u) — Za (u1)] + 2 (Bm) [Z3 (u) — Zs (ull 


with the following abbreviations 
u = Bo |(n/2u4)” + n] 
u = u + Bo (n + p?) (2.6) 


and the definition (compare Dyson 1956, eq. (66)) 


Ad BG Dab ef (2.7) 
p= 


The strong dependence of ine on u (in contradistinction to the analogous quan- 
tity for 4-wave processes discussed in next section) is worth to be noted here. 

Expression (2.5) may be simplified for the important case u <V27; we get 
then approximately 


5 = A (Bo)? (ale) (22a (u) + 2u Z (u) + u? Zo (u) 28) 


lad 


with u œ fpo (n/2u)?. For values of u small compared to unity, the leading term 
of the right-hand side of (2.8) can be found by means of asymptotic expansions for 
Z,, (x) given in the Appendix. The order of magnitude of Ho is thus 


5 ~ 2 (2) A (bo)? nju 


_¢2) J [D\? [kT\? na 
= (a) y 2) 
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In order to obtain deeper insight into our results, we insert into (2.9) the values of the 
parameters for nickel (S = 1/2, J/2 = 230 k . 1°K and, up to the order of magnitude 


me (g — 2)?, g = 2.2; compare reference I) getting approximately 


ne 108 a m (2.10) 
u 
na? ~ 0.58 x 10-° H/1 Oe. 
It must be noted here that to 1/7 may contribute also some terms of the form 
(2 SN%)~* 2 Ty at, 44%, %, Xia Z from the Hamiltonian I, (3.22), but they 
were eae here as they contain the small factor NT! X a9. 
A 


§ 3. Relaxation time as determined by 4-wave processes 


In this section we shall discuss a contribution 1$? to the inverse relaxation time, 


determined by mutal scattering of two spin waves. Results which will be obtained 
here extend for the case u œ 1/ôthe ones given previously for spin waves with zero 
wave vector (I, Section 6). 

In the nearest neighbour approximation and for long spin waves coefficients 


involved in yey are given by 
a Ja (5 > 2 heel Poe 
— 6t- Ñ + 60: A) — 3p (622 + 302 + 303 — 420s — Ot2Az + O0242)} (3-1) 


with p equal 4 D/J, 0, } D/J respectively for simple cubic, body centered cubic, 
and face centered cubic lattices. After performing summations with respect to one 
index, say % simplifying some expressions with use of properties of Dirac 6 functions, 
passing from summations to integrations (extended over all the reciprocal lattice), 
and ordering according to powers of the parameter p, we get finally 


1 1 1l 1 A 
— = | — ar Gams oe 
ry Ee: eal Fae So 


= -0-28 [a | 26r, ma- E- 83) 
BIT 


(=) = 3p (1 — 2p) s fa [aoe + T? — 3 (ve + T3) — Yz Ts + 
a/u ; 
+ Byte (ve + te — aa] FO, 1) (E — 7) E — E) 6-4 


with 


706 J. Morkowski 


2 
(+) - (3) B fa apr T? — 3 (v3 + 12) — 2, Ts + 
Th } II 2 


+ 6 pt, (Ye + Tr. — fs) [v2 4 t? — O 42.22) 
+ Apts (Vz + Tz — Mz)] F(v?, T°) 6 ((4 — 9) - (4 — 2) (3.5) 


and, by definition, 


1l l 
F (92. a) = ot Boy? + t — p’) Po il {sai E: 1 + 
1 1 
+ ett bar __ T SF l 4 eet Bur __ 1 eft Bwr* — 1 ? (3.6) 
3 J 8 i 2 

= 51,5 RO nye 

ns AS?” > (3.7) 

Jef pL = TE (3.8) 


The quantity (1 — 2p)? (1/rs), is of special interest, as it may be interpreted 
as an inverse relaxation time determined by exchange interactions only. Indeed, it is 
the sole contribution to (1.1) which does not vanish in the absence of pseudo-dipolar 
coupling (i.e., in the limiting case p = q = 0, coresponding to D = 0). Integration 
over all the directions of 2, leading to the mean value (1/t,), of (1/t%);, is performed 
with use of spherical coordinates with polar axis parallel to the vector » + 7; then, 
in a similar way integration over the directions of » and 7 is performed. The result 
of the calculations is 


(5), 0-20 apes i5 ela ba bed (39) 


Tp 


n = l, 2, 4 for simple, body centered and face centered cubic lattices, and 


1 
1 Migs 1 5 
ant farl say aa +7) 
0 


sam EA ai - [Z(e + oy) — Z (e + ah (3.10) 


L 


Sn Eea eaa) 


o ar 


5 i ; 
~ eGo. e S21 (€ + 0) + 52; (e + 9) -= 5 ö Za (e + o) — 


re 2 Za (e + 0) Z, (e+ a}, (3.12) 
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with the definition 


= i 
S2,1 (x) = > cP HO 3.13 
an BIER Co) 
14 
and abbreviation 
o =f wu’. (3.14) 


For u œ 1/6, not very low temperature (say kT/JS = 1071) and for usually employed 
magnetic field strengths (of order of several thousands of oersteds) the relation o < 
e <1 is fulfilled. By means of asymptotic expansions given in the Appendix, we 
may approximately evaluate (3.10)— (3.12) for o<e<]1 obtaining the result 


Dog -+ og toys > (—2ine +Inte + ..).- (3.15) 
Finally we get thus 
1 BN OF\2 3-2 ZUB H 2 (Ge $ 
(=), = (1 2p) (27) n h 52 (ua) JS x 
kT kT ne 
2 CERAS 3.16 
«(i hog + 2m gia) Set 


As can be easily seen from this formula, exchange interactions lead to long spin-spin 
relaxation time. This result is consistent with Dyson’s (1956) conclusion of weak 
interactions only between spin waves (in absence of non-isotropic coupling). 

Evaluating of (1/17);, and (1/t;)z7 we confine ourselves to small values of 
- ù only. We transform the integration variables in (3.4) and (3.5) by the substitution 
p= ù+, T =u +T and expand the function F (»?, t?) = F (v2 + Qu.” + we, 
v2+ 24.7 + u?) into powers of (2u.9 + u?) and (2u. T + u?) respectively, 
retaining quantities up to the order of u? only. We average then over directions of 
u (we denote by (1/t,);, and (L/T pir the averaged values); the values of the inte- 
grals over » and 7 are found by introducing polar coordinates, whereas the radial 
parts are integrated by the method mentioned in Section 2. The results are 


a ~ (2) 3p (1 — 2p) Bu? (bw)? 2 (2) + 


Th ian 


T , S2, (£) — 5 Zp (e) Z o] ’ (3.17) 


‘ 


1 in eB a 
ae a T (27)? p? Bu? (Boo) la Z3 (e) + 


Tp 


G28, tc ge lO 
4606) + 2 510-2 4040}. 6.18 
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ONO ee 


where 1/7, is the relaxation time for spins waves with zero wave-vector, already 
calculated in I (compare eq. (6.6) and the next ones). 
In the limiting case of s <1 we have (see Appendix) 


(+) = =~ ae A p(1 — 2p) Bu? (Bw)? e (m? e — 21n e), (3.19) 
II 


Tyu 


0. 285 


= x £4 ont @) p? Bu? (fa) * in (3.20) 
mur To 


Tu 


The inverse relaxation time Mare may be thus given combining Eqs. (3.16), (3.19) 
and (3.20): 


l 
1 


SD 


= Z Huana ea (0. a "HE (FE) ~ 


kT kT 72 J 
2 ces Fal Apne Sa 3.21 
5 [m 2ug H ak e 5 E2) p hS 3 G2) 


kT 2 Da 
JS 2usB H i 


Taking, for example, the numerical values for nickel as in Section 2, we get 


1 T ; TaN ai 
so E10" Gaal + (ua) {18% 104 Soe a E (74% 10°) 


i$ leii ik : 
+ 2 In (74 x 108 a + 4.5 x 106 En ln (22 x 108 ny sec} 
(3.22) 


In usual circumstances and for u œ 1/6 = 5 x107%. 1/a terms proportional to u? are 


not important. Thus we get also a test of correctness of our approximate evaluation 
of lyre 


§ 4. Concluding remarks 


The resulting spin-spin relaxation time 1/7, = 1/7 + 1/7 is thus determined 
by Eqs. (2.9) and (3-21). For spin waves with 1/6 — as can be seen immediately 
from (2.10) and (3.22) — 1/r® is large compared to l/t? leading to Fee. 

At room temperature T = 293°K and in the magnetic field H = 5200 Oe the 
skin-depth 6 estimated from parameters determined by ferromagnetic resonance 
measurements in nickel crystals is ô = 5 x 1076 cm (Mitchell 1957); thus wa = a/d = 
= 7 X107? (a = 3.52 x 1078 cm). The inverse of the spin-spin relaxation time appro- 
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priate for these conditions of ferromagnetic resonance in nickel 1/t = 1/t® with 
u œ 1/6, is equal in order of magnitude? to 1.41078 sec7}. 

This value is smaller by about one order of magnitude than the ferromagnetic 
resonance line width 3.8 x 10° sec” 1 observed for nickel crystals in the same conditions 
(Reich 1956). It is thus evident that the spin-spin relaxation processes cannot determine 
completely the observed resonance line widths in metals. Also the temperature depen- 
dence of l/t is — at least for low temperatures — different from the well-known 
dependence of the line width on temperature. 

For a very wide temperature range, up to very low temperatures, the other impor- 
tant relaxation processes are interactions of spin waves with conduction electrons, 
i.e., magnetic interactions of ferromagnetic” d-electron spins with orbital moments 
of conduction s-electrons (Abrahams 1955) and exchange interactions between d-elec- 
tron and s-electron spins (Kittel and Mitchell 1956, Mitchell 1957). These processes 
give a relaxation time for nickel — in the conditions considered above — of the order 
of 6 x 10-8 sec and 5 x 10-® sec respectively, with a quite insignificant temperature 
dependence. 

Thus the apparently most important processes in metals, i.e., the mutual interac- 
tions of d-electrons (’’spin-spin relaxation’”’) and the interactions between d-electrons 
and s-electrons can jointly explain only about 1/10 of the observed line width (with 
rough quantitative explanation of the temperature dependence). It is possible that 
the assumption of ideal crystal lattice with identical spins at each lattice point is a too 
simplified model for discussing the relaxation processes (c.f. Van Kranendonk and 
Van Vleck 1958). The author proposes therefore to take into account the influence 
of the lattice imperfections on the spin-spin relaxation time. 

We remark here that our discussion is confined to the relaxation time for ferro- 
magnetic metals. In materials such as ferrites, accurate calculations will be very difficult. 
On the other hand, however, ferromagnetic line width — at least for some such mater- 
ials — can be satisfactorily explained by known theories initiated by Clogston et al. 
(1956). A 
The author wishes to express his gratitude to Professor S. Szczeniowski for his 


helpful advice concerning this investigation. 


Appendix | 


We quote here some useful asymptotic expansions of functions Z, (x) and Sy, (x) 
for argument values small compared to unity. Starting from well-known expansion 
of function x/(e* — 1) (compare for instance Fikhtenholts 1948, page 520) we get an 
Sree cise, A I ie See stl iNT 


3) It must be emphasized once more that our result is only a rough estimate of the order of magnitu- 
de, as 1/t depends on (D/ J)” and the pseudo-dipolar coupling constant D is a quantity no known precisely. 
The estimation D ~ J(g — 2)? which we employ here following Van Vleck (1956) is a sog orade one. 
Note Added in Proof. — According to the recent estimation of S. Charap and P. Weiss (Phys. eal 
116, 1372 (1959)) D/J should be about 0.06, giving 1/73) (and other expressions proportional to (D/J)*) 
greater by a factor about 9 than considered above. 
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expansion for Z, (x) = 1/(e* — 1). By means of the obvious relation d Z,,,, (x)/dx = 
= — Z, (x) expansions for n = 1, 2, 3,... can be thus obtained in succesion by inte- 
grating the series term by term. The results, valid for 0 < x < 2m, are 


Ze e =i De = amet, (A.1) 
Z,(x) =— nx + Ban xen, (A.2) 
1 2 2n (2n)! 


1 Gr Ea : 
AGO) aa en Cael (A.3) 


ZA =O- CQ) 2+ Sat — 5 tne y — 


Bon 2n+2 A.4 
Tee FDI" ; S 
where B,, (Ba = 1/6, B, = — 1/30...) are the well-known Bernoulli numbers. 
In a similar way, on account of the identity 
d 
© Sala) = — Zy(2) 2l) 
we get 
af 2 
S2,1 (x) = — 2¢ (3) — ¢ (2) In x + 5 |x—xlx+... (A.5) 
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RAYLEIGH’S RATIO AND TURBIDITY OF IMPERFECT GASES 


By S. Krericu 
Institute of Physics, Polish Academy of Sciences, Poznan 
(Received April 28, 1960) 


General expressions for Rayleigh’s ratio and the turbidity are given, containing 
the molecular constants Sa and S of isotropic and anisotropic light scattering. For 
imperfect gases, the constants S> and pe can be expressed as a virial expansion in in- 


verse powers of the molar volume. The second virial coefficients B's and B's’ of isotropic 
and anisotropic light scattering are calculated for certain molecular models of dipole 
and quadrupole molecules. The respective quantities are discussed and evaluated numeri- 


cally for the imperfect gases CO,, NH;, CHF, CH,;CN and COS. 


l]. Introduction 


Kammerlingh-Onnes’ equation, which is the equation of state of imperfect gases 
in terms of an expansion in virial coefficients, represents an important method for the 
investigation of intermolecular forces. In particular, the second virial coefficient B (T) 
not only gives insight into the parameters determining the central forces (see, Lennard- 
Jones 1924), but, moreover, provides information on the non-central intermolecular 
forces (see, e.g., Stockmayer 1941, Hirschfelder et al. 1942, Rowlinson 1949, Pople 
1954, Buckingham and Pople 1955 a). 

In recent years, the method of expansion in virial coefficients was applied to the 
theory of other measurable properties of imperfect gases. Harris and Alder (1953), 
Buckingham and Pople (1955 b), Hill (1958), and Jansen (1958) proposed virial 
theories of the dielectric permittivity of imperfect gases. Moreover, Buckingham 
computed the second virial coefficient for the molecular Kerr constant (1955) and 
molecular refraction (1956). 

Similarly, information on the nature of the forces acting between the molecules 
can be derived from an investigation of the divergence between light scattering in 
a compressed and in an ideal gas. To this aim, a theory of the virial coefficients for 
Rayleigh’s ratio S and the turbidity A of imperfect gases is proposed in the present 
paper. In particular, the second virial coefficient for light scattering on different 
models of dipole and quadrupole molecules is computed numerically. 

(711) 
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2. General formulation of the theory 


By statistical-mechanical theory of classical light scattering, Rayleigh’s ratio 
S and the turbidity A are given by (see, Kielich 1960) 


2 (n? + 2)? =F 2)? is anis 21 
S= gy (Sm + Sm (2.1) 
lôr (n? + 2) 1S anis 
a P 22 
h A re + DH J (2.2) 


wherein V is the molar volume of the scattering medium of refractive index n. The 
molecular constants SÌ and S#"* account for the molecular mechanism of isotropic 
and anisotropic light scattering. In the case of a scattering medium consisting of N mol- 
ecules of one kind, whose linear dimensions are small as compared to the light 
wavelength A, we have 


N N 
: Bat y am? dm 
isoeo æ 23 
gerdae at. ea 
#=1 J=1 : 
NaN : i ; : 
Br pace P A Im om? am? my 
h EE 4 
Sm = F578 Oa Ps fs For E CE D C4) 
i= 


wherein. m is the a-component of the dipole moment induced in the i-th molecule 
of the medium by the electric field F, of the light wave. Here, the brackets < > are 
intended to denote the statistical mean value, as defined by the formula 


U(r) 


S FEG ia 


fi fcafet EF GTi diya. CON 


k — Boltzmann’s constant, T — the Kelvin temperature, U (t) — the total potential 
energy of the system in the configuration T. 

For an anisotropic molecule immersed in the medium, m is given by the follow- 
ing expansion due to Buckingham and Stephen (1957): 


mO = ali, (Ey + FP) + 469, (Ep + FP) (Ep + FO) + 


(<X)= (2.5) 


apy 
Farem (Ep + FR) (E, + FP) (By + FP) + 5 (2.6) 
with a} denoting the polarizability tensor of the i-th molecule, and Bo}, and y@ — its 


sy ENER tensors as discussed by Buckingham and Pople (1955 c). The tensors 
Gop» Papy ANd 75,5 account for the properties of the isolated molecule, and are sym- 
metrical in all indices a, f, y and ô. F is the a-component of the molecular field 
at the centre of the i-th molecule due to the charge distributions of all the others, 
in the presence of the electric field E,. 
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The application of the foregoing theory to imperfect gases consists in expanding 
the molecular constants S's and S= in powers of 1/V: 


Be 


st 2 ak —S ES 
Ss 7 55 Pe 


AS (2.7) 
Bo Geo 


cans ae es 
Sat V i z 


TE (2.8) 


The coefficients of the expansions 41, 43is, Bis, Baris, Cis, Canis, |.. are termed the 
first, second, third,... virial coefficient of isotropic and anisotropic light scattering, 
respectively. 

The first virial coefficients, A} and AS, characterize light scattering in a gas 
whose molecules do not interact (ideal gas), and, with respect to eqs. (2.3), (2.4) and 
(2.6), are given by 


is 824 
Je aa a2 N, (2.9) 
anis 52m4 
Ags = 4544 (3048 hap — Xaa pp) N, (2.10) 


wherein œ = }@,, is the mean polarizability of an isolated molecule. 

The second virial coefficients, BÌ$ and B's relate to pairwise interaction between 
the molecules. The remaining virial coefficients, ae Gane etc., account for the inter- 
action of three, four etc. molecules of an imperfect gas. The coefficients yield a meas- 
ure of the divergence between light scattering in a real gas and in an ideal gas. We 
shall restrict ourselves to a discussion of the virial coefficients Bs and Bws only, 
which, for a not too strongly compressed gas, account essentially for the divergence 


from the ideal gas. By eqs. (2.3), (2.4) and (2.9), (2.10), we have quite generally: 


is 16204 «2 N nxt N2 y _ Ui ‘ 
Bs = — A B(T) + ARR ae 9a?) e FT dr;;dw; dw; 
(2.11) 


i 4 y2 X p: U AAi 
poun Í 4 (BV Pep— VQ p) e FT drjdojda;. (212) 


Herein, the tensor 


7 1 ({[dm® Im? am® vi 0 G@) G) » i 
yu. a a AA E ee a ECA, 
mids RE ES {( Eg 3 Eg 3 Es 3 Es (Hap Ayt Map ais) a 


and the second virial coefficient of the equation of state of the real gas (Pople, 1954) 


_ U 
B(Z)= we Í Ey (FT — 1) drj do; doj, (2.14 
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have been introduced; U;; denotes the total potential energy of interaction of the 
i-th and j-th molecules, oh face the distance between their centres, the variables 
w,;, @; describing their aedi orientations; Q = f dw; isthe integral over angular 
Daa 


3. First Approximation to the Theory 


In the first approximation, a molecule of the imperfect gas under consideration 
is assumed to possess the polarizability of one isolated, i.e. a polarizability independent 
of the surrounding, neighbouring molecules. In this case, the expansion of eq. (2.6) 


yields 
ore a, (3.1) 
and the virial coefficients (2.11) and (2.12) assume the form 
pees ccs Bie e fff — 1) dr; dw; doj, 
Bes = n | Il R A aus) en H do, dwj. (3.2) 


Hence, it is seen that the second virial coefficient of isotropic light scattering, B's, is 
accounted for directly by the second virial coefficient B (T) of the equation of state 
of the imperfect gas only if the molecules exhibit a polarizability that is unaffected 
by the presence of their neighbours. 

The total potential energy U;,; consists of the energy U? (r,;) related to the exist- 
ence of central forces of attraction and repulsion acting between the molecules, 
and of the energy V;; = u (r,;, @;, @;) resulting from non-central electric intermolec- 
ular forces: 


Uy = U (ry) + V, 


The energy V;; will now be considered to represent a perturbation in U (r;;), 
so that the virial coefficients B'$ and B®" of eq. (3.2) can be expanded as follows: 


Uer) 

i is , 8ata? N? 1 1 ae ij 

Bs = Er centrBs = SSeS j ab ab Vie at drij, 
n=1 


(3.3) 

pes _ 52N? SA 1 IRAT o n, 

S wae 7544 va mi — ar) (Gap Zap Viz — 3a? Vive At drij, 
FETA 


wherein ~ 
UO 


4 y2 N2 pias Bee S153 À 
aAa a ai ew fe (e T —1)drj, (3.4) 
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is the second virial coefficient of isotropic light scattering, dependent only on 
the model of central intermolecular forces; the horizontal line is to denote isotropic 


= (w) do; dw; 9 


ie. integration over all possible configurations of either molecule. 
The virial coefficients, as given by eq. (3.3), will now be computed for some models 
of the potential energy V;,; of electrostatic interaction between the molecules. 
Dipole molecules. The potential energy of electrostatic interaction between two 
dipole molecules is 
V; = wh pb TS, (3.5) 


wherein 4 is the -component of the permanent dipole moment of the i-th isolated 
molecule, and 


$ l 
T =] $ {ri Oap — ija TRS iA], (3.6) 
1J 


is the tensor of dipole-dipole interaction. 
Substitution of (3.5) in eq. (3.3), with the isotropic mean values V7 and a) oh Vi 
of Appendix B, yields 


is is 8t q? pA N? ah < 3y4 =12 
dipBS = centrBs + 372 T fe a IST? KY SE ee (9 (3.7) 
nis 2624 N? =e 1204 —12 
dip BS = 3375 Ja aT? (3Xa8 [la Mp — Qaa MB)” | o ga STD Tory 
(3.8) 
wherein : 
uir) 
<r = f pee ar an (3.9) 


is a quantity that depends solely on the model of central forces. « 
In the case of a molecule having the axial symmetry, with the axis of symmetry 
directed along the 3-axis, the following tensor components of u, and yg are non-zero: 


Mg = os Qu = Mag =ð]; X33 = Qi; 


hence, 
3093 He Hp Xaa u3 = 6a 6, K’, (3.10) 
wherein the quantity 
gol ae Reais aN (3.11) 


Se ga ay 2% 
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accounts for the anisotropy of polarizability of the isolated molecule. 
By (3.10), eq. (3.8) reduces to 


anis 104rt æ? oe uNe =6 yt —12 
dipBs = EE Dit eT? Lip) ee fs (3.12) 


valid for axially symmetric molecules. 

We shall now proceed a step further, and shall take account of inductive inter- 
action between the dipole molecules. Now, for anisotropically polarizable molecules, 
we have (see, Appendix A) 


Vi; = u® ee TEP — 4 (af? pus Q) 4 cos pu ®) Te TD edie: (3.13) 


wherein, in addition to the first term accounting for dipole-dipole interaction, terms 

determining the inductive effect of dipole-induced dipole interaction appear. 
With the previous expression for V;, the virial coefficients are now 

l6z* g? u? N? 


iSS is =g 2u? 
se ae hake (e< > E Ta TAA ee 


+ 2 (Bap Ha bp — Caa Y] <r > + at (3.14) 


Bas anis 5274 N? 2 —6 
S =dpBs + 3375 2 T (3008 Ma Mp — Xaa Mp) \(3%y5 6 — Ayy Gas) <tyz >) + 


12u4 
os agate [4 (3a y5 Zye — Ayy Mae) Ma Me + (3Xy3 My Hö — yy M) Xes + 


+ (8ay¢ Xy — Ayy 5s) ue] ire a + ah $ (3.15) 


wherein the quantities dip?’ and dips are given by eqs. (3.7) and (3.8). 
By eq. (3.10) and the relations 


== 2 
3a ap Aap ra X za Q ge = 18 a? Ons 


(Batgy Wey — Can Kp) Hp Hy = 6a? 4, (1 + 28) H”, (3.16) 


eqs. (3.14) and (3.15) yield, for axially symmetric molecules, 


is is AIOT a? uN 2 
Bs = aipBs + mes i rae ker »+ Ta T (5 + 462) Cras =r a (3.17) 


20874 æ? 62 u? N? 


BS" Bs + eee a eee = AEREO =! 


(3.18) 


= Te 


Rayleigh’s Ratio and Turbidity of Imperfect Gases TIT 


whe the additional contributions to B§ and pani from inductive interaction between 
dipole molecules have been computed. 


Quadrupole molecules. For non-dipolar, anisotropically polarizable molecules 


possessing a quadrupole moment given by the tensor @,,, the potential energy of 
electrostatic interaction is (see, Appendix 4) 


V; = — 40%) OY Bon =. n (ats 0; J) op + AP Ov 0) TH T Gi) AL t (3.19) 


ayd ~ Ben 


the first term accounting for the energy of quadrupole-quadrupole interaction, whereas 
the second determines that of quadrupole-induced dipole interaction. 
Substitution of the foregoing in eg. (3.3) yields 


mo ere NE Ee T 
abs = = centrDs + E Ou Oup fisa Cr IN a ayy, Ov Cri os Sab ne 


(3.20) 


on 20824 N? Ee 
auadBS = 479512 kT (Saag Day Ogy — Xaa Opy Opy) Gan Xoe — Xas Kee) Tij $ a 


D (30 Oon On — X66 Oen Oen) ai °°) -+ p 0 (3.21) 


Herein, the terms in (75 *) result from quadrupole-induced dipole interaction, 
whereas those in <r; 10, _ from quadrupole-quadrupole interaction. 
In the case of leer exhibiting the axial symmetry, we have 


Osp Ou = 36°; 3 aT: OD, Op, ~ Kaa Op, O py =4& Ôa 6, (3.22) 


and the virial coefficients of eqs. (3.20) and (3.21) reduce to 


824 %2.@? N? Ss OA ares. , 

OES ae Sarma fisa Cra D + gr sE t \ (3.23) 
. is 20874 a? On Q02 N? ô —8 20? —10 (3 24 
quaaBS” = 475948 ET Ba a LT > + Gap <TH" ae © 24) 

with O = O, = —2 Oy, = —2 O» denoting the quadrupole moment of the axially 


symmetric molecule. © 
= Polar molecules. For molecules exhibiting both a dipole and quadrupole mo- 
ment, the potential energy V; is given by (see, Appendix A) 


V; = = te UP e TE H( uË oy ot uP) WA, = + OÑ) oR TS AP oes (3.25) 


wherein the terms account for dipole- -dipole (vj =3), dipole- eas (r°) and 
quadrupole- -quadrupole (735°) interaction, Fenech: 
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With the foregoing expression, and by the formulas of Appendix B, eq. (3.3) yields 


Ome a ly” 6 ee 
= = dipBs “iF eee fu Op Oas Cis S TER (Ozg Ha Ls)? [Tij > a 
(3.26) 
: T 20874 N2 bs 
BS" = aipBS + TAAT fauua [p—an 113) (3&3 Oys Oss—&yyOs Oa) Lr > — 


8 a 
— app las Pay Mp by — Xap Pap Hy — Xaa Ogy He py)? Cis Hy — 
2u” Xal Oas 2 A 3 PT 
E 5ĘT [3u Xy O; SF 4 (3z O yp E m Ayy Os.) Me Hel Crij N 4 e.s.f 9 ( . ) 


wherein the quantities 4,,B and 4;,,By"* are given by eqs. (3.7) and (3.8). 
For molecules exhibiting the axial symmetry, 


3 tap Oay Hg My — %ap Pap My — Zaa Opy Mp My = 3% Ô H? O, (3.28) 
and hence the virial coefficients of eqs. (3.26) and (3.27) reduce to 


87t a? u? O? N? Sfi | 
ea {o> —% 54T r< at ii 


20874 a2 62 u? O2 N2 
PARAE TA N AATE kra’ — 


BS = dipBS + (3.29) 


= dipBS” “fe 


17544 k? T eae rg > + ah (3.30) 


4, Further approximations to the theory 


Further approximations to the theory should account for the fact that, in a suffi- 
ciently condensed gas, the polarizability of a molecule is affected by the presence of 
its neighbours, and that the molecules exhibit the effect of hyperpolarizability. Bucking- 
ham and Stephen’s expansion as given by eq. (2.6) should now be applied, yielding 
the following expression for the differential polarizability of the i-th molecule in the 
absence of an external field: 


am Phe ie capes? 6 
( 3E, f a fe + BF + 5 sa FP FS + a lon g- ani 
5 aA 


(4.1) 


Neglecting herein the terms in =>, the virial coefficients of eqs. (2.11) and 


OF 
dL,’ 
(2.12) are obtained in the form 


i l6zc4 @? N ; 
Bs =— rg ars B(T)+ same fegt Pr EY + By Fy + 


+E YD FOFD + 47D FORO HL en ia dr,,do;do;, (4.2) 


— 
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Bais 5224 N? i i) i i i f. ] 
Beg Pa | ff 2B 2b — 28 off + 303 + 0) 6B, FP + 1B, FP + 


+4 Np ee joy FË +4 Mele F9 FY) +...) — (ai 4 xt) ( ) (63, FO + By, F9) =e 


+4 (i) FO FD + 4 (7) FO FD + -Wg . dw. dw: 4.3 
YBByò YB6yvö +...) +..}¢ r,,dw,dw;, (4.3) 


The foregoing expressions will now be discussed separately for quadrupole and 
dipole molecules. 

Quadrupole molecules. Since „gy = 0 for non-dipolar molecules, and since the 
electric field at the centre of the i-th molecule due to the quadrupole of the j-th mole- 
cule is (see, Appendix A) 


Fo = iT OB. a) 
eqs. (4.2) and (4.3) yield with respect to (3.19) 
a a N? 


BS = as anig ae TTO VaappOrs O; ar » at as (4.5) 
is so lrt N? 
Bs = quaa BS SF AF {7 (S&a Yapyy — Xaa YeBvy) Oss Oos + 
Sr AY aab (3a,¢ Oy. Ose — Ayy Os. O5.)} Ca F760 2 (4.6) 


the quantities quadBs and Bee being given by eqs. (3.20) and (3.21). 
In the case of axially symmetric molecules, we have the following non-zero tensor 
components of 7,,4 (see, Buckingham and Stephen, 1957): 


Vii = Yz2 = 3V2 = Vp» Yus = 72233 = 3 (Yi) +Y) Yss33 = YI 


whence 
3 ep Vary — Taa Vapyy = 21% Y Òa by, (4.7) 
the quantities 
b= MEM, y= (i+ ad (4.8) 


yielding a measure of the anisotropy of hyperpolarizability and the mean hyperpolari- 
zability of a non-dipolar molecule, respectively. With respect to eqs. (3.22) and (4.7), 
the virial coefficients are given by the expressions 


‘ : O2N* 32 
BE = queaBS + OE Cri) Ho (4.9) 


l ay ô, (20 + 496,) OPN? , 
pet — qhaBs e Samt ay d 20 + PA) N CrO epig (4.10) 
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which contain terms accounting for the effect of the hyperpolarizability of the mole- 
cules on light scattering in an imperfect non-dipolar gas. 

Dipole molecules. For simplicity, the terms relating to the tensor y,,,, will now 
be omitted in eqs. (4.2) and (4.3) and the problem restricted to one of dipole-dipole 
interaction. In the present case, in addition to V; as given by eq. (3.5), account should 
be taken of the molecular field F® given rise to dt the centre of the i-th molecule by 
the dipole of the j-th one: 


FO = — TH uP; (4.11) 


thus, by the formulas of Appendix B, eqs. (4.2) and (4.3) yield 


is is 6404 au? Baap Lp N? —6 —12 
Bs = aipBs 1 O10 Tae Cie Dyas ae Chip S Ffa (412) 


20874 N? 


pais zA Bai a eee ISALT com Busy — Xaa Beby) Hy us G mes 


6u4 l f 

+ SELITA 252 T2 [rij = Js a A 35 (3008 Ma Ug — Xaa Hp) Pyy Lo (<5 cis 
12u4 e 6u? 

ae “ans T? mj) Pie J A ae T2 (3Xap Ma Hp — Xaa LB) (3B y se Uy — 


— Byys Me) Ho be <r > + ia (4.13) 


the virial coefficients qip Bs and Bira being given by the expressions of eqs. (3.7) 
and (3.8). 


For molecules exhibiting the axial symmetry, 


Ha =M, Ps3 = By, Bis = Boog = By > 
whence 
3 Gap Pegy My — Zaa Pegy by = 18 ap Ôa Og H, 
(3 Bapy Ha — Baap Hy) Mp liy = 6 BS, u’, (4.14) 


wherein the quantities 


dai -ÎL ft, p= 5 Out 20D (4.18) 


yield a measure of the anisotropy of hyperpolarizability and the mean hyperpolariza- 
bility of a dipolar molecule, respectively. 
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Substitution of (3.10) and (4.14) in eqs. (4.12) and 4.13) yields the virial coefficients 
for axially symmetric, hyperpolarizable dipole molecules as follows: 


64rt afu? N? È 
SAAT koz" yet ae 95h2 T2 tie y at (4.16) 


Sais anis 4167* aĝôx u? N? fs 
Bs = dipp S 1 ny a + 2956p) CTi; 55 + 


is is 
Bs =aipBs + 


out 
ay moe qa (2 + 51 ô) <r 12S a (4.17) 


Dipole molecules of variable polarizability. In the present section the terms 
OF 
previously neglected, relating to the derivatives Ae of the molecular field and yield- 


ing additional contributions to the constant polarizability a,, of the isolated molecule, 
will be considered. However, in order to avoid further complicating the calcula- 
tions, the hyperpolarizability terms of eq. (4.1) will be dropped; Kirkwood’s method 
(1936) now yields, for anisotropically polarizable molecules, 


dm? 
( = So) — a) TH ae) + PR BSE pete (4.18) 


and the virial coefficients of eqs. (2.11) and (2.12) can be represented in the form 


is 164 a?.N 16a N? D i 
np = SEEEN peryg BAN fi (oat 9 


Eu a% aus Ty jet je T ir dr;; da; do; , (4.19) 
penis 5274 N? i ( i) .@  .@ OW) mi) 
pon = S220 IT (alt af — hoff — 2608 a — a 8 TH? — 
Uy 
=] 26a, ee a?) aG) of TP rei cee an dev; deny. (4.20) 


Substitution of (3.5) yields 


164 a N? 
BS = aipBS + gaT (Ot — 
— (Bttap Ma Hp — Xaa Hp)} (B%ya My Ha — Ayy u3) rap bones (4.21) 
104a4 N? 


B = gipBs + IOIA kIT e (3tap Gap — Xaa % pp) (Ayo Hy Mo — Xyy Ho) H? + 


ar [30% p? = Kaat la ug — Qaa up)l [3 Bays rye — Ayy “tt Ho hs — 


— (3a y5 Xyð — Ayy xaa) Ml} STi > + (4.22) 


722 S. Kielich 


Hence, by eqs. (3.10) and (3.16), we have, for axially symmetric molecules 


644 x? 6, (10 — ôa) ut N? 


PRAT Cr > +s (4.23) 


is is 
Bs = aipBs + 


41674036, (5 + 46, + 962) tN? , —9 
Sisal Te <i dear 


BS" = aipBS + (4.24) 


5. Discussion and conclusions 


In computing the virial coefficients BS and B2"* in sections 3 and 4, no assump- 
tions were made as to the form of the potential energy of central interaction between 
the molecules. It will be remembered that Lennard-Jones (1924) proposed the follow- 
ing general expression for U® (r,): 


U% (rj) = — 2 aF tn, m >n, (5.1) 
io Tij 


the first term accounting for the energy of attraction and the second — for that of 

repulsion between the molecules. The constants A, and A„ and the exponents n and m 

can be determined experimentally. i 
To simplify, the Lennard-Jones potential will henceforth be applied in the form 


(see, Pople 1954) 
n\n Om 
U(r;;) = 4e 1 — (2) \, (5.2) 


wherein e represents the negative value of the lowest potential energy and To is the 
intermolecular distance r; at which the attractive and repulsive energies exactly 
balance. 

By the Lennard-Jones method (1924), with respect to eq. (5.2), we have (see, 
Buckingham and Pople 1955 a) 


co Ul) Crj) eae 
f (FF — 1) dry = EF (Ha) — 28a) (5.3) 


co UO (..) 
SF =" Sg OA nr- 

lon Sm f ite ary = ET HG), (5.4) 

0 
wherein 
a7—n a 6m +n— 3 a 4e 
H, = 6 Bees eee oe bpieet & a Oe a 
Nar ea is a E 6.5) 


Rayleigh’s Ratio and Turbidity of Imperfect Gases 723 


are functions introduced by Pople (1954); these have been discussed and their values 
tabulated by Buckingham and Pople (1955 a). 

Eqs. (5.3) and (5.4) can serve for putting the virial coefficients B$ and pans 
as calculated previously, in a form adapted to numerical evaluations, provided the 
central force parameters € and rọ as well as other molecular parameters, such as & 
and u, whose values for some gases are assembled in Table I, are known. 


Table I. Numerical values of the molecular parameters and of Ae and ew (A = 4358 A) for some 
imperfect gases. 


Gases NH, 
u xX 1078 e.s.u. 1,47 
ay X 1024 cm? 2.42 ) 
ay x 1024 cm?® 2.18 © 
a x 1074 cm? 2.26 
Og 0.04 
Ta X 108 cm 2.60 5) 
elk °K 320.0 5) 
is | cm? 
Ags — |X 104 6.62 
mol 
0.03 


1) Stuart, H. A., Die Struktur des freien Molekiils, Berlin, 1952. 
2) Le Fevre, C. G. and Le Fevre, R. J. W., Revs. Pure and Appl. Chem. 5, 261 (1955). 

3) Computed from polarizability of the C = 0 and C = S bonds. 

4) Johnston, H. L. and Mc Closkey, K. E., J. Phys. Chem., 44, 1038 (1948). 

5) Keys F. G., J. Amer. Chem. Soc., 60, 1761 (1938). 

6) Buckingham A. D. and Pople J. A., Trans. Faraday Soc., 51, 1173 (1955). 

7) Lambert, Roberts, Rowlinson and Wilkinson, Proc. Roy. Soc., A 196, 113 (1949). 
8) Landolt-Bérnstein, Physikalisch-Chemische Tabellen, Springer. 
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(i) Quadrupole molecules. For anisotropic molecules possessing a quadrupole mo- 
ment, eqs. (3.23) and (3.24) yield for the Lennard-Jones model 


A 52 O? N? f60a 10? 
quaaBs = contr BS nt ey oz H; (y) + Fa Hio o), (5.6) 
is 26r a? 6g O? N? \6a0q 8% \ 
a EEEE A 


; 5 g2 2 - 16x4 a? N 
centrBS T aa {He Oiz 2 Hys (y)} S ak eee B® (T). (5.8) 
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Assuming a value of O = 5.29 x 107°% e.s.u. for the carbon dioxide molecule 
(see, Buckingham and Pople 1955 b) in addition of the data of Table I, the foregoing 
formulas yield, at T = 300°K, 


; P mals iat 
quaaBg = (31.79 + 1.86 + 9.59) x 10-2? Z = 43.24 x 107 20 = 
Banis — (0.14 + 14.4) x 1074 TÈ Lobi = 
quad” S SA mol? mol? ` 


Thus, in the case of isotropic light ZE the contribution from non-central 


intermolecular forces, (1.86 + 9.59) x 107 22 


is seen to amount to 26% of the 


T 4 
quaas, one fifth thereof resulting from quadrupole-induced dipole 
interaction, and the remaining four fifths ~ from the quadrupole-quadrupole interac- 


total value of 


tion predominating in the case under consideration. The second virial coefficient of 
anisotropic light scattering, PON a yields a contribution that is barely 0,4% of that 
of uadB5» and so plays practically no part at all-here. 

By eqs. (5.3) and (5.4), the effect of hyperpolarizability on light scattering as 
given by eqs. (4.9) and (4.10) assumes the form 


R is 80x5 ayO? N? 
Bs = quaaBs + “wi Ay (y} + «+5 (5.9) 
oy 


5275 wy5_ (20 + 496,) O2 N2? 


pa eae Bo = 31544 yt isp (y) + eee (5.10) 


Assuming a value of y = 25 x 107*8 e.s.u. for the CO, molecule, eq. (5.9) yields 


Be abs ae LU 2° “ai 
thus, the effect of the hyperpolarizability of non-dipolar molecules on. B} is quite 
unimportant in the present case. 


(ii) Dipole molecules. By (5.3) and G 4), the expressions of eqs. T 7) and (3.12) 


. now assume the form 


is me a 7AN“ SNE s3 7 
dipB'S = centr Stab ease 1844 r3 e? 2 | H, (y) = sot 5 H (y) +.. i 4 (5.11) 
anis dam? a2. UAN? 
~ eis T RRE i He 0) + soa Hiab) +» |. - G12) 


The expressions of eqs. (3.17) and (3.18) which, in addition to dipole-dipole 
interaction, account for dipole-induced dipole interaction, can now be written as 
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follows: 
a A Ag a? 2 N2 4,4 
B Z i 1s 4 4 X 
S dipB S 2 344 r? ey? g (y) =f 200r8 62 (5 IF 462) Hi (y) Sr a ? (5.13) 


anis anis 
Bs” = dipBs a - 


S27 Oey 


SIOAR TOE. 


pA spe 
de Hey) + HI 7 + 1189) Fh) +... 


(5.14) 


Table II. Theoretical values of virial coefficients BS and Bow of light scattering, in cm5/mol?. 


central-forces dipole induced dipole-dipole 
potential go pac interaction al 
action 
NH, Bs X 10? 7.36 6.63 32.11 46.10 
T= 320°K BŞ” x 104 0. 0.003 0.45 0.453 
CH,F Bs X 10? - 9.06 4.53 20.83 34.42 
T = 320°K Boae 0 0.02 3.04 3.06 
CH,CN BË x 10? 153.72 67.64 707.11 928.47 
T = 400°K BE! x 10 0 |l w% 113.0 113.20 


In Table II, the numerical values of the virial coefficients BS and BÈ computed 
from the foregoing formulas for the dipole gases NH}, CH,F and CH,CN are asembled. 
Each of the respective values is split in three parts, the first accounting for the contri- 
bution from light scattering due solely to the presence of central forces. The second 
and third parts are contributions arising from the fact that, in addition to central forces, 
also forces of electrostatic i.e. dipole-induced dipole and dipole-dipole interaction 
are present. The data in Table If prove immediately that the principal contribution 
to the total value of B} is related to the non-central forces, the inductive effect playing 
a lesser part as compared with the considerable dipole-dipole interaction. As to the 
second virial coefficient of anisotropic light scattering, its contribution to the total 
is negligible. The latter conclusion was reached earlier by Benoit and Stockmayer 
(1956) in their -theory of the degree of depolarisation of imperfect gases. 

Eqs. (4.16) and (4.17), which accounted for the effect of the hyperpolarizability . 
of the molecules on light scattering in a dipole gas, now assume the form 


ORE 
104r a Bbq, u? N? 
112544 rè £y? 


E N 


S $ 167° 3 N2 3u4 yt 
ens e fa, 0) + spore Ha 0) + el, (5.15) 


B — aB + la + 2564) He (y) + 
š 3y 


D lt 5.16 
© + ee (5.10) 
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To exemplify, in the case of CH,F, with the data of Table I and a value of B 
= 10 x 10-® e.s.u. (see, Buckingham 1956), eq. (5.15) yields 


i is = cm? 
BY = 4,B$ + 205 x 10°? 


hus, the contribution due to the effect of the hyperpolarizability of the dipole mole- 

Bi = 34.42 x 107? e This 

points to the possibility of a quite considerable effect of the hyperpolarizability of 

the molecules on light scattering in a dipole imperfect gas, and the situation resembles 
the one found to exist in a dipole liquid by Buckingham and Stephen (1957). 
By eqs. (5.3) and (5.4), eqs. (4.23) and (4.24) can be rewritten as follows: 


cules is seen to amount to 60% of the value of qip 


An" or? 8 (10 — ôa) u4 N? 


Bs = dipBS + — oF Ne 25 J4 76 g 78 z2 H; (y) + ...} (5.17) 
7 anis , 207° a? 6, (5 + 46, + 963) u4 N? 
po eines, i a oe ee £ $ 
S =dpBs + 1125 4175 e? H, (y) + -- (5.18) 


The contribution to isotropic light scattering resulting from the effect of the 
molecular field on the polarizability of the molecules is very small, amounting to 
barely 0,3% in CH,F, and to 1% in CH,CN. 

(ii) Dipole-quadrupole molecules. With the Lennard-Jones model, eqs. (3.29) 
and (3.30) can be rewritten as follows: 


až a? u? O? N 
BS = aipBS + eee n0- 2 5 re ~ Has (9) + iF (5.19) 


ma 1325 a2 ô? u2 02 N? 
BS = E ee mo-i zas tn) + af ee 


the quantities 4,,B§ and 4;,Be"* being given by eqs. (5.11) and (5.12). 
For COS, we have u = 0.72 x 10-18 es.u. and O = 2.9 x 10-28 e.s.u. (see, 


Gordy et al. 1953) in addition to the figures of Table I, and eq. (5.19) yields, for 
6S A 


Bit = (189.5 + 1.21 + 2.48) x 1072 _, 
mol? 


Here, the contributions of dipole-dipole interaction (the second figure in the brackets) 
and dipole-quadrupole interaction ae third value) to Big are very small, as compared 


with the value of 189,5 x 10- Fee 


mË resulting from the central forces alone. Clearly, 


this is due to the small numerical value of the dipole moment of the COS molecule. 
In the present case, dipole-induced dipole interaction plays a greater part; by eq. 


Rayleigh’s Ratio and Turbidity of Imperfect Gases TOA 


; 5 

(5.13), it contributes to B$ to the amount of 3.95 x 10-2 a . The same is true 
mo 

of quadrupole-quadrupole and quadrupole-induced dipole interaction; by eq. (5.6), 


5 
f ETRA cm 3 SRM R 
their contribution is 6.5 x 10-2 E: Thus, the entire contribution resulting from 


5 
EA : cm 

electrostatic interaction forces amounts to 14.14 x 10-2 2? 
mo 


which is 8% of that 


of the forces of central interaction. 

The immediate conclusion from the examples discussed above is that, in principle, 
divergences between light scattering by an imperfect gas and by a perfect gas should 
be expected to reside only in the isotropic light scattering, accessible to investigation 
in the quantity B's. For molecules having a constant polarizability, B} is given directly 
by the second virial coefficient B (T) of the equation of state of imperfect gases. In 
the general case the mean hyperpolarizability of the dipole or quadrupole molecules 
can be determined from the experimental values of B} and B (T). Hence, parallel 
theoretical and experimental investigation of light scattering in imperfect gases should 
provide new and valuable information on the intermolecular forces and, moreover, 
on the effect of hyperpolarizability of the molecules. 

The author wishes to express his great indebtedness to Professor Dr A. Piekara 
for his interest in the present investigation and for his helpful discussions. 


Appendix A. 


Potential energy of electrostatic interaction between two molecules. The poten- 
tial of the field at the centre of the i-th molecule due to the electric charges of the 
j-th molecule can be represented by the expansion 


ei) 1 l ¢ l Lie l 
P= ae — wp? Va (2 +5 deb VaVe (+) SA by Ve Ve Vy i. + wees 


z 
(A.1) 
wherein the quantities 
M =v, 7 
n 
=D P, (A2 
n 


ely Bs j) (4) 
B= Dre 12 1 
n 
Se 
ID, DPD QD, 
n 


account respectively for the electric unipole, dipole, quadrupole, octupole, 2*-pole 
moments; e — is the n-th electric charge of the j-th molecule, r9 — its radius 


ð 


vector, T = 1; — the vector connecting the centres of both molecules, and 7, = on 


Fi 
— Hamilton’s differential operator.. 


728 S. Kielich 


By (A. 1), the potential energy of electrostatic interaction between two neutral 
molecules is obtained in the form (see, Pople 1952; Jansen 1958) 


Vy = pe uP TP +4 (WO GP — GQ uD) TH) + 
+ dy (2UP IP, — 3qG} a + 21S, 09?) Tiyo + --- ven 


wherein the tensors 
= J j l 
TP = — Va V6 (>) » Taby = — Va VeVy (>) > 
1J tJ 


TO = = Vag Vals ae (A.4) 


Tij 


account respectively for dipole-dipole, dipole-quadrupole, quadrupole-quadrupole etc. 
interaction; these are to be found in explicite form in Pople’s paper (1952). 

With regard to Laplace’s equation, the following tensors of the quadrupole and 
octupole moments of the molecule appear (see, Buckingham, 1959): 


OG = + Gai — Ivy ap)» 
2, =h (519, = ID, Op, T 1s bya T Liss Oap)> (A.5) 
and eq. (A.3) can now be rewritten in the form 
Vy= Wu TP +10 OP — O uP) TH + 
+ Bu? QD, — 509 0 + 300, uP) TR. +, (A6 
where the multipole moment tensors are symmetric in all indices. 
By the expansion of eq. (A.1) and the definitions of eqs. (A.4), and (A.5), the 


electric field due to the charge distribution of the j-th molecule at the centre of the 
i-th molecule is 


FO = — 7,9 = TË WP £372 OY ATRO +... (AD 


For polarizable molecules, account should also be taken of the potential energy 
related to the effect of induction; this energy is, quite generally, given by 


Hip = — 409 FP EP + ofp FO FY, as 
Substitution of (A.7) herein yields 
p S DP wh? + 0} uP uP) TE TH + 
HERP OR + af POD TY TR- 
— 18 (aaa Oy8 O + aah Oy OS?) TSR They + 


+ ie (OG) KP OD + of? uO OD) TE TOD — + ... (A.9) 
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Appendix B. 


Isotropic averages. With the potential energy of electrostatic molecular interac- 
tion of eq. (A.6), the following mean isotropic values result: 


wie @) aC D 
V,=0, af} a} Vi; = 0, 


= 7 Wey DE ty alii l K ; 
Vi = = ut TEP TEP + Tye Ors Oas Tors Toes + ggg (Oan Oap)? Terre Torre» 
V3 = — se up Ha H? (TEP TP Tee + 2 Te? Tove TD), 
nos 1 = i i Gi ij j 
Vij = — 5g (Oep Ka te)? ES ES SE a PE IS 


V = ae Mè (TÈ Te? Tre? Trg? + 2T ex) Ter’ Tre? Tea) 


l 


O a Vi = yaar (450a? ut + (Bap Ha pip — Hen H)*} TE Tov + 


Xab Xap 


2 
T T75 {31502 u? Ogg Oag + 2(3%ap Oay Ogy — Xaa Opy Ogy) (Xoe Mo He — 


ee ME 1 ; 
— aoo M} TER TE? + gg 410 (ap Gap)” + BBa Oar Oz, — 
— Aaa Ogy Oz,)*} T T ? (B.1, 
Se E T l ‘ 
a aB Vig = — 71095 {441a (Oag Ha Up)” + 8(3%ap Oav HB My — 


Gj) mi) pmj) (ij) plij) mG) 
— %ap Oap by — bea Ogy He uy)’ TA Te Tovrea + 2T cr Tove Trvo} <a 


iF ae 7 [Baas Oas u? T 4 Gaas 0, — Ara Ogy) Ug Hy] Qis Oas p? = 
— 32 (Betas Ony Hp Hy — ap Pap VŽ — Fea Oy Mp HV?) TEP Ton Tone 


a C TOD POD PG GD Gi) pi) pë 
ava GR (TP To? Tee? Tre +2Tor Ter’ Tre Te) + 
4 gfs Gag Ha Hip — Cae UA)? (10 TE TH? T Tas) — To Toe T Tre): 


The products of the tensors TED, TE and 7 are: respectively, 
unter’ gi ROE: e A ERTE 
TOTO = or a T T S 


TO TA = 9075, T, Te = 2520r i, 


OTY otre 


TD TM TD = — 1215", TEP TEP Tore Rog (B.2) 
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With the foregoing values, eq. (3.3) readily yields eqs. (3.7) and (3.8) or (3.26) and 
(3.27). 

With the dipole-dipole interaction energy of eq. (3.5) and the molecular 
field in the form (4.11), the following non-zero mean values are finally obtained: 


lise BO) FOV, V = —? u? Beat Ls a 
BO, FP VE = — H 1° Baas Ho TG” 
on poy ya > 3 U? ap Papy My Tj 
a} BY, Fe V; = E 336 Baap Heg We + (3e,,5 By Ks — Z yy pA} 1 ij z (B.3) 
at) Ba pa, FOV? Vie= = a5 M’ Xap Doy Hy T, ij Piei 
@ 80) FO 73 — 24 r712 — 8u? 242 
Xah Baby Sty = EGT Xu 6 Baap pTi 1225 { Id Peab He + 


+ (3bapy He Ug Hy — Baay lig LÈ)} (3056 Ha He — Xaa HÈ) ry” 


Substitution of (B.1) for © = 0, and of the foregoing relations in eqs. (4.2) and 
_ (4.3) leads directly to the virial coefficients of light scattering of eqs. (4.12) and (4.13). 
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The effect of constant and acoustic frequency variable flexural stresses on the dielec- 
tric properties of BaTiO, polycrystals was investigated. The variations of the dielectric 
permittivity resulting from bending the polycrystalline sample were measured in their 
dependence on the strain, polarizing field strength and temperature. Within the ferroelec- 
tric region, flexural stress was found to raise the permittivity, independently of the direc- 
tion of bending. An inversion of the sign of the permittivity variation resulting from dis- 
tention was observed to occur above the Curie point. The relative variations of the permit- 
tivity as measured in the steady state some minutes after application of flexural stress are 
by one order of magnitude smaller than those resulting from homogeneous on e-dimensio- 
nal pressure. 


Introduction 


Investigation of the dependence of the dielectric properties of polycrystalline 
BaTiO, on mechanical stress resulting from the application of one-dimensional pres- 
sure has been continuing for some years. The earliest papers are due to Takagi, Sawa- 
guchi and Akioka (1948), who observed a decrease in the dielectric permittivity and 
a shift in the Curie point towards higher temperatures when pressures of the order 
of 250—750 at were applied. They also found the Curie-Weiss temperature to fall 
as one-dimensional pressure augmented. 

The effect of one-dimensional pressure on the dielectric permittivity of polycrystal- 
line BaTiO, throughout the pressure range of 100—800 at was also investigated by 
Sinyakov and Izhak (1950). They, too, found that the permittivity as measured in 
the direction of compression decreases and the Curie point shifts towards higher 
temperatures. Later papers by Izhak (1956, 1957) showed the relative variation 
of the permittivity due to one-dimensional compression to depend on the measuring 
field strength and on the temperature, and to amount to 6 x 10-4 at~? at room temper- 
ature and a field of I kV/cm. As the pressure varies, the time required for the attain- 
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ment of a steady value of the permittivity is not constant, but depends on the measuring 
field strength, the temperature, and the kind of admixtures present in the barium 
titanate solid solutions. 

A similar investigation of the effect of one-dimensional pressure on the dielectric 
permittivity of BaTiO, was carried out by Lurié and Myedovoy (1956, 1957). Contrary 
to Izhak’s results, they found that the permittivity measured in the direction of 
compression increases below the Curie temperature and decreases above the latter. 
Moreover, they found the sign of the variation of the permittivity to be independent 
of that of the pressure. Their measurements had been carried out at rather low pres- 
sures of the order of 10 at. 

X-ray investigation of structural changes in polycrystalline BaTiO, as arising 
from one-dimensional pressure is due to Ismailzade (1957), who found that, throughout 
the range of 700—1000 at and at room temperature, the volume of the elementary 
BaTiO, cell remains practically constant, whereas only the lattice tetragonality ratio 
c/a increases. Hence, the shift of the Curie point towards higher temperatures, as 
resulting from the application of one-dimensional pressure, is due to induced tetrago- 
nality of the elementary cell. 

Theoretical considerations by Kholodenko and Shirobokov (1951) prove one- 
-dimensional pressure to produce a change in the direction of the spontaneous polari- 
zation in ferroelectrics. In ferroelectrics exhibiting positive volume electrostriction, 
such as BaTiO,, one-dimensional pressure causes the vector of spontaneous polariza- 
tion to assume the direction parallel to that of distention and perpendicular to that 
of compression. The change in the direction of spontaneous polarization as produced 
by one-dimensional pressure brings about a change in the dielectric permittivity both 
in the direction of compression and in the one perpendicular thereto; this is corrobo- 
rated by the results due to Shirane and Sato (1951) and to Izhak (1956, 1957). 

It was felt that investigation -of the variations of the dielectric properties as produ- 
ced by mechanical stress arising in the process of bending of polycrystalline BaTiO, 
samples could contribute to elucidate the subject. The present paper deals with the 
effect, hitherto not investigated, of such stress upon the dielectric properties of poly- 
crystalline BaTiO}. 


§ 1. Preparation of Ferroelectric Material 


l The ferroelectric samples were synthesized from polycrystalline BaTiO, by the 

well known methods of Trzebiatowski, Wojciechowska and Damm (1952), and others, 
using barium carbonate BaCO, and titanium dioxide TiO,. Chemically pure materials 
were used. The compounds were mixed thoroughly in the stoechiometric ratio and 
sintered for 3 hours at 1300°C. The BaTiO, thus obtained was ground to a powder 
and compressed under a pressure of some 2—4 t/cm? to form plates of 2 x 10 x 35 mm. 
Subsequent to compression, the material was again sintered for 2 hours at 1200°C 
and for 3 hours at 1350°C. The BaTiO, samples were then polished down to a thickness 
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of 0,5 mm, and the plane surfaces were provided with electrodes by laying on silver 
paste and heating for 15 minutes at about 800°C. 

The material thus prepared was cut into plates of 0.5 x 3 x 30 mm, which 
were stuck together pairwise to form two-layer rods. The interior electrode of the rod 
(the one between the layers) was continued outwards with a thin sheet of aluminium. 


§ 2. Measuring Device 


The variations in the dielectric properties of polycrystalline BaTiO, arising from 
flexural stress were measured with a resonance method, the block diagram of which 
is shown in Fig. 1. The 275 kHz frequency generator was very weakly induction- 
-coupled to the receiver resonance circuit containing the ferroelectric sample C, the 


D.C. 
Amplifier 


Detector 


HF. 
Generator 


AZ 


Fig. 1. Resonance method for measuring the dielectric permittivity. 


block condenser C, for the permanent polarization of the sample, and the measuring 
condenser Cp. The voltage required to polarize the samples was applied from-a DC- 
source over a safety load of R = 5 MQ. The measuring voltage from the receiver 
circuit was fed over the detector to the DC-amplifier and thence to the horizontal 
plates of the oscillograph. The variations of the dielectric permittivity as due to 
flexural stress were measured on the right and left slope of the resonance curve by the 
method proposed by Jezewski (1928). At a capacity of 4000 pF in the receiver circuit, 
the sensitivity 4 C/C of the resonance method employed amounted to 5 x 10~*. 

The polycrystalline ferroelectric samples in the form of two-layer rods were set 
in a state of flexural vibration of acoustic frequency with the device described in 


another paper by the present author (1960). 


§ 3. Experimental Results 


d the effect of flexural stress of acoustic frequencies on the 


Investigation covere 
as well as the effect of time-indepen- 


dielectric properties of polycrystalline BaTiOg, atte 
dent stress produced by bending the sample once and keeping it bent. 
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Flexural stress was found to raise the dielectric permittivity of BaTiO, both 
in the case of compression and in that of distention. In consequence, the variations 
of the dielectric permittivity caused by the flexural vibrations in the acoustic range 
occur with a frequency that is the double of that of the vibration. The variations of the 
dielectric permittivity e in polycrystalline BaTiO, at flexural vibrations of 400 Hz 
are shown in the oscillogram of Fig. 2. The upper part of the oscillogram shows the 


Fig. 2. Permittivity variation frequency (upper oscillogram) versus frequency of flexural vibrations 
(lower oscillogram). 


variations of ¢ associated with bending in one layer of the rod, whilst the lower one 
shows the variations in the piezoelectric voltage which is given rise to in the linear 
piezoelectric effect that results from bending the other, polarized layer of the poly- 
crystalline rod. Similar variations of the dielectric permittivity, albeit resulting from 
an AC electric field, were observed by Kaczmarek and Kryezkowski (1957) in polyerys- 
talline BaTiO}. 

The increase in € of polycrystalline BaTiO, caused by flexural stress depends 


on the magnitude of the mean strain S, as given by the following formula: 


d 


ps (3.1) 


TT 
here, d denotes the thickness of a single layer of the rod, l — its length, and s — the 
deflection. The total increase of the dielectric permittivity consists of two contributions: 
the one, A €, is constant and probably results from delay effects; it brings about a dis- 
placement of the working point on the slope of the resonance curve of the measuring 
circuit that does not vary during the vibrations; the other, A ¢,, is a variable component 
whose frequency is twice that of the flexural vibrations. The variations of the permittiv- 
ity are shown in the oscillogram of Fig. 3. The horizontal line therein shows the initial 
value of e in BaTiO, previous to application of periodical flexural stress. The depend- 


ence of the relative variation 4e/e of the permittivity on the strain Ss is illustrated 
by the curves in Fig. 4. Curves 1, 2 and 3 yield, respectively, the relative variation 


The Effect of Flexural Stresses on the Dielectric Properties 135 


of the constant component A ¢,/e, that of the variable component A ¢,/e, and that of 


the total mean permittivity A ¢/e as versus the strain S}. The values Je/e were compu- 
ted numerically for the various strains from the relation: 


Acme eee) Ag. 
Se + — 


(3.2) 
E E E 

An external polarizing field F affects the variations of the permittivity as resulting 
from bending a polycrystalline BaTiO, rod. At first, an increase in the polarizing 
field intensity brings about a decrease in the variation of the permittivity caused 
by flexural stress. At a value of the field strength of the order of 2—4 kV/em and 


strain of S} = 2:x 1074, the variation in permittivity becomes dependent on the kind 


Fig. 3. Variable component Ae, and constant component Ae, of permittivity variations as resulting 
from flexural vibrations. 
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Fig. 4. Relative variation of constant component Ae, |e (graph 1), of variable component 4£,/£ (graph 2) 


and of total mean permittivity Ae/e (graph 3) versus the strain S,. 
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of flexural stress involved; the dependence is illustrated by the curves of Fig. 5. Above 
this range of fields, for flexural stress involving compression of the sample in a direc- 
tion perpendicular to that of the polarizing field applied, the permittivity variations 
are almost independent of the field intensity. For flexural stress causing distention 
of the sample perpendicularly to the polarizing field direction, the permittivity varia- 
tion decreases as the latter’s strength mcreases, and approaches 0 for fields of 16—20 


ey Sy, 


Fig. 5. Relative variation of permittivity as resulting from compression (dashes) and from distention 
(continuous), versus the polarizing field strength E. 


kV/cm, according to the polycrystalline material used. The curve (dashes) in Fig. 4 
gives A e/e versus the polarizing field intensity for a sample compressed by flexural 
stress. The continuous curve shows the corresponding dependence for a sample 
subjected to distention. The oscillograms of Fig. 6a, b, c illustrate the permittivity 
variations as arising from sine flexural vibrations in polycrystalline rods at polarizing 
field strengths of 0.8 and 16 kV/cm, respectively, applied to the sample for a time of 
about 80—100 minutes. 

` When investigating the effect of constant flexural stress, applied for several mi- 
“nutes, on the dielectric properties of polycrystalline BaTiO,, delayed phenomena typical 
of the process leading to the attainment of the steady state value of the permittivity 
in polycrystalline BaTiO, were observed to appear. Shirane ant Sato (1951), Izhak 
(1956) and other authors observed similar delayed phenomena arising in polycrystalline 
BaTiO, from one-dimensional homogeneous mechanical stress. These effects as caused 
by an external polarizing field were originally investigated in polycrystalline BaTiO, 
by Partington, Planer and Boswell (1949), Piekara and Pajak (1953), Pajak, Kacunieck 
and Jastrzębski (1956) and others. 

If a polycrystalline BaTiO, rod is bent once in a way to bring about distention 
` or compression of one of either ae of the rod, the permittivity € rises to a maximum 
value e,, that depends on the flexural stress applied, and then decreases slowly during 
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3—5 minutes until it attains a constant value exceeding the initial value of €, as meas- 
ured previous to bending the rod. Following the removal of flexural stress, the permittiv- 
ity again increases to a maximum value somewhat lesser than e,, and subsequently 


AAAA 


GRE LALA ELD LE a a 2 


€ 
Fig. 6. Oscillograms illustrating permittivity variations produced by bending, as versus the polarizing 
field strength, a) E = 0 kV/cm; b)' E = 8 kV/cm; c) E = 16 kV/cm. 


decreases during several minutes to its initial value before þending. Both compreanidn 
ess lead to the same values of the relative variation of 
e variations of the permittivity following bending and 
he oscillograms of Figs. 9 and 10. 


and distention by flexural str 
the dielectric permittiyity. Th : 
subsequent removal of the strain are shown mm t 
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Curves J and 2 (continuous) of Fig. 7 illustrate the dependence of the relative maximum 
variation of the permittivity and that of the relative stationary variation of the permit- 
tivity on the magnitude of flexural stress T, respectively, as given by the following 
formula: 


a 3 d 


Ty = 4, 5, = 4 a (3.3) 
here, cı denotes Young’s modulus. Samples compressed at a pressure of p = 2.3 
and 4 t/cm? previous to the final sintering were used in these measurements; the experi- 
mental values of the relative variation of the permittivity versus the flexural stress 


o (p=2t/cm?) 
o (p=3t/cm?) 
A (p=4t/em?) 


— I, (kg/m?) 


Fig. 7. Graphs of the relative variations of the permittivity versus one-dimensional stress (dashes) and 
versus flexural stress (continuous), for different values of the pressure applied in preparing the polycrys- 
talline samples. 


as shown in Fig. 7 (small circles for samples compressed at 2 t/em?, squares for 3 t/cm?, 
and triangles for a pressure of 4 t/cm?) are such as to eliminate the pressure applied 
in processing the sample as a possible factor influencing the dependence under consid- 
eration, within the range of 2 to 4 t/em®. The dependence itself is linear up to stress 
of T, = 300 kg/cm?, a figure not exceeded in the present investigation. 

From the materials used for preparing the polycrystalline samples subjected to 
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flexural stress, others were made whose permittivity was measured versus the one- 
-dimensional homogeneous stress resulting from compression parallel to the measuring 
field direction. The curves in Fig. 7 (dashes) illustrate the results obtained in collabora- 
tion with J. Klimowski and J. Pietrzak for samples that had been compressed at 2.3 
and 4 t/cm?. 

From Fig. 7, the variations in the permittivity due to one-dimensional stress 
arising in the process of compression are seen to depend on the pressure p applied 
for compressing the sample previous to final sintering. An increase in p as from 2 to 4 
t/cm? produces a rise in the variation of the permittivity due to one-dimensional stress. 
The relative variation of the stationary value of the permittivity due to flexural stress 
is by one order of magnitude smaller than the analogical variation as resulting from 
one-dimensional compression. 


Table I 
Results of measurements of density ọ, Young’s modulus ¢,,, dielectric permittivity € and relative varia- 
tions of the permittivity as arising from unit stress T} when samples manufactured by compression at 
2, 3 and 4 t/cm? are subjected to bending and one-dimensional pressure. 


= = aay: 
maximum variation steady value 
pressure dielectric from bending for one-dimen- 
ae j Young’s ittivi- steady sional compres- 
applied in density ERA iie a a C E nest o 
processing $ 
frequency bending T,= 200 kg/em? 
c E 7 
{ : A E a a eee 1 Ae, 
e WAT eee AT, e€ AT; e AT, 
[g/em®] | [dyn/cm°] [kg/em2]-1| [kg/em?]-+ | [kg/em*]~1)__ [kg/em*]™ 
jose cae 
4.92 USE ANG A PEP 1.6 x 10-4) 1.8x 10-* | 3x 10-5 1G Er 
5.05 7.6 X 1031| 1540 TOA 10 "180% ais 1075 20X J0 
5.17 8.0 x 1011| 1620 1.6 X10 418 x107 | 3 x 10-5 22. %.10+ 
The results of measurements of the dielectric permittivity €, the density 9, Young’s 


modulus c}, and the relative variation of the 
and one-dimensional stress due to compression, 


in processing the sample previous to sintering, 


In order to provide a comparative review of the resu 
and those of other authors relating to the effect of one 
dielectric properties of BaTiOg, all values of the relative va 
resulting from unit stress, hither 
in Table II. The relative variations of the permittiv 
directions parallel and perpendicular to that of the measuring 


linearity of the dependence 


1 Ae 


A 
a AT, and FI amp 


1 


, respectively. With respect to the non 


permittivity as arising from flexural 
as dependent on the pressure p applied 
are assembled in Table I. 

lts of the present investigation 
-dimensional pressure on the 
riation of the permittivity 
to known to have been published, are assembled 
ity arising from compression in the 


field have been denoted 
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HENT Se a 


of Ae/e on one-dimensional pressure, all figures have been referred to stress of ke 

= 200 kg/cm? on the basis of the graphs in the various papers, to the exception of the 

results of Lurié and Myedovoy (1956), who applied weak stress of the order of 10 kg/cm?. 

The various values of ia: and pede in Table II are seen to differ both in magni- 
e AT, e ATi 

tude and in sign. From the results for Ae/e versus the one-dimensional pressure of graph 

7 and of a non-published paper by J. Klimowski and J. Pietrzak, it follows that the 


experimental variation of Ae/e depends strongly on the pressure applied in compress- 


ing the samples previous to their final sintering. The effect of the other technological 


Table II 


Review of results hitherto obtained relating to the effect of one-dimensional pressure on the dielectric 
properties of polycrystalline BaTiO,. 


pressure 

applied in 1 Ae i Ae 

oe e AT, eA Aa 
P 

[t/cm?] [kg/cm?]~* [kg/em?]~> 

= d: loa —0.4 x 10-4 Shirane, Sato (1951) 
= <22 X 10 +3.0 x 10-4 Izhak (1956, 1957) 
— +2.5 x 10-4 — Lurié, Myedovoy (1956) 
— +0.8 x 10-4 +5.0 x 10-4 Ksendzov, Rotenberg (1959) 
2 —1.5 x 10-4 ee Pietrzak, Klimowski (not published) 
4 -+5.0 x 10-4 - — Pietrzak, Klimowski (not published) 
2 +1.5 x 10-4 = Krajewski, Klimowski, Pietrzak 
3 +2.0 x 10-4 — Krajewski, Klimowski, Pietrzak 
4 +2.2 x 10-4 = Krajewski, Klimowski, Pietrzak ` 


parameters, such as the time and temperature of sintering, has not been investigated 
hitherto. Some idea of the effect of the sintering temperature on the permittivity 
Ae 
e AT, 
at the end of Table IT obtained with samples compressed at the same pressures of 2 and 
4 t/cm? but sintered at different temperatures. Papers hitherto published fail to give 
a precise account of the conditions in which the polycrystalline material investigated 


had been manufactured. That is why no reliable comparative study of the various 
results can be made. 


variations due to one-dimensional pressure can be derived from the figures for — 


The variations of the permittivity arising from flexural stress in polycrystalline 
BaTiO, depend on the temperature of the sample. The temperature dependence of the 
relative maximum variation of the permittivity was investigated, and the experimental 
results obtained are shown in the graphs of Fi ig. 8. As the temperature rises Ae/e 
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at first decreases, falling to its lowest value at about 50°C—60°C and then increases, 
to attain its maximum value at a temperature by about 10 degrees centrigrade lower 
than that of the Curie point, after which it decreases rather steeply. The relative varia- 
tion of the permittivity in polycrystalline BaTiO, as arising from distention of the 
sample due to flexural stress is found to change its sign above the Curie point (dashed 
line in Fig. 8). No such effect was observed to arise from compression due to flexural 


AE 4n3 £ 

£ 10 ; 
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30 3000 
20 2500 

AE 3 
\ S~ (compression) 
10 é 000 
Ae 

(6) 1500 
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Fig. 8. Graph of £, and of the relative permittivity variations obtained by flexural distention (continuous) 
and compression (dashes) of polycrystalline BaTiO, at different temperatures. 


stress. For comparison, the same graph brings the values of € as measured in the same 
temperature range before flexural stress had been applied. The variations in-g, as 
dependent on the temperature, observed when the sample was bent once, are illustrat- 
ed by the oscillograms of Figs. 9 and 10. Hence, the effect consisting in the delayed 
attainment of its steady value by the permittivity is seen to occur. only at temperatures 
below the Curie point, and thus should probably be attributed to the presence. of 
domain structure. Above the Curie point, € attains its steady value corresponding 
to that of the flexural stress immediately. es . . = 

-The present investigation failed to reveal a shift in the Curie point arising from 
flexural stress. That due to one-dimensional compression as found by yatious authors 
resulted from stresses largely exceeding anything to be achieved n bending the sample. 
Although the method of investigation adopted here admits of car cel rasp 
ment of the permittivity variations due to compression and distention o t c sample, 
thus presenting a higher degree of universality than those used in similar ee enone 
it brings with it the disadvantage of yielding rather low mechanical stress. The highes 
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Fig. 9. Oscillograms of the permittivity variations resulting from bending (compression) and removal 


af flexure. a) 60°C; b) 140°C. 


é b 
Fig. 10. Qscillograms of the permittivity variations resulting from bending (distention) and removal 


of flexure. a) 60°C; b) 140°C. 


obtained successfully throughout the present investigation did not exceed 300 kg/em?; 
higher stresses lead to cracking of the samples. 
The author wishes to express his thanks to Professor Dr A. Piekara for directing 
the present investigation and for his valuable suggestions and advice. 
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The piezoelectric effect given rise to in polycrystalline barium titanate by flexural 
stress was investigated by the static and quasi-static methods. Samples im the form of 
two-layer rods of 1 X 3 X 30 mm were used. Flexural stress in non-polarized BaTiO; 
rods was found to induce piezoelectric charge on the surfaces perpendicular to the plane 
of bending. The effect thereon of strain and temperature was measured. Moreover, the 
present investigation covered the piezoelectric effect in polarized polycrystalline BaTiO, 
rods and its dependence on the strain, the time of polarization and the polarizing field 
strength. A mechanism accounting for the induction of piezoelectric charge is proposed. 


Introduciion 


In recent years the piezoelectric properties of polycrystalline BaTiO, have been 
the subject of investigation by various authors. Papers were published by Rzhanov 
(1946), Roberts (1947), Mason (1947), (1948), Roy (1950), (1956), Mesnard and 
Eyraud (1955), Marutake and Ikeda (1957), and others. The imvestigations dealt 
chiefly with longitudinal, transversal and radial resonance vibrations and on the effect 
thereupon: of the polarizing field strength and the temperature. The measurements 
were in most cases carried out by dynamic methods. The elastic and piezoelectric 
constants of polycrystalline BaTiO, were computed from the experimental results. 

Hitherto but scarce information is available on the piezoelectric effect exhibited 
by barium titanate rods when subjected to bending. The effect was first investigated 
by Bauer (1948), who used a quasi-static method. It was one of the aims of the present 
paper to investigate this piezoelectric effect by a static and a quasi-static method, and 
to explain the underlying mechanism. 


§ 1. Method of Investigation of Piezoelectric Properties 


Polycrystalline ferroelectric samples were prepared in the shape of two-layer 
rods of 1 x 3 x 30 mm. The piezoelectric charge induced on the surfaces of the sam- 


ple perpendicular to the plane of bending by lasting flexural stresses applied but once. 
(743) 
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was measured with a heterostatically connected Wulf electrometer. In the case of 
induced flexural vibrations of a frequency of up to 2000 Hz, the piezoelectric charge 
was determined from the piezoelectric voltage, as measured with a valve voltmeter 
of an input resistance of about 3 x 106 Q,' and from the dielectric permittivity, as 
measured by the bridge method at a frequency of 1 kHz and a measuring field strength 
of 10 V/cm. The piezoelectric voltage was measured separately for either layer of the 
polycrystalline rod. 

The polycrystalline ferroelectric samples in the form of two-layer rods were set 
in a state of flexural vibration of acoustic frequency ranging up to 2000 Hz with the 
device shown in Fig. 1. A vibrating coil stiffly connected with a concentrically cor- 


membrane 


holder 


Fig. 1. Device for producing flexural vibrations in polycrystalline rods. 


rugated fibre membrane was placed within the field of the slit of a loud-speaker magnet. 
The membrane together with the coil was fixed to the bulk of the magnet with two 
brass rings. The surface of the upper ring supported a holder of organic glass provided 
with electrodes and serving to fix one end of the polycrystalline rod. Its remaining, 
free end was stuck to a thin aluminium rod that was stiffly attached to the loud- 
-speaker membrane. The deflection of the polycrystalline rod was measured to + 5u 
with a microscope having a micrometer scale. 


§ 2. The Piezoelectric Properties of Non-polarized Barium Titanate 


The piezoelectric effect exhibited by crystals belonging to 20 crystallographic 
classes consists in electric charge being induced on the crystal faces as a result of 
. externally applied mechanical stress (Cady 1946; Mason 1950). A general criterium 
for the appertenance of a crystal to one of the piezoelectric classes is given by the absence 
of a centre of symmetry and the presence of polar axes of symmetry. The piezoelectric 
effect also appears in ferroelectric polycrystals after the latter have been polarized 
in a strong external electric field, thus producing non-zero electric polarization. 
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The piezoelectric effect has hitherto not been observed to be induced by mechani- 
cal stress in non-polarized polycrystalline BaTiO,. In the absence of an external electric 
field, the resultant polarization P in polycrystalline BaTiO, is zero because of the ab- 
sence of orientational ordering in the distribution of the crystallographic axes of the 
various ferroelectric domains. This is the reason why non-polarized polycrystalline 
BaTiO, fails to exhibit anisotropy, notwithstanding the fact that the monocrystals 
of which it consists are anisotropic. 

The present investigation proved that the mechanical stresses accompanying flex- 
ural strain of anon-polarized polycrystalline BaTiO, plate give rise to electric charges + q 
induced on the surfaces of the plate perpendicular to the plane of bending. In the absence 
of an external electric field, a positive charge, + q, is induced on the surface that under- 
goes distention by flexural stress, whereas a negative one, — q, is induced on the sur- 
face undergoing compression. If the direction of bending is reversed, so will be the 
sign of the charges induced. Fig. 2 shows the distribution of charge as it appears 
on the surfaces of a non-polarized polycrystalline BaTiO, plate subjected to bending. 


ee 


=O) 


Fig. 2. Distribution of electric charge on surfaces of non-polarized BaTiO, sample subjected to bending. 


Measurements by the present author proved that the piezoelectric voltage U, be- 
tween the electrodes of the polycrystalline sample, resulting from the presence of the 
charges induced by the flexural stresses, depends on the mean value of the strain 


S, as given by the following formula (see Appendix, eq. X): 


EET, 
S (2.1) 


~ 


here, d denotes the thickness of a single layer of the rod, J — its length, and s — the 


deflection. The voltage U, was found to be a linear function of Sy throughout the range 
of strains involved, which were of the order of 4 x 1074. The dependence is shown 
in Fig. 3. À 
If the piezoelectric voltage Up and the dielectric permittivity £ of polycrystalline 
BaTiO, are known, the resulting polarization P giving rise to the electric charges 
on the surfaces of the sample perpendicular to the plane of bending can be computed, 
Assuming the direction of the bending force to coincide with that of the z-axis of 


a Cartesian system of reference and the strain S, to occur lengthwise and in the direc- 


tion of the x-axis, the resultant polarization P, in the direction of the z-axis is given by: 


SED o e2 
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wherein E, is the electric field strength arising from the direct piezoelectric effect. 
For a strain of S, = 1074 at 20°C the resultant polarization amounts to P; = 0.8 x 1074 
uCb/em?. 


Yim) 


ea 


Fig. 3. Piezoelectric voltage U, versus strain Sh for a non-polarized polycrystalline BaTiO, sample. 


In order to give a more general description of the effect observed, use can be 
made of the equation accounting for the direct piezoelectric effect in crystals. The po- 
larization P, is then obtained as follows: 


P; = ds, T; (2.3) 


wherein 7, denotes the mechanical stress in the direction of the x-axis, and ds, — the 
piezoelectric modulus of the material under investigation. Substituting in eq. (2.3) 
the value of eq. (2.2), we have: 


zy E; = dy T} = dy yy S (2.4) 
Hence: 
E Ag 
3 = — Cy, dg, = hay (2.5) 
Sy e 


with hg, denoting the piezoelectric constant of the material. Using the value of Young’s 
modulus c}, = 0.8 x 101? dyn/cm? as computed from dynamic measurements, and 
the experimental value of the permittivity e = 1600, the numerical value of the piezo- 
electric modulus dg, and that of the piezoelectric constant hg, were computed. At room 
temperature these values amount to: 


ds, = 3.2 x 107° CGS and = Ag, = 20 CGS. 
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If the strain S} is kept constant, the value of the piezoelectric voltage Uy is found 
to depend on the temperature. Fig. 4 shows the dependence as measured for a strain 
of Sy = 2.5 x 1074. At first the piezoelectric voltage U increases with the temperature, 
attaining its maximum value at about 50°C. As the temperature continues to rise, 
the voltage Up decreases; at first, this occurs rather steeply, so that the tangent at 
an arbitrary point of this part of the curve intersects the temperature axis at a point 
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Fig. 4. Piezoelectric voltage Ug versus the temperature, for a non-polarized polycrystalline BaTiO, 
sample. 


corresponding to the Curie point temperature; the fall in voltage then becomes less 
and less steep, and it is only at about 170°C that the voltage disappears altogether. 
When the temperature was made to decrease, a similar dependence of the voltage Uo 
on the temperature was observed, the maximum being slightly displaced towards 
lower temperatures. The displacement is probably due to thermal hysteresis of the 
dielectric permittivity as observed in polycrystalline BaTiO, by Piekara and Pajak 
(1952) and in single crystals by Kanzig and Maikoff (1951). 

The piezoelectric voltage Up, especially at higher temperatures, cannot be meas- 
ured with exactitude because of the rather considerable electric conductivity o of BaTiO}. 
The latter was measured in BaTiO, single crystals by Busch and Merz (1948) and in 
polycrystalline material by Trzebiatowski and Pigon (1952) and by Pigon (1954), 
who proved o to rise with the temperature according to the formula: 

— Ea 
o = 0g € $T ulated) 


acteristic of the kind of material used in the meas- 


wherein oo denotes a constant char , 
urements, E, — the a tivation energy, k — Boltzmann’s constant, and T — the Kelvin 


temperature. Within the temperature range of 20°C to 200°C, the conductivity of 
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polycrystalline BaTiO, ø = 1076 + 1074 Q1 cm™!t. The rise in ø towards higher 
temperatures leads to a decrease in the experimental values of the piezoelectric voltage 
corresponding to given values of the strain Si rendering their exact measurement 
impossible. 


§ 3. Piezoelectric Properties of Polarized Barium Titanate 


The dependence of the piezoelectric voltage Up on the strain S, in polycrystalline 
two layer rods polarized by, applying an external electric field, and the effect on the 
value of U, of the polarizing field strength F and of the time during which the polar- 
izing field had been applied, were investigated. The piezoelectric voltage Uy was mea- 
sured separately for either layer of the polycrystalline rod. 

The piezoelectric voltage U, induced when a polycrystalline plate is bent once 
and kept so, decays exponentially over 10—40 seconds. The curves in Fig. 5 show the 


0 10 20 Soano O A E0 70 N, 


Fig. 5. Piezoelectric voltage U, versus the time, for a sample bent once and kept so. 


time dependence of the voltage U, for flexure of various amplitudes. The decay of 
the piezoelectric voltage is due to the rather high electric conductivity of polycrystal- 
line BaTiO}. 

OER by the quasi-static methods at a vibration of 400 Hz of the rod 
showed that the piezoelectric voltage U, depends on the value of the strain Sı. This 
dependence i is illustrated by the curve in Fig. 6. It is linear only for strains of about 
2x 10~*. The voltage Up was not found to depend on the frequency of flexural vibra- 
tions throughout the Tange of 60 to 2000 Hz. 
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Fig. 6. Piezoelectric voltage Up versus the value of the strain ise 
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The value of the piezoelectric voltage U, is, moreover, dependent on the field 
strength E of the electric field applied for polarizing the sample, and on the time 
during which the latter underwent polarization. The dependence of the piezoelectric 
voltage Up on the time of polarization, for a field of E = 10 kV/cm, is shown in Fig. 7. 
At first, as the polarizing field is applied during 10—20 minutes, the voltage rises 


t(min) 


Fig. 7. Preadélecthig voltage U versus the time of polarization, at E = 10 kV/cm. 
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rapidly; later on, it tends to saturation. If the polarizing field E is applied for more than 
70 minutes, no further time-dependent variations in Uù are observed. 

The dependence of U, on the polarizing field strength E, for long-time polariza- 
tion as determined from the measurements discussed above, is illustrated by Fig. 8. 


E (kV/cm) 


Fig. 8. Piezoelectric voltage U, versus the polarizing field strength E. 


As the polarizing field strength increases, the dependence remains practically linear 
up to fields of 6 kV/cm; for still stronger fields, there is a tendency to saturation, 
which occurs at 20 kV/cm and long time polarization. The foregoing results were 
obtained at 20°C. 


The values of the piezoelectric voltage U, the strain S}, Young’s modulus cz; 
and the dielectric permittivity € as obtained by measurements, were used for comput- 
ing the piezoelectric constants d3,, 251, h31 and eg, of polycrystalline BaTiO, polarized 
to saturation in a field of E = 20 kV/cm. The numerical values, in CGS units, are 
assembled in Table I. For the sake of comparison, numerical values of the piezoelectric 


Table I 


Numerical values of the dielectric permittivity £, Young’s modulus c}, and piezoelectric constants 
as obtained from measurements by the quasi-static method in bending polarized polycrystalline BaTiO, 
rods in the course of the present investigation, and by Bauer (1948). All figures are in CGS units. 


1600 | 0.8 x 1012 Present author - 


2000 | 0.8 x 1012 Bauer 


The Effect of Flexural Stresses on the Piezoelectric Properties 751 


constants as computed from results due to Bauer (1948) and obtained by him quasi- 
-statically when subjecting polycrystalline BaTiO, rods to bending are also given. For 
all values of the piezoelectric constants in Table I, there is agreement as to the order 
of magnitude. The divergences in the numerical values are probably due to different 
conditions of preparation of the polycrystalline material. 

It is seen that the piezoelectric moduli dą; of polarized and non-polarized BaTiO, 
differ by two orders of magnitude. Since the piezoelectric modulus dą is defined as 
the variation in resultant polarization P, due to unit stress, this variation in samples 
polarized to saturation is seen to exceed that in non-polarized samples by two orders 
of magnitude. 


§ 4. Mechanism of Piezoelectric Polarization in Bent Polycrystalline Barium 
Titanate Plates 


The mechanism underlying the piezoelectric effect resulting from flexural stress 
in non-polarized polycrystalline BaTiO, has not been fully elucidated. The mechanical 
stress that results when the rod is bent causes one of its layers to become distended 
with respect to the neutral plane of symmetry n (the plane containing the centre of 
the rod and perpendicular to the plane of flexure), whilst the other layer undergoes 
compression (Fig. 9 a). The strain is the greater as the distance from the plane of sym- 
metry increases, and hence a gradient of strain arises parallel to the curvature radius R. 
Considering maximum values S; and S; of the strain for two elementary fibres of 
the bent rod of length l situated at distances AR, and AR, from the neutral plane 
of symmetry n, we have (see Appendix, eqs. I and IX): 


aT TAR OTAR: 
Si = = = Ge wet O S ans (4.1) 
The difference of these expressions yields the strain gradient in the form: 
A Sy e 3 S 5 (4.2) 


AR ork 

with AS, = S; — S; , and AR = AR, — AR,. The strain gradient existing in the 
bent rod, and the inhomogeneous mechanical stresses related thereto, may give rise to: 

1) variation of the orientation of the ferroelectric domains, so that a direction , 
set by the flexural stresses now becomes the prevailing one, or 

2) variation in the shape of the elementary cell and, hence, of the-domain, or else 

3) ordering, by the flexural stresses, of the distribution of defects present in 
the crystal. 

Assuming that the inhomogeneous mechanical stress resulting from flexure when 
a polycrystalline rod is subjected to bending involves a change in the shape of the 
elementary cell, the following model of this change together with a mechanism of 
the induction of the electric charge, which is positive, +g, on the distended and nega- 
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tive, — q, on the compressed surface, can be put forward. The elementary cell of 
BaTiO, can be schematically represented as in Fig. 9b (Shirane, Jona and Pepinsky, 
1955). The same cell, according to the model proposed here, assumes. the shapes shown 
in Fig. 9 c, d, the former relating to the distended layer and the latter to the one sub- 
jected to compression. This change in the shape of the elementary cell should cause 


Fig. 9. Model showing change in shape of elementary BaTiO, cell, as arising from flexural stresses. 


displacement (squeezing) of the titanium ion Titt from its position of equilibrium 
in a direction coinciding with that of flexure; such stress-induced polarization will 
generally give rise to electric charge on the surface layer of the ferroelectric, of posi- 
tive sign on the convex side and negative on the concave one. The charge in the sur- 
face layers will produce electrostatic induction on the electrodes, and the free charges 
= will come to be observed as piezoelectric charge. 

From investigations by Fienkel (1926), Wagner and Schottky (1930) and various 
other authors, a crystal lattice at thermal equilibrium is known always to present 
a number of ions or atoms situated at interstitial points (Frenkel defects) and a number 
of non-invested node points (Schottky defects). Other kinds of defects are also possible, 


é 
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thus e.g. various types of dislocation. The number of defects present is given by an 
exponential function of the activation energy and the temperature (Mott and Gurney, 
1950). It is also known that, in an external electric field, the defects exhibit a well- 
-defined mobility accounting, in some cases, for the electric conductivity of the crystal. 
The assumption that the electric charge observed is due to defects or dislocations and 
to their motion in the inhomogeneous mechanical stress field, and that the direction 
of their motion and the magnitude of displacement is given by the direction of the 
stress gradient, would lead to predict an ordered distribution of defects within the 
crystal. 

The present state of the investigation as yet admits of no decision as to which 
of the interpretations proposed is the correct one in accounting for the mechanism 
of the effect observed. Investigation on polycrystalline BaTiO, and on single crystals 
is proceeding at this Laboratory with the aim of elucidating the mechanism in question. 

The distribution of the piezoelectric charge induced by flexural stress on the 
surfaces of a polarized sample perpendicular to the plane of bending is dependent 
on the residual polarization P, from the external polarizing field F and on the direc- 


tion of the strain Se Such charge results from the variation in residual polarization P_, 
caused by the flexural stresses. The static method used for investigating the piezo- 
„electric effect yields the exact relationships between the sign of the piezoelectric 
charge induced on a surface perpendicular to the plane of bending by variation of 
the residual polarization P,, on the one hand, and the original direction of the re- 
sidual polarization and that of strain, on the other. The results obtained by the static 
method are assembled in Fig. 10. 

According to Mason (1948) and Smolensky (1951), BaTiO, belongs to the group 
of ferroelectrics exhibiting positive volume electrostriction. By theoretical considera- 
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10. Distribution of piezoelectric charge as dependent on direction of residual polarization P, and 


Fig. i 
direction of the strain S,. 
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tions due to Kholodenko and Shirobokov (1951) and from measurements carried out 
by Mason (1955), resultant polarization in a ferroelectric exhibiting the positive 
effect of volume electrostriction assumes the direction parallel to that of the distention 
brought about by the external mechanical stresses, and perpendicular to that of the 
compression. Accounting for the way the direction of resultant polarization in BaTiO, 
depends on those of distention and compression, and accounting for the distribution 
of piezoelectric charge in its dependence on the direction of the residual polarization 


P, and that of the strain Si as observed in the paper cited, the following mechanism 
of the induction of these charges in polarized BaTiO, may be proposed. Polarization 
gives rise to surface charges of density ọ = +P, in the surface layer of the ferroelectric; 
these induce charge in the electrodes. If external stress T, = 0, the bound charges 
induced in the electrodes fully compensate the polarization charges, whereas the free 
charges are compensated as a result of the electric conductivity of the sample. If the 
residual polarization P, is directed towards the plane of junction (denoted by the 
dashed line in Fig. 10) as shown in Fig. 10a, then bending associated with distention 
of the layer induces positive charge on the convex surface and negative charge on 
the concave side. The mechanical stress arising in the process of bending, and causing 
distention of the upper layer of the rod, reduces the polarization P, by the amount 
AP,, arraying the latter in the direction of distention. Due to the decrease in the 
polarization P_, the positive charge of the electrode, which hitherto compensated the 
negative charge of polarization P,, now becomes a free charge and is available for 
observation in the electrometer as piezoelectric charge. If the direction of bending 
is reversed so as to compress the layer of the rod in the direction perpendicular to 
that of the polarization P,, the latter will be increased by. the amount AP,. This, 
in turn, gives rise to induction on the electrodes, as a result of which the free charge, 
negative in the case under consideration, comes to be observed as piezoelectric charge. 
On the other hand, if the residual polarization is directed from the plane of junction, 
as shown in Fig. 10 b, distention diminishes it by the amount AP_, as a result of which 
charges hitherto compensating the polarization are set free on the surfaces of the layer 
subjected to flexure. Since the direction of the polarization P, is now reversed with 
respect to the foregoing case, the signs of the charges observed will be reversed, too. 

Thus, from the experimental distribution of piezoelectric charge, distention of 
a polycrystalline sample perpendicular to the residual polarization P, is seen always 
to diminish this polarization, whereas compression of the sample perpendicular to 
the direction of P, always results in an increase of the polarization P,. 


Conclusions 


Mechanical stress appearing when a non-polarized polycrystalline BaTiO, plate - 
is subjected to flexure causes electric charge to be induced on its surfaces perpendic- 
ular to the plane of flexure. The effect may be due to variation in the orientation 
of the ferroelectric domains occurring in the direction favoured by flexural stress, 
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or to variation in the shape of the elementary cell and domain, or to ordering of the 
distribution of defects in the crystal lattice. 

The distribution of piezoelectric charge induced by flexural stress on the surfaces 
of the polarized BaTiO, sample perpendicular to the planè of bending depends on 
the direction of the residual polarization resulting from the polarizing field, and on 
the direction of strain. Distention of the polycrystalline sample in the direction perpen- 
dicular to that of the residual polarization always produces a decrease whereas compres- 
sion in the same direction always leads to an increase in the latter. 


Appendix 
Method of Measuring Flexural Stresses 


Bending of a rod, of constant cross section along its length J, causes distention 
of one of the layers and compression of the other one with respect to the neutral 
surface n through the centre of the rod and perpendicular to the plane of bending. 
The strain S, arising from distention of e.g. the upper layer of the rod in Fig. 11 is 


: ar83 
Fig. 11. Section through bent rod, oag derivation of eq. S, = vert s. 
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Or 


given by the relative variation of the latter’s length A///. With the notation of Fig. 11, 
we have: 


os I 


wherein R is the radius of curvature of the neutral surface n at point b, and d — the 
thickness of the layer subjected to bending. Since the length / of the rod on the neu- 
tral surface remains unchanged in the process of bending, the mean strain of a layer 
of thickness d is given by: 


II 
The curvature at any point of the neutral surface, in the coordinates x, z of Fig. 11, 
is given by: 


Reg eae celts Ill 


wherein M(x) denotes the bending moment acting upon the ie section of the rod 
at a distance 1—x from the point of application of the bending force F, and aS the = 


stiffness modulus of the rod for bending. To all purposes, at small doflezions, J ns ae E 


Hence: 
1 dz MOL ee) Windi » 
o dt fis <i giriga a 
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The angle of bending p of the neutral surface is given ni 
f WF ao gG $ (he AS Npp > ai 


i 
i 
) 
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For small values of the angle of bending ¢g, the bending curve n of Fig. 11 can be 
replaced by a circular arc of radius R. For an arc of length /, the angle ø is now 
given by: 


Comparing eqs. VII and VIII, the curvature of the neutral surface of a rod of length Z 
is obtained in the form: i 
1 3 os 


= 2 


Te ie i 
Substituting the value of 1/R of eq. IX in eq. II, the mean strain S, is obtained as a func- 
tion of the deflection: 


5 3 d 


ary ee i 
According to Hooke’s law, the mechanical stresses arising in the rod from bending 
are given by: 


=, 3 d A 


wherein c,, is Young’s modulus in the direction of the x-axis. Eq. XI, which was de- 
rived according to Southwell’s monograph (1941), served here for computing the 
flexural stresses arising in polycrystalline barium titanate. 

The author wishes to express his great indebtedness to Professor Dr A. Piekara 
for his valuable advice and encouragement throughout the present investigation and 
for his numerous helpful suggestions and discussions. 
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In the present paper an attempt is made to discuss the magnetic superstructures 
from the point of view of spin waves propagation. The classification of superstructures pro- 
posed by the author in earlier papers is here developed more precisely, introducing sym- 
metric and non-symmetric superstructures and dividing the former into regular and non- 
-regular orders. The present paper deals with the regular order only. The discussion com- 
prises mainly antiferromagnetic translational lattices, though the general conclusions hold 
also for ferri- and antiferrimagnetic superstructures. 

As the most characteristic property of regular superstructures, the author assumes 
the presence of alternating parallel ferromagnetic planes. Their existence allows to pre- 
dict the phenomenon of superstructural anisotropy of spin waves dispersion. The presumable 
behaviour of spin waves in particular superstructures of cubic lattices is discussed in 
detail. The theoretical proofs of these hypotheses are announced. 

Moreover, a simple geometrical interpretation of Luttinger’s method of description 
of superstructures is given. 

I. Introduction 


In the present paper, the classification and some general properties of magnetic 
superstructures are considered from the point of view of the theory of spin waves. 
The subject covers all magnetic translational lattices consisting of two kinds of spins: 
A and B, each kind forming a ferromagnetic sublattice. These sublattices are labelled 
A and B respectively. The numbers of the lattice points in A and in B are assumed 
to be equal; consequently, both interpenetrating sublattices will be identical. Considera- 
tions deal chiefly with antiferromagnetic lattices, but the classification and general 
conclusions hold also for ferri- and antiferrimagnetic ones. The detailed discussion 
concerns those superstructures, which will here be termed “regular”. 

In the subsequent sections, certain theorems concerning the various categories 
of superstructures introduced in Chapter II will be given. It should be pointed out 
that some of these theorems present the character of empirical rules. It is to be hoped 
that in the future they may be derived theoretically. 

(759) 


760 H. Cofta 


IT. Classification of superstructures 


In previous papers (Cofta 1958 and 1959, Szczeniowski and Cofta 1959) a classi- 
fication of superstructures which appeared convenient for studying spin waves in 
ferri- and antiferrimagnetics was proposed. This classification will now be developed 
more precisely. 


1) Symmetric and non-symmetric superstructures 


The first necessary step consists in dividing all superstructures considered into 
two classes, which we shall define as follows. A superstructure will be said to be sym- 
metric if each lattice point is a center of symmetry. In other words, we are dealing 
in this case with such orders of points A and B for which the following rule holds: 


pn em WWE ony 
Fig. 1. Example of non-symmetric Fig. 2. Example of non-symmetric superstructure in fcc 
superstructure in bcc lattice lattice 


ifthe lattice point r belongs to A, then —r also belongs to A, and analogically for the 
sublattice B. This is clearly valid for any origin of the coordinates system (fixed at 
an arbitrary lattice point). 

All superstructures that are not symmetric will be termed non-symmetric. From 
the above definition it follows immediately that, in the case of a non-symmetric super- 
structure, for each lattice-point there exists at least one such relative position r, belong- 
ing to A, for which —r belongs to B and vice versa. The appearence of such non-sym- 
metric linear chains ... AABBAABB... of nearest or second neighbours distingu- 
ishes essentially the non-symmetric superstructures from symmetric ones with respect 
to the propagation of spin waves in antiferromagnetics, since the ordering ... A A BB 
AABB... of the spins probably constitutes a kind of “obstacle” for spin waves 
(Cofta 1959). Two examples of unsymmetric superstructures are shown in F igs. 
1 and 2. A detailed analysis of unsymmetric orders will be given in a subsequent 
paper. Let us but-state that the existence of unsymmetric linear chains implies directly 
the impossibility of a translational character of the sublattices A and B, which there- 
fore must be lattices with basis. 
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2) Regular and non-regular superstructures 


In the case of symmetric superstructures two possibilities exist: 
a) both sublattices A and B are translational lattices; 
b) the sublattices A and B are not translational lattices (i.e. they are lattices with 
a basis). i 

For the first of these two categories, let us retain the name “regular superstruc- 
tures”, introduced previously (Cofta 1958 and 1959) in a not entirely exact manner. 
We shall examine the regular orders in the ensuing Chapters. 

The second category of symmetric superstructures will be termed “non-regular”. 
In this case, calculation of spin wave problems is more complicated than for regular 
superstructures. Indeed, the equations for the spin wave amplitudes must be written 
separately for each translation sublattice. Therefore, if the basis of sublattice A (or B) 


Fig. 3. Example of non-regular superstructure in fcc lattice 


consists of n points, we have to solve 2n simultaneous equations, leading in general 
to 2n branches of the energy spectrum. Instead, in the case of regular superstructures 
we have n — 1 and therefore always only two branches are obtained. 

The. class of non-regular superstructures does not comprise a large number of 
examples!, as compared with the other two classes (i.e. regular and unsymmetric ones). 
A typical example of non-regular superstructure is that of the hypothetic “order 
of the second kind” of the fcc lattice, considered often in the literature (see e.g. Van 
Vleck 1951 and Roth 1958). This superstructure is shown in Fig. 3. 


III. Characteristic property of regular superstructures 


1. Ferromagnetic lattice planes 


It is a well known fact that in some cases of superstructures, ferromagnetic 
sheets occur. Sets of parallel lattice planes, each consisting exclusively of A points 
or exclusively of B points, exist in the case of non-symmetric as well as in that of 
Je S sete tae hee, eS 


1 See the note at the last page of the present paper 
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regular superstructures (their existence in non-regular superstructures does not seem 
possible). Among the manifold types of sets of parallel ferromagnetic sheets, the most 
interesting is the case of alternatingly arranged planes: ... A B AB ... In this case, 
either of both sublattices consists of identical equidistant lattice planes, and conse- 
quently is a translational lattice. This proves that alternating ferromagnetic sheets 
can occur only in the case of regular orders. We believe that they must occur in each 
regular superstructure. Hence, we will consider the presence of alternating ferromagnet- 
ic sheets as the characteristic property of regular superstructures. 


2. Three types of regular orders 


The foregoing property enables us to distinguish three different types of regular 
superstructures. As a criterion of such a classification we shall use the number of 
different sets of alternating parallel ferromagnetic planes, taking into account only 
those with the lowest Miller indices (i.e. 0 or 1). These three types are the following: 

a) Monoplanar type. Only one set of parallel ferromagnetic sheets exists. (As 
to examples, see Ch. V). 

p) Biplanar type. Two differently directed, equivalent sets of parallel ferromag- 
netic sheets exist. (For examples, see Ch. V). | 

y) Multiplanar type. More than two equivalent sets of alternating ferromagnetic 

_sheets, each having a different direction, exist. This type is realised only in the ortho- 


Fig. 4. Natural (multiplanar) superstructure of the simple orthorhombic lattice 


rhombic, tetragonal and cubic Bravais systems. It is readily seen that the symmetry 
degree determines strictly the number of sets in the case of “multiplanar” type: 
there can be only 4 different equivalent sets, namely the 4 sets of planes (111) (Fig. 4). 
However, it should be mentioned that in the case of body centered lattices there exist 
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also 3 sets (100), (010), (001) with lower Miller indices (see Fig. 5), whereas in the case 
of face centered lattices the multiplanar type cannot be realised. 

As can be seen from Fig. 4, a. 5, the “multiplanar” type of regular superstructure 
is identical with the well-known order in which every nearest neighbour of any given 
lattice point belongs to the other sublattice. As proposed previously (Cofta 1958), 
we prefer to use for this case the term “natural order”. 

In the subsequent sections it will be shown that all the natural orders must pos- 
sess the same physical properties; moreover, they will be seen to present the least 
interesting cases of regular superstructures. 


Fig. 5. Natural (multiplanar) superstructure of the body centered orthorhombic lattice 


Our observations on the existence and role of sets of alternating ferromagnetic 
sheets provide a simple way of marking the individual regular superstructures. This 
consists in giving, in addition to the lattice symbol, the Miller indices of the unique 
(for uniplanar) or two (for biplanar orders) sets of ferromagnetic sheets. Natural 


superstructures require no additional marking. 


IV. Luttinger’s description of superstructures 


The results of Luttinger’s (1951) considerations relating to the energy of Ising’s 
model provide a very simple method of description of individual superstructures. 
It is sufficient to find such a vector V, appropriate for the given superstructure, for 


which: 
wr is an even number, if r belongs to A, 
wr is an odd number, if, r belongs to B, 
or vice versa. Then, for each point of A, the number 
= (ee ee (1) 


or vice versa (Ziman 1952).. Let us term V 


is equal to 1, and for each point of B, to —1, 
ctor of the appropriate reciprocal lattice. 


the Luttinger vector. Clearly, v is a lattice ve 
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Luttinger has calculated the vectors » for the following 5 superstructures (in 


the notation introduced above): sc-(110)-(110), fcc-(001), fec-(111) and for both 
natural superstructures of cubic lattices (sc and bcc). The Luttinger vectors for other 
superstructures — if existent — could be found by guessing. When Luttinger vectors 
are used for calculations by the spin-wave method (Ziman 1952), no information is 
available on the range of validity of the formulas obtained, since nothing is as yet known 
on the existence of Luttinger vectors for the various superstructures. It seems that 
e.g. Ziman’s dispersion law for antiferromagnetics is often held to be valid for few 
orders only (see e.g. Elliot and Lowde 1955). 

The situation is entirely clarified on considering the connection between the 
Luttinger vectors and the alternating ferromagnetic sheets. The connection is a very 
simple one: the Luttinger vector is a vector perpendicular to the ferromagnetic 
planes, so that their equation is of the form 


wr =n (2) 


n being an integer. Choosing the origin e.g. at a point of A, we obtain, for n even, 
the equations of the A planes, and, for n odd, those of the B planes. In other words, 
the number 


Creole (3) 


changes its sign on transition from any ferromagnetic plane to the neighbouring one. 
Now, it is obvious that for the biplanar superstructures we can find two equivalent 
Luttinger vectors having different directions, and for the natural orders there exist 
as many as 4 different Luttinger vectors which can be used in the description of such 
superstructures. 

The geometrical interpretation of Luttinger vectors proposed here yields immedi- 
ately the important conclusion that Luttinger’s method of the description of super- 
structures is applicable to all regular orders, and to these only. Consequently, the 
scope of validity of Ziman’s (1952) dispersion law for spin waves in antiferromagnetics 
can be strictly determined: this law is valid exactly for regular superstructures of 
translational (magnetic) lattices (see also Cofta 1959, page 222). 

Instead of Miller indices we may, of course, employ the Luttinger vectors for 
marking indivdual superstructures. 


V. Superstructural anisotropy of spin waves dispersion 


1. Isotropic and anisotropic superstuctures 


The division of regular superstructures into three types has not only a formal 
meaning but also corresponds to physical differences in the propagation of spin — 
waves. Let us first consider the monoplanar orders. The very characteristic arrangement 
of spins undoubtedly privileges the » direction perpendicular to the ferromagnetic 
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sheets. This suggests that the same direction will be the preferred one with regard 
to the propagation of spin waves. For the same reason, it should be expected that 
in the biplanar superstructures the direction » = », x V, along the line of intersection 
of two nonparallel ferromagnetic planes will be the privileged one for the spin waves. 

By the foregoing, the monoplanar as well as the the biplanar orders can be consid- 
ered to be anisotropic superstructures. In such superstructures, the spin waves 
should obey an anisotropic dispersion law, independently of the presence of the struc- 
tural anisotropy. According to an earlier proposal of ours (Szczeniowski and Cofta 
1959), the possible effect predicted here will be termed “superstructural anisotropy 
of spin wave dispersion”. A more detailed discussion is inserted in Chapter VI. 

In the natural superstructures the situation is quite different. For reasons of lattice 
symemtry all four v directions of (111) ferromagnetic sheets are equivalent and the 
same holds for the six directions of intersection of ferromagnetic planes. Therefore 
no direction is preferred. This conclusion expresses the well-known fact that for natu- 
ral superstructures the isotropic dispersion law holds (naturally, if there is no structural 
anisotropy). Hence, the expected behaviour of spin waves in mono- and biplanar 
orders should be essentially distinct from their behaviour in natural orders. 

It is to be noted, however, that the conclusions of the present section are valid 
only when the exchange integral between two atomic spins depends solely on their 
distance i.e. when the direct as well as indirect exchange interactions are isotropic. 
If in positions A and B there are different ions, our suggestions must be slightly correc- 
ted. If, however, the chemical lattice is so highly complicated that the indirect ex- 
change coupling is different between different spins of the same sphere of neighbours 
or, moreover, when semicovalent exchange interactions between coplanar orbitals 
exist, then our conclusions must be radically revised. 


2. Effects of second interactions 


When the second interactions i.e. interactions between next-nearest neighbours 
are present, then an important role is played by the substructures of the given lattice. 
Let us recall that each Bravais lattice of a rectangular system can be divided into 
2 or 4 interpenetrating lattices of the same system. These identical component lattices 
are called substructures or submotifs (Roth 1958). In this way: 


1° a simple lattice is composed of 2 face centered ones; 
2° a body centered lattice is composed of 2 simple ones; 
3° a face centered lattice is composed of 4 simple ones. 


The role of substructures is clear with respect to the well-known fact that each sub- 
structure consists solely of points which are second neighbours in the whole lattice. 
A natural suggestion results from this fact. Namely, in the dispersion law the term 
due to second interactions should show a character corresponding to the type of spin 
order of the given substruetures. Thus, when substructures have mono- or biplanar 
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order, an anisotropic second interactions term may be expected. Instead, when the 
substructures are ferromagnetic or present the natural order, an isotropic second 
interactions term may be anticipated (when not taking into account the purely struc- 


tural effects). 


VI. Regular superstructures in cubic antiferromagnetic lattices 


We shall now proceed to review the particular regular superstructures which 
can appear in cubic lattices. The cubic superstructures are the most convenient 
examples, not only due to their simplicity but also because of the absence for them 
of structural anisotropy. Although here we confine our considerations to antiferro- 
magnetics, the conclusions will be valid probably also for ferrimagnetic translational 
lattices. , 

Let the lattice constant in the structures considered be a. All conclusions of this 
Chapter hold exactly only if the condition of isotropy of exchange interactions is 
satisfied. 


1. Review of cubic regular superstructures 


a) Natural superstructures 


First of all we shall settle the matter of natural superstructures. Such orders can 
appear in sc and bcc lattices only; in the fcc lattice their existence is impossible, since 
in this structure the nearest neighbours of a given lattice point are also nearest neigh- 
bours to one another. In contradiction to the opinions met with sometimes, it should 
be stated that the natural superstructure is no typical order in antiferromagnetics 
(as well as in ferrimagnetics with antiparallel spins). It can appear only when negative 
i.e. antiferromagnetic coupling between nearest neighbours predominates. On the 
other hand, in most substances with antiparallel spins, strong negative second inter- 
actions seem to exist. Only three cases of the natural superstructure in cubic lattices 
have been revealed hitherto. These are the antiferromagnetic perovskites LaCrO, 
(Koehler and Wellan 1957), LaFeO, (l.c.) and CaMnO, (Wollan and Koehler 1955) 
with simple cubic magnetic lattices. All these compounds owe their natural super- 
structure to the fact that here each two nearest magnetic ions are colinearly separated 
by an intervening oxygen ion. Let us mention that, outside the cubic system, also 
only three cases of natural order have been observed, namely in the antiferromagnetics 
MnF, FeF, and CoF (Erickson 1953) with tetragonal magnetic structure. 


b) Superstructure sc-(001) 


_ This monoplanar superstructure, determined unambiguously by our symbol, 
is shown in Fig. 6. It can be described by aid of the Luttinger vector a 


v=a*e o 


= 
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Only one case of such superstructure seems to have been discovered hitherto, namely, 
that of the antiferromagnetic perovskite LaMnO, (Wollan and Koehler 1955). Unfor- 
tunately, this case does not satisfy our restricting condition. Probably we have here 
coplanar hybridized orbitals in magnetic ions, coupled ferromagnetically by semico- 


Fig. 6. Superstructure sc-(001) 


valent exchange within the planes (001) and antiferromagnetically between two neigh- 
bouring (001) planes (Shull and Wollan 1956, page 217). In the hypothetical lattice 
in which our condition would be satisfied, the superstructure sc-(001) could pre- 
sumably exist when the negative second interactions predominate over the positive 
(i.e. ferromagnetic) nearest interactions. ` 


c) Superstructure sc-(110)-(110) 


This biplanar superstructure is shown in Fig. 7. Here we have the choice between 


two Luttinger vectors: 
vı =a? (a + b) 
V, = a`? (a — b) (5) 


each of them being suitable for describing this superstructure. The order under 
consideration has been observed in (0.2 La — 0.8 Ca), MnO (Wollan and Koehler 
1955), but it is not sure that our condition of isotropic interactions is fulfilled.in the 
case of this compound. If this condition is satisfied, is seems that the superstructure 
se-(110)-(110) can exist when the negative second interactions predominate over the 


negative nearest ones. 
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Fig. 7. Superstructure sc-(110)-(110) 


d) Superstructure bcc-(110) 


This monoplanar superstructure has not been considered in the literature. It 
may be seen in Fig. 8. The corresponding Luttinger vector is 


v=a*(a+ Db) (6) 
No realisation of such order in nature is known at present. The possibility of its ex- 


istence is rather questionable. A more detailed discussion on the problem of its 
existence will be given in a subsequent paper. 


e) Superstructure fcc-(001) 


This superstructure, shown in Fig. 9, may be described by the Luttinger vector 


vV = 2a? e (7) 
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It is known in the literature as the “order of the first kind” of fcc lattice (Van Vleck 
1951, Smart 1952). This superstructure is realised in the antiferromagnetic substance 
MnTe, (Corliss et al. 1958). It seems very probable that in this case the nearest inter- 
actions are negative and the second ones most probably positive, though other rela- 
tions may also be possible. 


f) Superstructure fcc-(111) 


This superstructure, shown in Fig. 10, has often been discussed in the literature. 
It may be described by the Luttinger vector 


v = (a + ò + e)a? (8) 


Smart (1952) and Labhart (1953) use the term: “order of the second kind” for this 
superstructure of fcc lattice. The superstructure fc-(111) has been found to occur in 


Fig. 10. Superstructure fcc-(111) 


the well-known case of MnO, in other monoxides of transition metals, such as FeO, 
CoO, and NiO. (Shull et al. 1951) and also in a-MnS (Corliss et al. 1956). It seems 
certain that such a superstructure is a consequence of considerable preponderance 
of negative second interactions. 


2. Anisotropic effects to be expected in cubic regular superstructures 


The possible superstructural anisotropic effects, of which a general discussion 
was given in Chapter V, can be predicted in detail for each of the individual cubic 
superstructures described in the previous section. The particular hypotheses listed 
below will be labelled P}, Pp,... and so on. These hold only if the interactions in all 
directions are the same at the same distance. 


a) Natural superstructures 


Here the situation is well-known: the propagation of spin waves is governed by 
an isotropic dispersion law. 
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b) Superstructure sc-(001) 


Let us note that, within each ferromagnetic plane, the spins present in this case 

a quadratic arrangement implying rotational symmetry of the superstructural aniso- 

tropy. Moreover, the fcc substructures occuring here have the (001) order; 

therefore, the second interactions privilege the same direction [001] as the nearest 
ones. Hence, the following expectations on the superstructural dispersion anisotropy 
may be formulated: 

P,: The nearest interactions term in the dispersion law should present anisotropy 
with the preferred direction [001]. This anisotropy should present rotational 
symmetry about [001] i.e. should be uniaxial. 

P,: The second interactions term should present anisotropy with the same preferred 
direction [001] as the nearest interactions term. Likewise, this anisotropy will 
be presumably uniaxial (see P, below). 


c) Superstructure sc-(110)-(110) 


The privileged direction is here also the [001] one. In planes (001) the spins are 
arranged antiferromagnetically in quadratic layers and show the natural (two-dimen- 
sional) order. Each of two fcc substructures presents the (001) order, the same as in 
the previous case. Hence, we can predict the following behaviour of spin waves: 
Ps: The nearest interactions term should present the uniaxial anisotropy with preferred 

direction [001] f 
P,: The second interactions term should present the uniaxial anisotropy in the same 

direction [001] 


d) Superstructure bec-[110] 


This superstructure undoubtedly privileges the direction [110]. The spin in each 
ferromagnetic plane (110) forms a centered rectangular array, so that the [110] 
axis is only a two-fold one. Therefore no rotational symmetry can be expected. Each 


of the two sc substructures shows the biplanar order (110)-(110) discussed above _ 


as case c). Hence, the following hypotheses may be made: 

Ps: The nearest interactions term should be anisotropic with the most preferred 
direction [110]. 

P: The second interactions term should show uniaxial anisotropy preferring the 


direction [001]. 


e) Superstructure fcc-(001) 


The privileged [001] direction is a four-fold axis. Each of the four sc substruc- 
tures is a ferromagnetic one, and therefore privileges no direction whatever. Hence, our 
hypotheses in this case will be the following: 
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P,: The nearest interactions term in the dispersion law should present uniaxial 

anisotropy with the preferred direction [001] 

Py: The second interactions term: should be isotropic. 
f) Superstructure fec-(111) 

The favoured [111] direction represents a three-fold axis, suggesting rotational 
symmetry of anisotropy. Each of the four sc substructures exhibits a natural order. 
Hence, the following hypotheses result: 

P,: The nearest interactions term should show uniaxial anisotropy in the direction 

[111]. 

Pio: The second interactions term should be isotropic. 

The hypotheses P,, P, and P, have already been confirmed by our previous 
approximate calculations (Cofta 1957). 

The hypotheses as to spin waves behaviour presented here will be proved in 
a subsequent paper. 
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Notes added in proof. 

Note I. 

Our statement as to the small number of different non-regular superstructures, expressed in Ch. II, 
§ 2 of the present paper, is not true. Recently, the author has discovered 5 other cases of such type. 
Note II. ' 

Quite recently the author has read the paper by Gersch and Koehler, published in the J. Phys. Chem. 
Solids, 5, 180 (1958). As consistent with the Ising model, accurately the same 7 cubic superstructures 
presented here; has been obtained by Gersch and Koehler. 
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ERRATA 


„On the Theory of Magnetic Anisotropy of Ferromagnetic Crystals“ — L. 
lewski (Acta Phys. Polon. 19, 59 (1960)). 
B (Rm) on page 63 should read: 


B (Rp) = AY (Ry) S—2P (Ry) (S—3 (G By) Rig + 20ST Sip 
+20 (Rp) [S?— 35-2 (LS) (Si Rig) Riri + (S — 4) (Sy- Rn) 


The line 7 on page 65 (“We confine our attention to the case of IKI 
long spin waves“) should be printed at the end of page 71. . 


page line incorrect 


| Rea ec) ISPS 


63 (14 oe where 
Ta a APR S—OP (Rig) =: 
AE oo 


W 


